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Linear  analysis  and  direct  numerical  simulation  (DNS)  are 
used  to  study  the  interaction  of  a  shock  wave  with  a  turbulent  flow.  Effects  central 
to  shock  wave/boundary  layer  interaction  are  identified  and  explained.  The  quantitaciv 
importance  of  the  turbulence  anisotropy  and  the  shock  oblique  angle  is  demonstrated. 
Drop  in  Reynolds  shear  stress  across  a  normal  shock  is  shown.  The  implication  for 
Reynolds  stress  models  is  discussed.  The  role  of  dilatatlonal  turbulent  fluctuations 
on  turlrulence  amplification  is  studied  in  detail.  An  explanation  is  suggested  for 
lower  amplification  reported  in  wind-tunnel  studies.  Linear  analysis  &  DNS  are  used 
to  demonstrate  the  considerable  influence  of  entropic  fluctuations  on  turbulence 
evolution  across  a  shock.  Positive  correlation  between  u'  and  T’  suppress  amplificatio 
while  negative  correlation  enhances  it.  Invalidity  of  Morkovin's  hypothesis  is 
demonstrated.  A  high  order  scheme  is  developed  to  allow  strong  shock/turbulence 
interaction.  Results  are  presented  for  isotropic  turbulence/normal  shock  interaction. 
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OVERVIEW 


This  report  sumniEirizes  our  study  of  the  interaction  of  a  shock  wave  with  turbulent 
flows  over  the  period  1992-1994.  Two  doctoral  students  were  supported  by  this  program; 
Dr.  Sangsan  Lee  obtained  his  degree  in  1992  (his  thesis  was  pubhshed  as  Report  TF-52; 
a  copy  was  sent  to  Dr.  Sakell)  and  Mr.  Krishnan  Mahesh  is  expected  to  complete  his 
work  in  1994.  The  salient  achievements  of  this  program  are  presented  below.  Details  are 
provided  in  the  attached  enclosures  of  relevant  reports  and  journal  pubh cations. 

A  combination  of  linear  analysis  and  direct  numerical  simulation  was  used  by  Lee, 
Lele  and  Moin^^~^^  to  study  the  interaction  of  isotropic  turbulence  with  a  normal  shock 
wave.  Subsequent  to  graduation,  Dr.  Lee  continued  work  on  this  problem  as  a  post¬ 
doctoral  fellow  at  the  Center  for  Turbulence  Research.  His  principal  accomplishments  axe 
as  follows; 

1.  A  numerical  scheme  that  made  possible  the  first  direct  ntimerical  simulation  of  shock  / 
turbvilence  interaction  was  developed.  The  significant  achievements  include  the  devel¬ 
opment  of  a  method  to  specify  ‘turbulence’  at  the  inflow  boxmdaries  of  computational 
domains  and  the  adaptation  of  a  high  order  shock-capturing  scheme  to  allow  accurate 
computation  of  strong  shocks. 

2.  Linear  analysis  was  conducted  in  parallel  to  the  simulations.  The  emalysis  was  seen  to 
reproduce  most  features  of  the  interaction,  an  eispect  that  suggests  its  incorporation 
into  turbulence  models.  The  analysis  also  offers  strong  support  for  the  decrease  in 
tiurbvJence  length  scale  across  the  shock  wave  -  a  fact  that  points  to  inadequacies  in 
some  of  the  current  experimental  work  on  this  problem. 

3.  The  database  generated  was  \ised  to  examine  a  model  for  the  evolution  of  the  shock 
front  in  a  timbulent  flow  by  Kevlahan, Mahesh  and  Lee^^^  Fxirther  use  to  evaluate 
turbiilence  models  is  anticipated. 

The  project  was  extended  by  Mahesh,  Lele  and  Moin^®“®^  to  study  the  interaction  of 
a  turbulent  shear  flow  with  a  shock  wave.  Both  analysis  and  direct  numerical  simulation 
were  used.  Mr.  Mahesh  is  currently  in  the  final  stages  of  his  dissertation.  His  principal 
achievements  are  s\immarized  below: 

1.  The  method  to  specify  inflow  turbulence  was  considerably  improved  upon.  The  mod¬ 
ified  method  made  possible  the  shear  flow  computation  which  we  believe  is  the  first  of 
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its  kind.  Other  direct  numerical  and  large  eddy  simulations  at  Stanford  and  elsewhere 
have  now  adopted  this  method. 

2.  Linear  analysis  was  used  to  identify  and  explain  the  considerable  importance  of  the 
anisotropy  of  the  incident  turbulence  in  shock/turbulence  interaction.  The  inability 
of  current  models  to  reproduce  these  effects  was  documented. 

3.  Analysis  and  computation  were  used  to  demonstrate  the  considerable  influence  of 
entropy  fluctuations  and  the  nature  of  their  correlation  with  the  velocity  field  in 
shock/turbulence  interaction.  The  results  were  explained  and  a  scaling  developed 
to  allow  modeling  of  these  effects  whose  importance  in  shock  wave  /boundary  layer 
interaction  seems  to  have  gone  unnoticed. 

4.  The  effect  of  compressibility  on  the  evolution  of  compressible  turbulence  across  a 
shock  wave  was  studied  in  detail.  The  results  suggest  a  possible  explanation  for  wind 
tunnel  experiments  on  shock/turbulence  interaction  reporting  different  results  from 
shock  tube  experiments. 

5.  The  shear  flow  computation  is  currently  in  progress  and  reveals  effects  that  were 
hitherto  unreported  in  the  literature.  The  preceding  analyses  explain  the  observed 
behavior.  A  database  is  currently  being  generated  to  allow  evaluation  of  turbulence 
models. 

This  report  is  organized  as  follows.  An  introduction  to  Mr.  Mahesh’s  Ph.D  thesis  is 
provided  in  Part  1.  Parts  2  and  3  deal  with  the  influence  of  the  anisotropy  of  the  turbu¬ 
lence  and  tmbulent  entropy  fluctuations  respectively  on  shock/turbvdence  interaction.  The 
influence  of  compressibility  on  shock/turbulence  interaction  is  described  in  Part  4.  Part 
5  concludes  the  report  with  the  interaction  of  isotropic  turbulence  with  a  strong  shock 
wave. 
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INTRODUCTION 


1.1  Motivation  and  background 

Recent  interest  in  high  speed  civil  transport  aircraft  and  hypersonic  propulsion  sys¬ 
tems  has  resulted  in  a  resurgence  of  research  into  compressible  turbulent  flows.  At  high 
speeds,  compressible  effects  become  very  important  to  the  design  of  both  external  and 
internal  surfaces  such  as  wings,  flns,  control  surfaces  and  nozzle  intakes.  The  compressibil¬ 
ity  of  the  medium  results  in  large  gradients  in  the  mean  thermodynamic  properties,  large 
fluctuations  in  the  thermodynamic  Vciriables  and  high  levels  of  mean  dilatation.  Simple  ex¬ 
tensions  of  incompressible  turbulence  models  appear  unable  to  predict  these  phenomena, 
particularly  those  associated  with  mean  and  fluctuating  dilatation.  Incorporating  these 
effects  of  compressibility  into  turbulence  models  has  motivated  fundamental  research  into 
the  physics  of  these  flows. 

The  presence  of  shock  waves  is  an  important  featme  that  distinguishes  high  speed 
flows  from  low  speed  ones.  Interaction  with  shock  waves  is  known  to  significantly  alter 
the  hydrodynamics  and  acoustics  of  high  speed  shear  layers.  Shock  wave/boundary  layer 
interaction,  a  common  flow  configuration,  can  induce  flow  separation  and  strong  fluctuating 
pressure  loads  on  external  surfaces.  The  presence  of  shock  waves  in  high  speed  unadapted 
jets  is  known  to  significantly  increase  the  emitted  noise.  On  the  constructive  side,  shock 
waves  offer  the  potential  to  increase  the  mixing  in  high  speed  flows.  Current  understanding 
of  the  physics  of  the  interaction  is  hmited. 

This  work  is  a  fimdamental  study  of  the  interaction  of  a  shock  wave  with  a  turbulent 
shear  flow.  We  use  linear  analysis  and  direct  numerical  simulation  to  study  the  interac¬ 
tion  of  a  normal  shock  wave  with  a  shear  flow  having  transversely  uniform  mean  velocity 
gradient  and  Mach  number. 

1.2  Review  of  past  work 

This  review  deals  with  the  hydrodynamics  of  the  interaction  of  shock  waves  with  non¬ 
reacting  turbulent  flows  in  the  non-hypersonic  regime.  Phenomena  such  as  detonation 
waves,  Richtmeyer-Meshkov  instability,  regimes  of  shock-reflection  and  shock-related  noise 
(eg.  ‘screech’)  are  not  discussed.  Experimental  work,  linear  analysis  and  numerical  compu¬ 
tation  are  reviewed  here.  Experimental  work  and  analysis  of  shock  turbulence/interaction 
span  about  five  decades.  Computational  work  on  the  other  hand  is  very  recent. 
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1.2.1  Experiments 

A  large  body  of  experimental  work  exists  on  the  interaction  of  shock  waves  with  a  tur¬ 
bulent  boimdary  layer.  Most  experiments  have  examined  the  flow  in  a  compression  corner 
with  some  experiments  studying  normal  shock/boundary  layer  interaction  and  supersonic 
shear  layers  at  reattachment. 

Green  (1970),  Korkegi  (1971)  and  Hankey  and  Holden  (1975)  provide  reviews  of  earlier 
experimental  work  on  shock  wave/boundary  layer  interaction.  The  AGARDograph  by 
Fernholz  and  Finley  (1981)  provides  a  comprehensive  summary  of  experiments  prior  to 
1981.  More  recent  reviews  are  written  by  Delery  (1985)  and  Marvin  and  Delery  (1987). 
A  lot  of  experimental  data  has  been  gathered  in  the  last  decade.  Current  experiments 
provide  data  on  mean  velocity  and  wall  pressure,  streamwise  Reynolds  stress,  shear  stress 
and  wall  pressiure  fluctuations.  Some  information  on  two  point  correlations,  spectra  and 
PDFs  is  also  available  {eg  Smits  and  Muck,  1989). 

The  flow  regime  covered  by  experiments  on  compression-comer  flows  extends  from 
attached  to  separated.  Amplification  of  the  streamwise  Reynolds  stress  (u^),  structure 
parameter  (mJ/u^)  and  temperature  fluctuations  across  the  shock  wave  was  observed  by 
Debieve,  Gouin  and  Gaviglio  (1982a,  1982b),  Dussauge,  Muck  and  Andreopoulous  (1986), 
Smits  and  Muck  (1989),  and  Sehg,  Andreopoulous,  Muck,  Dussauge  and  Smits  (1989). 
Kuntz,  Amaducci  and  Addy  (1987)  observed  that  the  boundary  layer  after  the  interaction 
showed  an  acceleration  of  the  mean  flow  near  the  wall  as  it  recovered  to  equilibriiim  and 
that  the  mean  streamwise  velocity  profile  downstream  of  the  comer  was  wavy. 

The  unsteadiness  of  the  shock  wave  in  the  separated  regime  was  studied  by  An¬ 
dreopoulous  and  Muck  (1987)  who  analyzed  wall  pressure  signals  and  concluded  that  the 
frequency  of  the  shock  wave  motion  was  controlled  by  the  bursting  frequency  of  the  inci¬ 
dent  boundary  layer.  Subsequent  experiments  by  Erengil  and  Dolling  (1991)  also  show  a 
direct  correlation  between  the  motion  of  the  shock  wave  and  pressure  fluctuations  in  the 
upstream  boundary  layer.  Further  work  aimed  at  correlating  the  motion  of  the  shock  wave 
to  pressure  fluctuations  in  the  bovmdary  layer  is  conducted  by  Erengil  and  Dolling  (1990), 
McClure  (1992)  and  Gramann  and  Dolling  (1992). 

One  of  the  results  of  shock-induced  separation  is  fluctuating  pressure  loads  on  the 
wall.  Dolling  and  coworkers  have  conducted  several  investigations  into  the  wall  pressure 
fluctuations  produced  by  shock/wave  boundary  layer  interaction.  Both  two-dimensional 
and  three  dimensional  flows  are  considered.  An  extensive  review  of  experimental  data  is 
provided  by  Dolling  (1993). 
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An  experiment  by  HayaJcawa  et  al  (1984)  examined  a  reattaching  supersonic  flow. 
Strong  longitudinal  compression  was  found  at  reattachment  resulting  in  a  three-fold  rise 
in  the  peah  level  of  mass  flux  (^).  A  comparable  subsonic  reattachment  experiment  by 
Bradshaw  and  Wong  (1972)  showed  a  two-fold  drop  in  the  level  of  u. 

In  recognition  of  the  complex  nature  of  shock  wave/boimdary  layer  interaction,  recent 
experiments  have  studied  the  interaction  of  a  shock  wave  with  ‘free’  turbulence.  Debieve 
and  Lacharme  (1986)  examined  the  interaction  of  a  shock  wave  with  the  separated  shear 
layer  over  a  compression  ramp.  They  concluded  that  turbulence  levels  were  amplified  and 
Taylor  microscales  increase  through  the  interaction.  An  intermittency  effect  in  turbulent 
statistics  at  the  shock  front  location  was  also  described. 

The  interaction  of  grid  generated  turbulence  with  a  shock  wave  was  studied  by  Keller 
and  Merzkirch  (1990).  Amplification  of  kinetic  energy  was  observed  along  with  an  ap¬ 
parent  increase  in  the  Taylor  microscale  of  density  across  the  shock  wave.  Honkan  and 
Andreopoulous  (1990)  examined  the  response  of  grid  generated  turbulence  to  a  normal 
shock.  Considerable  amplification  of  turbulence  levels  was  noted.  The  amplification  was 
seen  to  depend  on  the  length  scale  and  intensity  of  the  incident  flow.  Another  experimental 
investigation  of  the  interaction  of  a  shock  wave  with  grid  generated  turbulence  was  con¬ 
ducted  by  Jacquin,  Blin  and  Geffroy  (1991)  who  noticed  negligible  ampHfication  of  kinetic 
energy  across  the  shock  wave. 

Detailed  experiments  on  the  interaction  of  a  weak  shock  wave  with  random  fluctuations 
in  density  were  conducted  by  Hesselink  and  Sturtevant  (1988).  Considerable  distortion  of 
the  shock  front  was  observed  by  them  which  they  explained  in  terms  of  focusing  and 
defocusing  of  the  shock  front  in  response  of  density  gradients  in  the  medium. 

As  noted  by  Jacquin,  Cambon  emd  Bhn  (1993),  experiments  on  grid  turbulence/shock 
wave  interaction  are  not  easy  to  conduct  and  certain  points  of  contention  exist.  The  in¬ 
crease  in  length  scale  that  experiments  report  is  both  counter-intuitive  and  in  contradiction 
to  computation  and  linear  analysis.  Lele  (1993)  observes  that  the  inconsistency  is  currently 
unresolved.  Another  feature  is  the  lower  amplification  of  kinetic  energy  reported  by  wind 
tunnel  experiments  as  compared  to  shock  tube  experiments.  Mahesh  ei  al  (1993)  suggest 
that  the  presence  of  acoustic  waves  in  the  wind-tunnel  experiments  is  responsible;  however 
further  data  is  needed  to  resolve  this  issue. 

1.2.2  Linear  analysis 

The  analysis  of  shock/turbulence  interaction  is  based  upon  the  decomposition  of  a 
compressible  flow  field  into  vorticity,  acoustic  and  entropy  modes  as  suggested  by  Kovasz- 
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nay(1953).  The  independence  of  these  modes  for  inviscid,  tmiform  mean  flow  is  used  to 
study  the  interaction  of  each  of  these  modes  with  a  shock  wave.  Both  free  shock  waves  and 
wedge-attached  shocks  have  been  considered  in  the  hterature.  The  analyses  by  different 
workers  differ  in  detail  but  are  conceptually  similar.  The  linearized  Euler  equations  are 
used  to  describe  the  flow  upstream  and  downstream  of  the  shock  wave  which  is  modelled 
as  a  discontinuity.  The  linearized  Rankine-Hugoniot  equations  yield  boimdary  conditions 
at  the  shock  front.  The  resulting  boimdary  value  problem  for  the  flow  downstream  of  the 
shock  wave  is  then  solved  to  yield  the  displacement  of  the  shock  front  and  the  flow  behind 
the  shock  wave. 

The  earliest  studies  seem  to  be  those  by  Blokhintzev  (Landau  and  Lifshitz,  1982), 
Burgers  (1946)  and  Kantrowitz  (1947),  who  examined  the  one- dimensional  interaction  of 
an  acoustic  wave  with  a  normal  shock  wave.  Carrier  (1949)  used  linear  analysis  to  study 
the  stability  of  supersonic  flow  past  a  wedge.  Interest  in  the  conflguration  of  shock  waves  in 
supersonic  channel  flow  prompted  Adams  (1949)  to  study  the  steady  interaction  between  a 
shock  wave  and  an  acoustic  wave  incident  from  upstream.  Linear  analysis  was  used  by  Chu 
(1952)  to  examine  the  interaction  between  a  wedge-generated  shock  wave  and  an  acoustic 
wave  incident  from  downstream.  The  interaction  between  sound  and  a  shock  wave  was 
also  studied  by  Lighthill  (1949). 

Prompted  by  the  problem  of  ‘shock-noise’,  Ribner  (1953)  studied  in  detail  the  inter¬ 
action  between  a  vorticity  wave  and  a  shock  wave.  He  (1954)  subsequently  extended  his 
analysis  to  consider  a  spectrum  of  incident  vorticity  waves  (in  three  dimensions)  and  com¬ 
puted,  for  an  isotropic  incident  spectrum,  detailed  statistics  of  the  downstream  flow  field 
with  emphasis  on  the  generated  noise.  Ribner  appears  to  have  been  the  first  to  predict 
the  interaction  of  isotropic  turbulence  with  a  shock  wave.  He  noted  (1953,1969)  that  weak 
levels  of  incident  turbiilence  produce  high  levels  of  noise  downstream  of  the  shock  and  later 
(1987)  extended  his  analysis  to  predict  one-dimensional  power  spectra  of  the  fluctuations 
behind  the  shock  wave. 

Around  the  same  time  as  Ribner’s  pioneering  study,  Moore  (1954)  performed  detailed 
analysis  of  the  unsteady  interaction  of  an  obliquely  incident  acoustic  wave  with  a  shock 
wave  of  infinite  extent.  He  studied  acoustic  waves  that  were  incident  from  downstream 
as  well  as  upstream  and  also  outlined  extension  of  his  analysis  to  study  the  interaction  of 
a  vorticity  wave  with  a  shock  wave.  Other  workers  that  have  studied  the  interaction  of 
unsteady  disturbances  with  a  shock  wave  include  Kerrebrock  (1956),  Johnson  and  Laporte 
(1958),  Lowson  (1968)  and  McKenzie  and  Westphal  (1968). 

Most  of  these  workers  considered  the  interaction  of  vorticity  and  acoustic  waves  with  a 
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shock  wave.  Chang  (1957)  used  linear  analysis  to  examine  the  interaction  of  a  plane  entropy 
wave  with  an  oblique  shock.  In  addition  to  considering  the  general  case  of  an  infinite 
oblique  shock,  Chang  also  considered  the  case  where  the  oblique  shock  was  produced  by 
a  wedge  and  accounted  for  refiections  from  the  wedge.  Chang’s  theoretical  results  were 
subsequently  used  by  Cuadra  (1968)  to  perform  a  numerical  parametric  study  of  the  infinite 
oblique  shock  problem. 

These  earlier  studies  were  largely  motivated  by  issues  of  shock-stability  and  the  ‘shock- 
noise’  phenomenon.  Recent  interest  in  the  interaction  of  shock  waves  with  turbulent  fiows 
has  led  to  renewed  interest  in  hnear  analysis.  Lee,  Lele  and  Moin  (1992,  1993,  1994) 
repeated  Ribner’s  (1953,1954)  analysis  and  extended  it  to  examine  some  turbulent  statistics 
that  Ribner  had  not  considered.  Mahesh,  Lee,  Lele  and  Moin  (1994a)  extended  Moore’s 
analysis  to  examine  the  interaction  of  cin  isotropic  field  of  acoustic  waves  with  a  shock  wave. 
They  also  outhned  description  of  the  simultaneous  interaction  of  a  field  of  acoustic  and 
vortical  disturbances  with  a  shock  wave.  Mahesh,  Lele  and  Moin  (1994c)  have  examined 
the  simultaneous  interaction  of  vorticity  and  entropy  fluctuations  with  a  shock  wave  of 
infinite  extent.  The  interaction  of  both  single  modes  and  an  isotropic  field  of  fluctuations 
with  the  shock  wave  was  studied.  Interest  in  receptivity  led  Duck,  Lasseigne  and  Hussaini 
(1994)  to  study  the  processing  of  free  stream  disturbances  by  a  wedge- generated  shock 
wave.  The  different  regimes  of  the  pressure  field  behind  the  shock  wave  were  discussed. 

Recently,  Lee,  Lele  and  Moin  (1992)  have  suggested  that  idealization  of  a  shock  wave 
as  a  homogeneous  one-dimensional  compression  and  examination  of  the  temporal  response 
of  turbulence  to  this  compression  might  approximate  the  spatial  evolution  of  turbulence 
through  a  shock  wave.  Linear,  inviscid  equations  or  Rapid  Distortion  Theory  (RDT)  were 
used  to  study  this  problem.  The  results  showed  reasonable  agreement  with  Ribner’s  anal¬ 
ysis  for  shock  waves  whose  mean  Mach  number  was  less  than  2.  Significant  disagreement 
was  seen  for  stronger  shock  waves.  This  conclusion  is  reiterated  by  Jacquin,  Cambon  and 
Blin  (1993)  who  in  addition  compare  the  homogeneous  compression  of  compressible  tur¬ 
bulence  to  Ribner’s  results  and  note  that  the  difference  is  larger  them  observed  by  Lee  et 
al. 

It  is  difficult  to  systematically  outline  the  differences  between  the  two  analyses;  how¬ 
ever  certain  general  comments  can  be  made.  Idealization  of  a  shock  wave  as  a  one¬ 
dimensional  compression  prevents  it  from  distorting  in  response  to  the  incident  fluctu¬ 
ations;  as  a  result  RDT  cannot  represent  the  effect  that  shock  front  distortion  has  on 
the  evolution  of  the  turbulence.  Also  RDT  has  no  knowledge  of  the  Rankine  Hugoniot 
equations;  as  a  result  solenoidal  fluctuations  do  not  produce  all  three  kinds  of  modes  upon 
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compression.  Another  important  difference  between  the  two  analyses  is  the  behavior  of 
acoustic  fluctuations.  The  acoustic  fluctuations  in  the  spatial  problem  can  either  prop¬ 
agate  unattenuated  or  decay  exponentially  behind  the  shock  wave  depending  upon  the 
angle  of  incidence  of  the  disturbance.  RDT  cannot  reproduce  this  behavior.  Mahesh  ei 
al  (1994b)  compare  their  analysis  of  a  shock  wave  interacting  with  acoustic  fluctuations 
to  Jacquin  et  aVs  (1993)  compression  of  compressible  turbulence  and  note  that  the  results 
axe  different  over  the  entire  range  of  Mach  numbers.  It  appears  that  the  homogeneous 
problem  can  approximate  the  evolution  of  solenoidal  fluctuations  across  shock  waves  of 
moderate  strength;  however  it  is  quite  inappropriate  as  the  incident  turbulence  becomes 
increasingly  compressible. 

There  has  been  some  recent  work  on  the  evolution  of  shock  waves  in  an  unsteady 
fleld  of  disturbances.  Lele  (1992)  has  used  Rapid  Distortion  Theory  to  calculate  changes 
in  the  shock- speed  caused  by  the  interaction  of  isotropic  turbulence  with  a  shock  wave. 
Kevlahan,  Mahesh  and  Lee  (1992)  modelled  the  shock  front  as  a  propagating  surface  and 
compared  the  resulting  evolution  to  numerical  simulation.  Good  agreement  was  obtained 
for  the  cases  considered.  A  recent  paper  by  Giddings,  Rusak  and  Cole  (1993)  has  developed 
a  transonic  small  disturbance  model  to  analyze  the  interaction  of  a  weak  shock  wave  with 
random  pertmrbations.  The  model  allows  effects  such  as  nonlinear  steepening,  focusing 
and  the  formation  of  caustics  to  be  represented. 

1.2.3  Numerical  computation 

Numerical  computation  of  shock/ turbulence  interaction  is  a  fairly  recent  undertaking. 
The  flrst  such  computation  appears  to  be  a  detailed  study  of  isotropic  turbulence  inter¬ 
acting  with  a  weak  normal  shock  (Mi  <  1.2)  by  Lee,  Lele  and  Moin  (1992,  1993).  Lee  et 
al  solved  the  three  dimensional  compressible  Navier  Stokes  equations,  using  a  non-miiform 
mesh  to  resolve  the  thickness  of  the  shock  wave.  Subsequently,  they  (1994)  used  a  high 
order  shock  capturing  scheme  to  study  stronger  shock  waves  (Mi  =  2, 3).  Their  results 
corroborated  Ribner’s  (1954)  linear  analysis.  Amplification  of  vorticity,  kinetic  energy  and 
thermodynamic  fluctuations  was  observed  along  with  decrease  in  Taylor  microscale. 

Subsequently,  Hannappel  and  Friedrich  (1994b)  have  used  shock  capturing  to  compute 
compressible  isotropic  turbulence  interacting  with  a  shock  wave  with  Mi  =  2.  The  interac¬ 
tion  of  compressible  fluctuations  (kinetic  energy  equally  distributed  between  solenoidal  and 
dilatational  modes)  was  contrasted  with  that  of  solenoidal  fluctuations.  The  amplification 
of  kinetic  energy  was  seen  to  be  lower  when  the  incident  fluctuations  were  compressible. 
Hannappel  and  Friedrich’s  computation  might  suffer  from  insufficient  resolution  of  the 
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shock  front  (Lee  et  al,  their  results  however  axe  in  qualitative  agreement  with  the 
linear  analysis  of  Mahesh  et  al  (1994b). 

Prior  to  Lee  et  aV s  (1992)  study,  computations  were  largely  restricted  to  the  interaction 
of  a  shock  wave  with  single  waves  of  vorticity,  sound  and  entropy  or  deterministic  flow 
inhomogeneities  such  as  hot  spots  and  vortices.  There  appear  to  be  only  two  detailed 
studies  (in  two  dimensions)  of  the  interaction  of  a  shock  wave  with  a  random  field  of 
disturbances. 

Zang  et  al  (1984)  evaluated  the  linear  analysis  of  McKenzie  and  Westphal  (1968) 
by  comparison  to  their  numerical  solution  of  the  two-dimensional  Euler  equations.  They 
examined  the  effect  of  incident  angle,  shock-strength  and  the  amplitude  of  the  incident 
disturbance  in  the  interaction  of  acoustic  and  vorticity  waves  with  a  shock  wave.  Results 
were  presented  only  in  the  freely  propagating  regime.  The  dependence  on  incident  angle 
was  examined  in  the  interaction  of  disturbances  of  amplitude  0.1%  and  10%  with  a  shock 
wave  of  Mach  number  8.  Good  agreement  was  seen  away  from  the  critical  angle;  divergence 
from  the  linear  analysis  prediction  was  seen  within  about  20°  of  the  critical  angle.  The 
linear  analysis  predictions  were  quite  robust  in  terms  of  dependence  on  shock  strength  and 
disturbance  amplitude.  For  an  incident  angle  of  30°,  the  linear  predictions  were  valid  for 
disturbance  amplitudes  as  high  as  25%  for  acoustic  waves  and  100%  for  vorticity  waves 
and  shock  waves  whose  Mach  number  was  close  to  unity. 

More  recently,  Mahesh,  Lee,  Lele  and  Moin  (1994)  have  computed  the  interaction  of 
a  Mach  1.5  shock  wave  with  an  acoustic  wave  incident  from  upstream.  The  dependence 
upon  the  incident  angle  was  examined  and  the  results  were  compared  to  the  linear  anal¬ 
ysis  of  Moore  (1954).  Good  agreement  was  seen  except  within  25°  of  the  critical  angles. 
This  disagreement  was  explained  by  noting  that  the  region  around  the  critical  angle  cor¬ 
responded  to  transonic  flow  in  transformed  coordinates  and  hence  linear  approximation 
Would  be  inconsistent  in  that  regime. 

Hannappel  and  Friedrich  (1994a)  have  computed  the  interaction  of  a  Mach  8  shock 
wave  with  vorticity  and  acoustic  waves  incident  at  45°  from  upstream.  They  contrasted 
the  evolution  of  these  two  waves.  The  terms  in  the  averaged  equations  for  the  kinetic 
energy  were  also  examined  in  their  study. 

Numerical  calculations  of  the  interaction  of  coherent  flow  disturbances  with  a  shock 
wave  have  also  been  performed.  Hussaini  et  al  (1986)  studied  the  effect  of  upstream  eddy 
motion  and  temperature  inhomogeneity  on  the  enhancement  and  production  of  turbulence. 
Meadows  et  al  (1991)  examined  the  effect  of  vortex  strength  in  a  study  of  two-dimensional 
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shock/ vortex  interaction.  Several  workers  have  computed  the  interaction  of  a  shock  wave 
with  a  density  inhomogeneity  to  examine  the  nature  of  the  vorticity  field  that  is  produced. 
A  list  of  references  may  be  found  in  a  recent  paper  by  Yang,  Kubota  and  Zukoski  (1994). 

Rotman  (1991)  performed  numerical  solution  of  the  two-dimensional  Euler  equations 
to  study  the  response  of  an  isotropic  random  field  of  disturbances  to  a  moving  shock  wave. 
Amplification  of  kinetic  energy  and  decrease  of  length  scale  were  observed  in  his  study. 
Comparison  of  Rotman’s  results  to  linear  analysis  by  Lee  (Private  communication)  showed 
good  agreement. 

Prior  to  their  three  dimensional  computation,  Lee,  Lele  and  Moin  (1991)  solved  the 
two-dimensional  compressible  Navier  Stokes  equations  to  study  the  interaction  of  two- 
dimensional  isotropic  turbulence  with  a  shock  wave.  A  non-uniform  mesh  was  used  to 
resolve  the  thickness  of  the  shock  wave.  Amplification  of  vorticity  and  kinetic  energy  was 
observed  accompanied  by  a  drop  in  Taylor  microscale. 

1.3  Overview 

Our  primary  objective  is  to  study  the  interaction  of  a  shock  wave  with  a  turbulent 
shear  flow.  Both  Unear  analysis  and  direct  numerical  simulation  are  used  for  this  purpose. 
To  our  knowledge,  DNS  of  this  problem  has  not  been  conducted  before.  The  DNS  provides 
a  database  that  can  be  used  to  elucidate  the  flow  physics  as  well  as  to  test  turbulence 
models.  The  Unear  analyses  were  conducted  to  supplement  the  simulations  in  predicting 
the  interaction. 

The  principal  contributions  and  findings  of  this  work  are  listed  below. 

•  Rapid  Distortion  Theory  was  used  to  study  the  response  of  turbulent  shear  flows 
to  rapid  one-dimensional  compression.  Both  normal  and  oblique  compressions  were 
considered. 

•  The  anisotropy  of  shear  flows  (as  determined  by  i?ii(/ci)  and  u^/g^)  and  the  obliquity 
of  compression  were  shown  to  strongly  affect  the  evolution  of  the  Reynolds  stresses 
upon  compression. 

•  Reynolds  stress  modelling  of  the  normal  compression  problem  was  considered  in  the 
RDT  limit.  The  inabiUty  of  popular  pressure  strain  models  to  predict  the  evolution 
of  kinetic  energy  and  the  decrease  in  the  magnitude  of  shear  stress  was  demonstrated. 

•  Linear  analysis  of  the  combined  interaction  of  vorticity  and  entropy  fluctuations  with 
a  shock  wave  was  carried  out.  The  interaction  of  a  single  wave  as  well  as  an  isotropic 
field  of  fluctuations  was  studied. 
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•  Significant  dependence  of  the  downstream  solution  on  the  amplitude  ratio  and  phase 
difference  between  the  vorticity  and  entropy  waves  was  demonstrated. 

•  The  presence  of  negatively  correlated  velocity  and  temperature  fluctuations  was  shown 
to  increase  the  amplification  of  kinetic  energy  and  pressure  fluctuations  across  the 
shock  wave.  Likewise,  positively  correlated  temperature  and  velocity  fluctuations 
resulted  in  suppression  of  amplification. 

•  The  inapplicability  of  the  Strong  Reynolds  Analogy  (SRA)  across  a  shock  wave  was 
demonstrated  and  explained. 

•  The  interaction  of  a  Mach  1.5  shock  wave  with  a  single  acoustic  wave  was  numerically 
computed  and  the  dependence  upon  the  incident  angle  examined.  The  results  were 
compared  to  linear  analysis  and  limits  placed  on  the  range  of  validity  of  linear  analysis. 
The  requirement  of  transonic  small  disturbance  equations  to  model  the  interaction 
around  the  critical  angle  was  demonstrated. 

•  Linear  analysis  was  used  to  study  the  interaction  of  an  isotropic  field  of  acoustic  waves 
with  a  shock  wave.  The  evolution  of  kinetic  energy,  sound  level  and  thermodynamic 
fluctuations  were  examined  in  detail. 

•  The  spatial  evolution  of  kinetic  energy  behind  the  shock  wave  was  examined  and  the 
rapid  non-monotonic  variation  immediately  downstream  of  the  shock  wave  explained. 

•  The  kinetic  energy  of  acoustic  fluctuations  was  shown  to  decrease  across  the  shock  for 
Mach  numbers  varying  from  1.25  to  1.8.  For  Mach  numbers  exceeding  3,  the  kinetic 
energy  was  seen  to  amplify  by  levels  that  far  exceeded  those  foimd  in  the  interaction 
of  vortical  fluctuations  with  a  shock  wave.  These  trends  were  explained. 

•  The  linear  analysis  was  used  to  suggest  acoustic  fluctuations  in  the  upstream  turbu¬ 
lence  as  a  possible  explanation  for  wind-tunnel  experiments  on  grid  turbulence/shock 
wave  interaction  reporting  lower  levels  of  kinetic  energy  amplification. 

•  Comparison  of  the  linear  analysis  to  Jacquin  et  aFs  analysis  of  the  homogeneous 
compression  of  compressible  turbulence  revealed  that  the  homogeneous  problem  was 
quite  inappropriate  as  the  the  incident  turbulence  became  increasingly  compressible. 

•  Direct  numerical  simulation  of  the  interaction  of  a  shock  wave  with  a  turbulent  shear 
flow  was  performed.  The  method  of  specifying  inflow  turbulence  developed  by  Lee, 
Lele  and  Moin  (1992)  was  modified  to  reduce  the  evolution  length  of  the  inflow  dis¬ 
turbances  by  a  factor  of  4-5. 
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This  report  is  organized  as  follows.  Ra,pid  Distortion  Theory  is  used  to  study  the  re¬ 
sponse  of  a  homogeneous  turbulent  shear  flow  to  rapid  compression  in  Chapter  2.  Chapter 
3  discusses  linear  analysis  of  the  combined  interaction  of  vorticity  and  entropy  fluctuations 
with  a  shock  wave.  Details  of  the  direct  numerical  simulation  are  presented  in  Chapter  4. 
The  significant  results  and  suggestions  for  future  work  are  outlined  in  Chapter  5.  Of  the 
three  appendices,  Appendix  A  describes  the  linear  analysis  of  the  interaction  of  a  shock 
wave  with  an  isotropic  field  of  acoustic  waves.  Details  of  the  analysis  of  interaction  of 
vorticity  and  entropy  waves  with  a  shock  wave  are  presented  in  Appendix  B.  Finally,  Ap¬ 
pendix  C  outlines  the  modified  method  for  specification  of  inflow  turbulence  in  numerical 
simulations. 
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PART  TWO 


THE  RAPID  COMPRESSION  OF  A  TURBULENT  SHEAR  FLOW 


2.1  Background 

This  chapter  uses  rapid  distortion  theory  (RDT)  to  examine  the  response  of  turbulent 
shear  flows  to  rapid  one-dimensional  compression.  Answers  are  sought  to  the  following 
questions:  How  does  the  response  of  a  turbulent  shear  flow  to  compression  differ  from 
that  of  isotropic  turbtJence  ?  How  does  a  shear  flow  respond  to  oblique  compression  as 
compared  to  normal  compression  ? 

One-dimensional  compression  is  used  to  model  the  shock  wave  in  this  chapter.  The 
turbulence  is  constrained  to  be  homogeneous  and  solenoidal.  The  time  scale  of  com¬ 
pression  is  assumed  small  as  compared  to  the  time  scale  of  the  turbulence  (a  reasonable 
assumption  in  shock/turbulence  interaction).  RDT  is  then  used  to  predict  the  response  of 
the  turbulence  to  compression.  As  discussed  in  Chapter  1,  this  approach  is  expected  to 
yield  reasonable  prediction  of  the  interaction  of  solenoidal  turbulence  with  shock  waves  of 
moderate  strength. 

The  roots  of  RDT  applied  to  flows  under  mean  compression  lie  in  the  work  of  Prandtl 
(1933),  Taylor  (1935),  Ribner  and  Tucker  (1953)  and  Batchelor  and  Proudman  (1954). 
Motivated  by  the  passage  of  turbulence  through  a  wind-tvmnel  contraction,  these  work¬ 
ers  developed  RDT  for  turbulent  flows  subjected  to  rapid  irrotational  distortions.  More 
recently,  Lee  (1989)  has  performed  a  detailed  analysis  of  the  response  of  the  turbulence 
to  axisymmetric  strain  and  dilatation.  The  above  studies  assumed  the  turbulence  to  be 
incompressible;  recent  studies  by  Durbin  and  Zeman  (1992)  and  Cambon,  Coleman  and 
Mansour  (1992)  have  applied  RDT  to  the  compression  of  compressible  turbulence.  While 
Durbin  and  Zeman  have  examined  the  pressure  fluctuations  and  the  pressure-dilatation 
correlation  in  the  limit  of  vanishing  turbulent  Mach  Number  (nearly  solenoidal  turbu¬ 
lence),  Cambon  et  al.  have  considered  finite  turbulent  Mach  numbers  and  shown  the 
negligible  effect  of  pressure  fluctuations  at  high  turbulent  Mach  numbers. 

All  the  above  mentioned  studies  considered  the  compression  of  isotropic  turbulence. 
To  the  best  of  our  knowledge,  the  response  of  shear  flows  to  rapid  compression  has  not 
been  studied  before.  Precedence  in  the  use  of  RDT  to  examine  the  rapid  straining  of 
anisotropic  turbulence  may  be  foimd  in  the  work  of  Townsend  (1980),  Sreenivasan  and 
Narasimha  (1978),  Maxey  (1982)  and  Cambon  (Private  communication). 

This  chapter  is  organized  as  follows.  A  brief  outline  of  RDT  in  Section  2.2  is  followed 
in  Section  2.3  by  the  response  of  turbulent  shear  flows  to  normal  compression.  Section 
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2.4  deals  with  the  oblique  compression  of  sheared  turbulence.  Application  of  the  analysis 
to  predict  shock/turbulence  interaction  is  described  in  Section  2.5  and  Reynolds  stress 
modeling  of  the  normal  compression  problem  is  examined  in  Section  2.6.  Finally,  the 
salient  results  are  summarized  in  Section  2.7.  Most  of  the  work  in  this  chapter  is  described 
in  Mahesh,  Lele  and  Moin  (1993)  and  Mahesh,  Lele  and  Moin  (1994a). 

2.2  Theoretical  Procedure 


Rapid  Distortion  Theory  combines  linearization  of  the  governing  equations  with  sta¬ 
tistical  averaging  to  describe  the  statistical  evolution  of  turbulence  imder  rapid  mean 
distortion.  The  formal  development  of  the  RDT  approximation  is  outlined  in  reviews  such 
as  those  by  Savill  (1987)  and  Hunt  and  Carruthers  (1990).  When  the  time-scale  of  the 
mean  distortion  is  much  smaller  than  that  of  the  turbulence,  then  the  turbulence  has  no 
time  to  interact  with  itself.  This  allows  the  neglect  of  all  terms  in  the  governing  equations 
that  involve  viscosity  or  the  product  of  fluctuations  yielding  a  set  of  evolution  equations 
that  are  linear  in  the  fluctuations.  An  alternative  formulation  (Pearson,  1959;  Deissler, 
1961)  corresponding  to  low  Reynolds  number  or  ‘weak’  turbulence  requires  retention  of 
the  viscous  terms.  We  neglect  the  viscous  terms  in  om:  analysis. 


Linearization  of  the  continuity  and  momentum  equations  yields  the  following  set  of 
equations. 


dui 

dxi 


=  0 


dui  dui  dUi  1  dp 


(2.1a) 

(2.16) 


where,  Ui  and  p  are  the  mean  velocity  and  density  respectively,  and  u,-  and  p  are  the 
fluctuating  velocity  fleld  and  pressme  respectively.  For  irrotational  mean  distortions,  it  is 
convenient  to  solve  the  linearized  vorticity  equation  which  reduces  to  Cauchy’s  equation. 


Note  that  the  fluctuations  are  assumed  to  be  solenoidal.  If  the  mean  distortion  is 
solenoidal,  then  the  above  equations  correspond  to  linearization  of  the  incompressible 
Navier  Stokes  equations.  However,  if  the  mean  field  is  dilatational  (as  for  a  one- dimensional 
compression),  then  the  governing  equations  are  the  compressible  Navier  Stokes  equations 
and  the  above  set  of  equations  therefore  describe  evolution  of  the  solenoidal  component 
of  the  compressible  flow  field  under  the  assumption  that  its  evolution  is  independent  of 
the  dilatational  component.  (Alternatively,  equations  (2.1a)  and  (2.16)  describe  the  evolu¬ 
tion  of  a  compressible  flow  field  with  spatially  uniform  but  time-dependent  density.)  The 
solenoidal  and  dilatational  component  are  coupled  even  in  the  linear  limit  for  rotational 
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mean  flows  (Blaisdell,  Mansoiir  and  Reynolds,  1991).  However  for  irrotational  mean  dis¬ 
tortions,  the  solenoidal  component  may  be  assumed  to  evolve  independently  in  the  limit 
Am  =  ToL/a  <C  1  where,  To  is  the  mean  strain  rate,  X  is  a  turbulence  lengthscale  and  a  is 
the  mean  speed  of  sound  ;  i.e.,  the  turbulence  is  nearly  incompressible  (Cambon,  Coleman 
and  Mansour,  1992). 

The  assumption  of  homogeneity  constrains  the  mean  velocity  gradient  to  be  uniform; 
i.e.,  the  mean  velocity  is  of  the  form  Ui  =  Aik{t)xk.  For  incompressible  turbulence,  this 
is  the  only  reqmrement  for  homogeneity;  since  as  discussed  above,  we  consider  evolution 
of  the  solenoidal  component  of  of  a  compressible  flow  field,  we  constrain  the  mean  field 
to  satisfy  homogeneity  for  compressible  tmrbulent  fluctuations.  As  a  result,  in  addition  to 
a  imiform  velocity  gradient,  the  mean  field  satisfies  the  compressible  Euler  equations  and 
has  uniform  pressure  and  density  (Blaisdell,  Mansour  and  Reynolds,  1991). 

The  procedure  for  solution  of  equations  (2.1a)  and  (2.1b)  is  fairly  well  established.  One 
method  of  solution  (Batchelor  and  Proudman,  1954)  involves  using  a  Fourier  representation 
where  the  wavenumber  changes  with  time  as 


k 

p(x,t)  = 

k 


where. 


dka 

dt 


kjAjct  —  0 


(2.2) 


An  alternative  equivalent  method  of  solution  (Rogallo,  1981)  is  to  transform  coordinates 
to  a  system  that  deforms  with  the  mean  field;  i.e.. 


(i  =  Bik{t)xk  T  =  t 


where, 

+  AjkBnj  —  0  (2-3) 

The  transformed  equations  are  then  solved  using  conventional  Fourier  representation. 
Knowledge  of  the  Fourier  coefiicients  enables  computation  of  the  energy  spectrum  ten¬ 
sor  which  is  then  integrated  over  all  wavenumbers  to  determine  the  Reynolds  stresses. 

Homogeneous  RDT  being  an  initial-value  problem,  requires  specification  of  the  initial 
energy  spectrum.  As  mentioned  previously,  earlier  studies  assume  isotropic  initial  spec¬ 
trum.  Our  analysis  differs  in  that  the  compression  of  anisotropic  turbulence  is  considered. 
Details  of  the  analysis  are  presented  in  the  following  sections. 
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2.3.  Normal  Compression  of  Sheared  Turbulence 


The  response  of  a  turbulent  shear  flow  to  normal  compression  is  discussed  in  this 
section.  Description  of  the  problem  and  the  solution  procedure  is  outlined  in  Section  2.3.1. 
Results  of  the  analysis  are  presented  and  discussed  in  Sections  2.3.2  and  2.3.3  respectively. 

2.3.1  Problem  Formulation 


Figure  2.1  shows  a  schematic  of  the  normal  compression  of  a  turbulent  shear  flow. 
One- dimensional  compression  is  characterized  by  the  following  mean  fleld. 


Ui 


To 

l  +  Fot"^^’ 


172  =  173  =  0 


(2.4a) 


P  = 


Po 

1  -h  Tot 


P  = 


Po 
[1  -1- 


(2.46) 


(2.4c) 


where  Ui,U2,  Uz  are  the  mean  velocity  components  in  the  xi ,X2,Xz  directions  respectively; 
p  and  P  are  the  mean  density  and  pressure  respectively,  and  are  uniform  in  space.  To  is 
negative  for  compression  and  positive  for  expansion. 


Consider  isotropic  turbtdence  that  is  subjected  to  rapid  homogeneous  shear  z.e., 


17i  =  5oX2,  P  =  Po,  P  =  Po 

(2.6) 

At  a  non-dimensional  time  0q  =  Soto  during  application  of  shear, 
dimensional  compression  i.e.,  for  t  >  to, 

we  introduce  the  one- 

To  So 

l  +  To{t  —  to)  ^  '  1 -I- ro(t  —  to)^^ 

(2.6a) 

_  Po 

(2.66) 

1  +  ro(<  —  to) 

p  - 

[1  +  ro(t  —  to)]"^ 

(2.6c) 

The  shear  rate  changes  with  time  during  the  application  of  compression  to  satisfy  the 
compressible  Euler  equations.  Note  that  the  shear  is  along  X2  while  the  compression  is 
along  xi ;  hence  the  terminology  ‘normal  compression’. 


Under  RDT,  the  Reynolds  stresses  before  compression  depend  only  upon  the  total 
shear  /3o;  their  subsequent  evolution  depends  upon  ^o,  the  ratio  of  shear  rate  to  the  rate 
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of  compression  (Sq/Tq)  and  the  total  volumetric  strain  (p/po).  Since  our  interest  is  in 
the  compression  of  shear  flows,  we  consider  the  regime  where  the  shear  sets  up  the  initial 
anisotropic  field  and  is  negligible  as  compared  to  the  subsequently  applied  compression; 
*.e.,  5o/ro  ^  1.  In  this  report  Sq/Tq  is  0.1  for  all  cases  presented.  Lower  values  of  S’d/Fq 
(for  eg,  Sq/Tq  =  0.01)  yielded  results  identical  to  those  shown  here.  Thus,  we  effectively 
consider  the  rapid  one-dimensional  compression  of  sheared  turbulence. 

The  RDT  equations  are  solved  by  transforming  coordinates  to  a  system  that  deformed 
with  the  mean  field.  Fourier  representation  in  the  transformed  coordinates  is  used  to  derive 
ordinary  differential  equations  that  describe  the  evolution  of  the  energy  spectrum  tensor. 
The  equations  are  numerically  advanced  in  time  and  at  each  time  step,  the  energy  spectrum 
tensor  is  numerically  integrated  over  all  wavenumbers  to  compute  the  Reynolds  stresses. 

As  indicated  above,  the  mean  field  is  pure  shear  for  t  <  to.  The  RDT  equations  for 
pure  shear  have  been  solved  analytically  (Moffatt,  1965;  Townsend,  1970)  and  hence  we 
do  not  give  the  evolution  equations  for  the  energy  spectrum  tensor.  We  only  note  that  the 
coordinate  transform  used  is  as  follows, 

(i=xi-SoiX2,  ^2  =  ^2,  6=2:3,  T=t  (2.7) 

At  time  to  =  /^q/Sq,  the  coordinate  transformation  changes  to  accommodate  the  compres¬ 
sion.  The  new  transformation  is  given  by 

6  =  ^  ~  ~  6=2:2,  6=2:3,  r  =  t  (2.8) 


For  t  >  to,  the  RDT  equations  are  transformed  to  the  above  coordinate  system.  Using 
the  Fourier  representation. 


K 

(2.9a) 

K6'r)  =  Y^p{K,T)e''^i^> 

(2.96) 

the  RDT  equations  tahe  the  following  form  for  t  >  to: 

dui  1  r  ip 

^  =  T— pr-^  -ro«i-5oU2  +  «i- 
dr  1-1-  For'  p 

(2.10a) 

du2  5or'  1  ip 

dr  1  -f  For'  p 

(2.106) 

II 

(2.10c) 
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(2.10J) 


^  _  /^iro  u\  +  [/ci5o  +  («i5o  —  K2ro)  (1  +  Fot')  ]  U2  —  /csFo  (1  +  For')  uz 
p  «1  [1  +  -  2  kiK2  Sqt’  (1  +  For')  +  (nj  +  k^)  (1  +  For')^ 

where  the  variable  t'  =  t  —  tq.  Note  that  the  expression  for  pressure  is  obtained  through 
elimination  of  the  continuity  equation.  Also,  Fo  =  0  in  the  above  equations  would  describe 
the  evolution  for  t  <to  (Pure  shear). 

The  RDT  equations  for  a  pmre  compression  allow  analytical  solution  as  do  the  equa¬ 
tions  for  pure  shear.  However,  as  seen  above,  the  combination  of  shear  and  compression 
yields  a  formidable  set  of  equations  which  we  have  not  been  able  to  solve  analytically. 
Instead  equations  (2.10a)  through  (2.10d)  were  used  to  derive  evolution  equations  for  the 
energy  spectrum  tensor  Eij{K,T),  defined  as  Eij  =  where  the  superscript  refers 
to  the  complex  conjugate.  For  reasons  of  brevity,  the  equation  for  Eij  is  not  reproduced 
here;  it  has  the  following  form: 

dEii 

+  CjiEii  (2.11) 

The  above  system  of  equations  is  integrated  numerically  to  compute  Eij{K,  r).  Eij  is  then 
integrated  over  all  wavenumbers  to  compute  the  Reynolds  stress  tensor  J2,y(r)  defined  as 

^iij  —  XXiXJLj^  1.6., 

=  jEij{K,T)d^K  (2.12) 

The  integration  is  carried  out  in  polar  coordinates 


Ki  =  Kcos^,  K2  =  Ksm(f>cos6,  Kz  =  K  sin  k  =  sin <f>  d<f> d9 dn 

where,  k  varies  from  0  to  oo,  (f),  from  0  to  tt  and  6,  from  0  to  27r.  Note  that  since  at  r  =  0, 
the  energy  spectrum  tensor  is  assumed  to  be  isotropic;  i.e., 

it  can  be  shown  that  Rijir )  is  independent  of  the  initial  three-dimensional  energy  spectrum 
tensor,  E{k)  and  the  magnitude  of  the  wavenumber  vector. 

2.3.2  Results 

The  analysis  described  above  essentially  applies  RDT  to  the  compression  of  a  shear 
flow  where  the  energy  spectrum  tensor  of  the  shear  flow  is  itself  the  RDT  solution  to  initially 
isotropic  turbulence  subjected  to  homogeneous  shear.  The  use  of  RDT  to  obtain  the  energy 
spectrum  tensor  that  characterizes  a  shear  flow  is  motivated  by  Townsend  (1970)  where  the 
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RDT  predictions  of  Reynolds  stress  correlations,  spectra  and  non-dimensional  Reynolds 
stresses  in  turbulent  shear  flows  were  compared  to  experiment  and  good  agreement  was 
observed. 

This  is  illustrated  in  Figure  2.2  where  Rij/q^  as  predicted  by  RDT  applied  to  the 
rapid  shear  of  initially  isotropic  turbulence  is  shown.  The  quantity  plotted  on  the  abscissa 
is  the  total  shear  ^  defined  as  /3  =  Sot.  Note  that  after  a  moderate  amount  of  total  shear  (/? 
between  2  and  3),  the  nondimensional  Reynolds  stresses  are  quite  close  to  values  obtained 
in  shear  flows  such  as  homogeneous  shear  flow,  turbulent  channel  flow  (Moin,  1988)  and 
boundary  layers  (Townsend,  1970).  For  example,  in  homogeneous  shear  flow  (Moin,  1988), 
7211/5^,^22/9^5-^33/$^  and  i2i2/$^  are  0.54,  0.16,  0.31  and  0.14  respectively.  For  a  total 
shear  of  2.5,  RDT  yields  values  of  0.51,  0.12,  0.36  and  0.16  respectively.  Townsend  (1970) 
notes,  ‘The  good  agreement  of  the  non-dimensional  Reynolds  stresses  suggests  that  the 
large  eddies  in  turbulent  shear  flows  have  very  nearly  the  shape  suggested  by  RDT’.  This 
conclusion  is  further  supported  by  direct  numerical  simrdation  (Lee,  Kim  and  Moin,  1990) 
of  tmbulence  subjected  to  very  high  shear  rates,  where  structures  similar  to  near  wall 
‘streaks’  were  observed. 

The  evolution  of  sheared  turbulence  upon  compression  is  next  examined.  The  abscissa 
in  all  plots  shown  is  the  total  volumetric  strain,  p/po.  normalized  with  its  value  at 
the  start  of  compression  is  plotted  against  the  total  volumetric  strain  in  Figure  2.3.  The 
different  curves  correspond  to  different  values  of  total  initial  shear  {^0).  /^o  =  0  corresponds 
to  the  compression  of  initially  isotropic  turbulence,  q^  is  seen  to  amplify  upon  compression 
with  the  ampHfication  ratio  increasing  as  the  initial  total  shear  increases.  The  amplification 
of  turbulent  kinetic  energy  upon  normal  compression  is  thus  higher  for  a  shear  flow  than 
it  is  for  isotropic  turbulence. 

The  evolution  of  the  components  of  turbulent  kinetic  energy  is  plotted  in  Figure  2.4.  /?o 
was  set  to  3  for  this  case.  All  three  components  are  seen  to  amplify  with  the  xi  (direction 
of  compression)  component  being  amplified  the  most.  The  reason  for  the  preferential 
amplification  of  the  xi  component  is  due  to  the  fact  that  it  is  directly  ‘produced’  by  the 
compression  while  the  other  components  ampHfy  through  the  redistributive  nature  of  the 
pressure-strain  correlation  in  the  Reynolds  stress  equations.  The  importance  of  the  initial 
anisotropy  on  the  amplification  of  uf  is  gauged  by  comparing  the  amplification  of  Uj  of 
initially  sheared  turbulence  to  that  of  initially  isotropic  turbulence  (Figure  2.5).  Note  that 
the  amplification  of  u\  in  the  shear  flow  is  much  higher  than  that  of  isotropic  turbulence. 
For  example,  at  a  density  ratio  of  3,  the  amplification  ratio  is  3.4  for  initially  isotropic 
turbulence  and  4.85  for  initially  sheared  turbulence  -  an  increase  of  about  43  %. 
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As  mentioned  earlier,  shear  flows  in  equilibrium  are  characterized  by  typical  values  of 
Rij/(^.  It  is  of  interest  to  see  how  these  values  change  upon  normal  compression.  Figure 
2.6  shows  the  the  diagonal  terms  of  the  tensor.  Compression  increases  the  contribution 
of  to  the  turbulent  kinetic  energy  while  decreasing  that  of  U2  and  uz.  Also  note  (over 
the  range  of  total  volumetric  strain  shown)  that  the  ordering  of  kinetic  energy  components 
(ui  >  Uz  >  U2)  is  retained  upon  normal  compression. 

The  evolution  of  Ru/g.^  upon  compression  is  plotted  in  Figme  2.7.  The  three  curves 
correspond  to  different  values  of  total  initial  shear.  Upon  normal  compression,  Ri2/q'^ 
decreases  in  magnitude  and  for  sufficiently  large  total  volumetric  strains,  it  changes  sign. 
This  trend  was  flrst  observed  by  Gambon  (Private  communication),  and  is  hastened  upon 
increasing  the  initial  total  shear.  The  terms  in  the  Reynolds  shear  stress  evolution  equation 
provide  an  explanation  for  this  behavior.  The  evolution  equation  for  R12  is. 


To 

1  +  Fot 


Ri2  — 


So 

1  +  Tot 


R22  +  7!’12 


(2.13) 


where  -Kij  is  the  pressure-strain  correlation  defined  as  7r,y  =  p{uij  +  Uj^i)/p.  Note  that 
both  the  strain  and  shear  production  terms  tend  to  increase  the  magnitude  of  R12  (make 
it  more  negative).  The  tendency  of  |i2i2|  to  decrease  must  therefore  be  due  to  the  pressure 
strain  correlation.  Figure  2.8  illustrates  the  evolution  of  terms  in  the  budget  of  R12  for  the 
case  with  /Sq  =  3.  We  see  that  the  tendency  of  R12  to  decrease  upon  normal  compression  is 
due  to  amplification  of  the  pressme-strain  correlation  and  the  consequent  upsetting  of  the 
initial  balance  between  ‘production’  and  the  pressure-strain  correlation  in  the  shear  flow. 

2.3.3  Interpretation  of  the  results 


We  have  seen  that  the  response  of  a  shear  flow  to  normal  compression  is  quite  different 
from  that  of  isotropic  turbulence.  A  noticeable  feature  is  the  higher  amplification  of  and 
in  the  shear  flow.  In  this  section,  we  attempt  to  explain  this  observation  by  posing  the 
following  question,  ‘What  aspect  of  the  initial  anisotropic  field  is  important  in  determining 
the  evolution  of  kinetic  energy  during  a  one- dimensional  compression  ?’ 

The  evolution  equations  for  uf  and  during  rapid  compression  are  given  by 


duj 

dr 


-2 


To 

1  -f  For 


uj  +  TTn 


(2.14a) 


dq^ 

dr 


-2 


To 

1  -f  For 


a; 


(2.146) 
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The  pressure-strain  correlation  ‘taJces’  energy  from  and  redistributes  it  among  the  other 
components.  The  amplification  of  u\  and  hence,  ^  through  equation  (2.146),  would  in¬ 
crease  if  the  pressure-strain  correlation  were  to  decrease.  Also,  since  in  the  RDT  limit,  the 
pressure-strain  correlation  is  the  only  term  in  the  budget  apart  from  the  production  term, 
a  higher  amplification  rate  automatically  implies  a  lower  pressure-strain  correlation. 

The  relative  magnitude  of  the  pressure-strain  term  in  the  energy  budget  may  be 
gauged  from  its  short  time  behavior.  It  is  easily  shown  that  in  turbulence  subjected  to 
rapid  one-dimensional  compression  tth  evolves  as, 

x„  =  -2r.c  j  (2.18) 

where,  c  =  p/po-  ■'I’li  thus  strongly  depends  upon  the  spectral  distribution  of  the  ui 
component  of  velocity.  Note  that  as  the  energy  distribution  in  uj  moves  to  smaller  «i ,  the 
magnitude  of  ttu  decreases.  The  small  time  evolution  of  tth  may  be  estimated  by  replacing 
Ell  in  equation  (2.15)  with  its  initial  value,  E^^.  This  small  time  evolution  shows  how 
the  spectral  distribution  of  in  the  initial  field  affects  the  rate  of  change  of  uf  when 
anisotropic  turbulence  is  compressed. 

Since  mean  shear  tends  to  stretch  the  turbtdence  in  the  streamwise  direction,  ui  in 
a  shear  flow  has  energy  at  smaller  ki  as  compared  to  isotropic  turbulence  and  hence 
according  to  equation  (2.15),  a  lower  initial  value  of  ttu  .  This  is  shown  in  Figure  2.9  where 
the  ratio  of  tth  to  the  production  term  in  the  Rn  equation  is  plotted.  The  pressiire-strain 
correlation  is  indeed  seen  to  be  smaller  when  shear  flow  is  compressed. 

As  one  continues  the  compression,  the  energy  in  ui  is  moved  to  larger  ki  and  hence,  the 
pressure-strain  correlation  would  progressively  become  more  important.  Analysis  shows 
that  in  the  limit  of  infinite  p/po,  the  ratio  of  tth  to  the  production  term  is  independent 
of  the  initial  energy  spectrum;  it  is  equal  to  1.  This  is  shown  as  follows.  For  turbulence 
subjected  to  one-dimensional  compression,  under  RDT, 


Bii(k)  = 


C^(«i/k)^  +  («2/«)^  +  («3/«)^ 


(2.16) 


where,  c  =  p/po- 

Denoting  the  denominator  in  the  above  expression  by  D,  the  rate  of  change  of  En  is 
given  by, 

dEii  E?,r^dc^  .dD\  ^217) 


11 


dr 


D2 


^dc^  ndD' 
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Noting  that  dc/di  =  -Foc^  and  considering  very  large  compressions,  it  can  be  shown  that, 


lim 

C— ►OO 


dr 


=  0 


(2.18) 


Since  En  when  integrated  over  all  wavenumbers  yields  iln,  the  above  equation  implies 
that  in  the  limit  of  very  large  c,  duf/dr  =  0.  Substitution  into  equation  (2.14a)  shows 
that  asymptotically,  tth  =  Production.  Note  that  this  relation  is  independent  of  the 
initial  energy  spectrum.  The  asymptotic  behavior  shows  the  diminishing  effect  of  initial 
anisotropy  upon  the  pressure-strain  correlation  as  time  progresses.  This  is  supported  by 
Figure  2.9  where  the  two  curves  tend  towards  each  other  as  c  increases.  Integration  to 
larger  time  shows  that  the  two  curves  do  indeed  asymptote  to  1. 

We  see  that  when  anisotropic  turbulence  is  compressed  (in  the  xi  direction),  the 
feature  of  the  initial  field  that  determines  kinetic  energy  evolution  is  E^^(ki).  As  the 
energy  in  u®  moves  to  smaller  ki,  the  relative  importance  of  ttu  decreases  resulting  in 
larger  amplification  of  Uj.  The  amplification  of  in  addition  to  being  influenced  by  F/jj, 
is  also  influenced  by  the  initial  value  of  Uj/ .  This  is  shown  as  follows.  The  kinetic  energy 
equation  for  a  one- dimensional  compression  (Equation  2.146)  may  be  rewritten  as 


1  __  To  ^  ^ 

29o'  dr  l  +  Fori??, 


(2.19) 


Note  that  as  Rii/q^  increases,  the  amplification  of  q^  increases.  This  explains  the  trend 
seen  in  Figure  2.3.  In  addition  to  having  energy  at  lower  ki,  a  shear  flow  that  is  sub¬ 
jected  to  normal  compression  also  has  a  larger  fraction  of  its  energy  along  the  direction  of 
compression,  resulting  in  a  higher  amplification  of  5^  as  compared  to  isotropic  turbulence. 
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2.4.  Oblique  Compression  of  Sheared  Turbulence 


Thus  far,  we  have  examined  the  effect  of  compression  on  sheared  turbulence  normal 
to  the  direction  of  shear.  Since  shear  flows  are  anisotropic,  one  would  expect  the  direc¬ 
tion  of  compression  relative  to  the  shear  to  be  an  important  parameter;  z.e.,  the  oblique 
compression  of  a  shear  flow  would  yield  results  different  from  normal  compression.  In  the 
section  that  follows,  we  discuss  the  procedmre  (Section  2.4.1)  and  results  (Section  2.4.2)  of 
RDT  applied  to  the  problem  of  oblique  compression. 

2.4.1.  Problem  Formulation 

Under  RDT,  the  obliquity  of  compression  may  be  characterized  by  the  angle  6  between 
the  direction  of  compression  and  the  direction  of  the  upstream  shear  flow.  Note  that  ^  =  0 
corresponds  to  the  normal  compression  discussed  in  the  previous  section.  Recall  that  the 
normal  compression  problem  was  formulated  by  considering  sheared  turbulence  subjected 
to  simultaneous  shear  and  compression  in  the  regime  where  the  compression  rate  was  much 
higher  than  the  rate  of  shear.  The  sole  reason  for  retention  of  the  shear  during  compression 
was  to  avoid  the  destruction  of  vorticity.  However,  one  might  adopt  the  point  of  view  that  if 
the  mean  shear  is  not  important  during  the  process  of  compression,  then  it  may  be  ignored 
except  for  its  effect  on  the  initial  spectrum.  This  is  the  approach  adopted  in  this  section 
in  applying  RDT  to  the  oblique  compression  of  sheared  turbulence.  We  essentially  repeat 
RDT  for  turbulence  subjected  to  a  one-dimensional  compression.  However  the  initial 
energy  spectrum  instead  of  being  isotropic  is  the  RDT  solution  to  isotropic  turbulence 
subjected  to  homogeneous  shear.  Also  as  mentioned,  the  initial  mean  shear  is  assumed  to 
be  at  angle  6  to  the  subsequently  applied  compression.  The  validity  of  this  approach  was 
verified  by  comparison  of  the  results  for  ^  =  0  to  the  more  completely  formulated  normal 
compression  problem. 

Figure  2.10  shows  a  schematic  of  the  the  oblique  compression  problem.  Note  that 
the  shear  is  assumed  to  be  in  the  X2  direction.  For  reasons  of  convenience,  the  time- 
dependent  wavenumber  approach  (Batchelor  and  Proudman,  1954)  was  used  to  solve  the 
RDT  equations.  The  RDT  solution  to  the  oblique  compression  problem  requires  knowledge 
of  the  evolution  of  the  energy  spectrum  tensor  when  isotropic  turbulence  is  subjected 
to  homogeneous  shear  and  the  transfer  function  of  the  vorticity  spectrum  tensor  when 
turbulence  (not  necessarily  isotropic)  is  subjected  to  a  one-dimensional  compression.  As 
indicated  earlier,  these  problems  have  been  solved  analytically.  Assuming  knowledge  of 
these  transfer  functions,  the  procedure  for  RDT  applied  to  the  oblique  compression  problem 
is  as  follows: 
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1.  Consider  isotropic  turbulence  that  is  subjected  to  homogeneous  shear  in  the  direc¬ 
tion.  Denote  the  wavenumber  vector  at  <  =  0  by  ko-  For  a  given  total  shear  (defined 
in  the  previous  section)  compute  the  energy  spectrum  tensor  at  the  end  of  shear. 

2.  With  reference  to  Figure  2.10,  rotate  the  energy  spectrum  tensor  and  the  wavenumber 
vector  by  the  angle  9.  This  aligns  the  field  with  the  compression. 

3.  From  the  rotated  energy  spectrum  tensor  and  the  wavenumber  vector,  obtain  the 
vorticity  spectrum  tensor  in  the  rotated  coordinates. 

4.  Using  the  transfer  functions  for  the  compression  problem,  obtain  the  vorticity  spec¬ 
trum  tensor  and  wavenumber  vector  after  compression. 

5.  Using  the  inverse  of  the  relation  in  3,  obtain  the  energy  spectrum  tensor  after  com¬ 
pression  in  the  rotated  coordinates.  Integrate  the  energy  spectrum  tensor  over  ko  to 
obtain  the  Reynolds  stress  tensor  in  rotated  coordinates. 

6.  Rotate  the  Reynolds  stress  tensor  back  by  the  angle  9.  This  yields  the  Reynolds  stress 
tensor  after  compression  in  the  original  coordinate  system. 

2.4.2.  Results 

We  describe  in  this  section  some  of  the  results  of  RDT  applied  to  the  oblique  com¬ 
pression  of  sheared  turbulence.  The  effect  of  the  oblique  angle  9^  is  gauged  by  comparison 
to  the  previously  discussed  normal  compression  problem.  The  amplification  of  the  stream- 
wise  component  of  turbulent  kinetic  energy  upon  compression  is  shown  in  Figure  2.11  for 
different  angles  of  obliquity. 

The  initial  condition  corresponds  to  sheared  turbulence  with  /?o  =  3.  Recall  that  ^  =  0 
corresponds  to  normal  compression.  Note  the  decrease  in  amplification  ratio  with  increas¬ 
ing  magnitude  of  oblique  angle.  The  effect  of  oblique  compression  on  is  shown  in  Figure 
2.12.  Once  again,  the  amplification  of  decreases  as  the  oblique  angle  increases.  In  the 
previous  section,  we  emphasized  the  importance  of  the  £’ii(ki)  and  u\/q^  on  the  amplifi¬ 
cation  of  kinetic  energy.  If  oblique  compression  is  viewed  in  a  coordinate  system  aligned 
with  the  compression,  we  see  that  the  effect  of  oblique  compression  may  be  explained  by 
the  initial  value  of  Egging)  and  u^/q^  ,  where  9  refers  to  the  direction  of  compression. 

In  the  previous  section,  we  noted  how  normal  compression  coiild  significantly  change 
the  nondimensional  Reynolds  stresses  of  a  shecir  flow.  The  effect  of  oblique  compression 
on  the  diagonal  elements  of  Rij/q^  is  shown  in  Figures  2.13,2.14  and  2.15.  It  is  clear  that 
obliquity  of  compression  has  a  significant  effect  on  the  evolution  of  the  Reynolds  stresses. 
The  qualitative  difference  in  the  curves  for  different  oblique  angles  is  striking  to  note. 
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Figure  2.16  shows  the  evolution  of  the  nondimensional  Reynolds  shear  stress  upon 
oblique  compression.  Recall  the  effect  of  normal  compression  to  decrease  and  even 

change  its  sign.  In  comparison,  we  see  that  this  trend  decreases  as  the  compression  becomes 
more  oblique  and  for  large  oblique  angles,  slightly  amplifies. 

2.5  Application  to  shock/turbulence  interaction 

We  illustrate  apphcation  of  the  analysis  to  shock/turbulence  interaction. 

Recall  that  the  Reynolds  stress  tensor  during  compression  of  a  shear  flow  is  of  the 
form: 

Rij  =  Rij  {pi po;  ,^0 , 0)  (2.20) 

Note  that  ^o/Fo  is  assumed  neghgible  in  writing  the  above  expression.  The  parameter  6 
m  the  above  expression  is  determined  by  the  obhque  angle  of  the  shock  wave.  Figure  2.17 
illustrates  a  couple  of  examples  that  correspond  to  positive  and  negative  9.  The  choice 
of  /3o  depends  upon  the  extent  of  shear  (quantified  by  Sq^je)  in  the  flow  upstream  of  the 
shock  wave.  As  discussed  before  =  2.5  seems  appropriate  for  shear  flows  in  equilibrium. 
Higher  values  of  may  be  chosen  for  strongly  sheared  flows.  The  variable  p/p^  in  the 
above  expression  is  the  total  amount  of  compression  the  shear  flow  experiences.  It  is 
equated  to  the  density  ratio  across  the  shock  wave  which  the  Rankine-Hugoniot  equations 
relate  to  the  shock-normal  Mach  number,  Mn.  i.e. 

P  _  {■y  +  l)Ml 

Po  2-t-(7-l)M2  (2-21) 

Thus,  once  the  strength  and  inclination  angle  of  the  shock  wave  is  known,  its  effect  on  the 
shear  flow  can  be  predicted. 

Consider  the  foUowing  two  examples:  the  interaction  of  a  shear  flow  with  a  Mach 
1.5  normal  shock  and  shock/turbulence  interaction  in  a  compression  corner  with  free- 

stream  Mach  number  2.9  and  wedge  angle  8°.  The  compression  corner  corresponds  to  an 
experiment  by  Smits  and  Muck  (1989). 

The  evolution  of  Rij  and  Rijjq^  is  predicted  for  both  cases.  The  incident  turbulence 
IS  assumed  to  be  in  equilibrium  and  hence  /3o  was  taken  as  2.5.  For  the  normal  shock 
problem,  it  is  clear  that  and  ^  are  1.5  and  0°  respectively.  In  the  case  of  the  compression 
corner,  these  parameters  are  obtained  by  assuming  that  the  flow  is  turned  across  a  single 
shock  wave.  Solution  of  the  corresponding  oblique  shock  equations  yields  =  1.29  and 
e  =  63.65^  Another  point  to  be  noted  is  that  the  coordinate  system  in  our  analysis  was 
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aligned  with  the  direction  of  mean  shear.  However  experimental  coordinate  systems  are 
aligned  with  the  local  wall  direction.  As  a  result,  the  compression  comer  problem  requires 
that  the  Reynolds  stresses  predicted  by  our  analysis  be  rotated  by  the  wedge  angle  to  be 
consistent  with  experiments.  The  results  are  tabulated  below. 


Rii 

R22 

R33 

RuIq^ 

Ri2fq^ 

Normal 

2.60 

1.30 

1.60 

0.64 

-0.07 

Oblique 

1.05 

2.22 

1.40 

0.40 

-0.21 

2.6  Application  to  Reynolds  Stress  modeling 


The  analysis  described  in  the  preceding  sections  reveals  significant  effects  of  anisotropy 
on  the  response  of  shear  flows  to  compression.  The  spectral  distribution  of  energy  prior 
to  compression  is  seen  to  strongly  influence  the  Reynolds  stresses  through  the  pressure- 
strain  correlation.  In  this  section,  we  consider  Reynolds  stress  modeling  of  the  normal 
compression  problem  to  see  if  the  model  equations  reproduce  the  evolution  as  predicted  by 
RDT.  Since  the  pressure- strain  correlation  is  the  only  term  that  needs  to  be  modelled,  the 
RDT  limit  is  a  natural  choice  to  evaluate  models  of  the  rapid  pressure-strain  correlation. 


We  evaluate  models  of  the  form:  7r,j  =  T^ij{hij,  S*j,Q,ij).  Linear,  quadratic  and  cubic 
forms  of  these  models  exist.  We  only  consider  the  linear  version  of  these  models.  The 
objective  here  is  to  see  how  these  models  perform  inspite  of  the  absence  of  explicit  spectral 

information.  Consider 


TTtj 


+  C2  ( Stibij  +  S*ibii  - -St^bmiSij  J  +  C3 1  +  ^jibii 


(2.22) 


The  models  of  several  workers;  eg  Launder,  Reece  and  Rodi  (),  Naot,  Shavit  and  Wolf- 
shteinO  and  Zeman()  are  of  the  above  form,  differing  only  in  the  constants  C2  and  C3. 


During  simultaneous  shear  and  compression  (Section  2.3.1)  S*j  and  Q,ij  are  given  by, 

0 

r.  I  .  t  \  ^  I'o 

5*  = 


*  _  ^  0  I  1  n  I  Q  -  = - — 

l  +  rot'l^ro  ^3  ‘J  i  +  Tot' 


O' 

2  To  ^ 
0 
0 


where  t'  =t -  to.  In  the  limit  5o/ro  <  1, 

To 


fiij  =  0 


(2.23) 


(2.24) 
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We  consider  the  evolution  of  /In  and  Ri^-  Substitution  of  S*j  and  Qij  into  equation  (2.22) 
yields 


7!-ii  = 


3^  l  +  Fof 


+  ^^2611 


(2.25) 


Note  that  the  trace-free  property  of  bij  is  used  to  obtain  the  above  expression.  Also 

To 


1  2 

7ri2  = 


3^  14- Tot' 


6*2612 


(2.26) 


The  modeled  equations  may  be  solved  analytically  in  the  RDT  limit  to  obtain  the  Reynolds 
stresses.  Consider 

(2.27a) 


dt 


1  -|-  Tot' 

^  -  -9  To 
dt  ^  ^ 


T,  .  -1.  (2-276) 

l+Tot'  ^  ■ 

Defining  a  variable  /  =  Ru/q^  =  bn  -f  1/3,  the  above  equations  may  be  manipulated  to 
obtain 


#  _  o  To 

dt'  1  Tot' 


(2.28) 


Substituting  for  TTn/q'^, 


To 


dt'  I  +  Tot' 


1,2  6*2- 


(2.29) 


Denoting  /  at  t'  =  0  by  /o ,  the  above  equation  may  be  integrated  to  yield 

T\  —  Kr2 


/  = 


1  —  « 


where,  k  is  defined  as 

JO  -  r2 

ri  and  r2  are  the  roots  of  p  -\-  /(^  ~  l)  +  3  (I  ~  =  0  and  are  given  by 


(2.30) 


(2.31) 


^1,2 


(2.32) 


Once  /  is  known,  may  be  calculated  as  follows.  We  rewrite  equation  (2.276)  as, 

(2.33) 


dt' 


1  Tot' 
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Substituting  for  /  and  integrating,  we  get, 


9 


2 


?0 


(2.34) 


where,  the  subscript  ‘O'  corresponds  to  time  t'  =  0.  Note  that  k  may  be  expressed  in  terms 
of  the  total  volumetric  strain  ;  i.e.. 


/o  -r2  \pq) 


(2.35) 


We  next  compare  the  predicted  values  of  Rn  and  with  the  RDT  evolution.  Our  choice 
of  C2  in  the  above  equations  is  dictated  by  the  asymptotic  behavior  of  the  model  solution. 
Equations  (2.32)  and  (2.35)  show  that  in  the  lirmt  of  infinite  volumetric  strain,  /c  — >  00  and 
hence  f  r^.  Siince  r2  depends  only  upon  C2,  this  suggests  that  the  asymptotic  value  of 
Rii/q^  as  predicted  by  the  model  is  independent  of  the  initial  anisotropy.  This  behavior 
is  in  agreement  with  RDT.  This  independence  of  Rn/q^  of  the  initial  condition  allows  the 
asymptotic  values  to  be  determined  from  the  compression  of  isotropic  turbulence.  It  can 
be  shown  (Lee,  1989)  that  Ru/q^  asymptotes  to  1/2.  Equating  r2  to  1/2  and  solving  the 
resulting  equation  for  C2,  we  get 


(2.36) 


Figure  2.18  compares  the  predicted  evolution  of  R\\  to  the  RDT  solution.  The  com¬ 
pression  of  isotropic  turbulence  (/?o  =  0)  eind  sheared  turbulence  (/?o  =  2.5)  are  both 
plotted.  Note  that  the  model  predicts  Rn  quite  well  when  the  initial  state  is  isotropic. 
However,  the  amplification  of  Rn  in  the  shear  flow  is  considerably  underpredicted.  Similar 
behavior  is  observed  in  Figure  2.19  where  the  evolution  of  q^  is  plotted.  The  evolution 
of  isotropic  turbulence  is  predicted  well  while  the  amplification  of  kinetic  energy  in  the 
shear  flow  is  underpredicted.  This  behavior  of  the  model  solution  is  attributed  to  the  lack 
of  spectral  information  which  prevents  it  from  reproducing  the  suppression  of  ttu  in  the 
shear  flow. 


Similarly  the  model  is  unable  to  reproduce  the  drop  and  change  in  sign  of  Reynolds 
shear  stress.  This  is  shown  as  follows.  The  evolution  of  R12  is  given  by  equation  (2.13)  as 


To 

1  -f-  Fot' 


+  '^12 


(2.37) 
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Using  the  model  expression  for  tti2  els  given  in  equation  (2.26)  we  get, 

In  the  limit  SqITq  <C  1,  this  becomes 

Integration  of  the  above  equation  yields 

i2i2  =  i2;2(l  +  roO"''^'’'  (2.40) 

Since  Fq  is  negative,  the  above  equation  predicts  an  increase  in  the  Reynolds  shear  stress 
upon  normal  compression.  Thus,  presence  of  spectral  information  seems  to  be  particularly 
important  in  predicting  the  evolution  of  i?i2. 

2.7.  Summary 

Homogeneous  rapid  distortion  theory  was  used  to  examine  the  response  of  incom¬ 
pressible,  turbrJent  shear  flows  to  rapid  one-dimensional  compression.  Both  normal  and 
oblique  compressions  (with  respect  to  the  shear)  were  considered.  The  response  of  shear 
flows  to  compression  is  found  to  be  quite  different  from  the  previously  studied,  compression 
of  isotropic  turbulence.  The  differences  are  interpreted  in  a  more  general  framework  and 
the  relevant  parameters  influencing  kinetic  energy  amplification  are  identified. 

£^ii(/ci)  of  the  initial  fleld  is  found  to  determine  the  evolution  of  the  streamwise  (the 
direction  of  compression)  component  of  kinetic  energy.  Flows  with  ui  at  lower  ki  have 
a  reduced  effect  of  pressure  during  compression  and  hence  a  higher  amplification  of  ui. 
The  evolution  of  upon  compression  is  influenced  by  the  initial  fraction  of  kinetic  energy 
in  the  direction  of  compression  (ui/?^),  in  addition  to  the  initial  jEJii(«i).  Flows  with  a 
larger  value  of  initial  u\j(^  and  ui  at  lower  ki  have  a  larger  amplification  of  ^ . 

Upon  normal  compression,  all  components  of  turbulent  kinetic  energy  of  sheared  tur¬ 
bulence  are  amplified  with  the  streamwise  component  being  amplified  the  most.  The 
amplification  of  uf  and  is  higher  than  that  in  isotropic  turbulence.  Normal  compression 
decreases  the  turbulent  shear  stress  and  for  large  enough  compressions  changes  its  sign. 
Examination  of  the  terms  in  the  shear  stress  evolution  equation  show  that  amplification 
of  the  pressure-strain  correlation  upon  compression  is  responsible  for  this  behavior. 
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The  oblique  angle  between  the  directions  of  shear  and  compression  is  seen  to  signifi¬ 
cantly  affect  the  response  of  sheared  turbulence  to  compression.  The  effect  of  oblique  angle 
on  the  evolution  of  kinetic  energy  may  be  explained  by  the  initial  distribution  of  Eee{Ke) 
and  ug^  jq^  ^  where  6  is  the  direction  of  compression.  Over  a  range  of  oblique  angles  vary¬ 
ing  from  —60°  to  60°  the  amplification  of  u\  and  is  seen  to  decrease  with  increasing 
magnitude  of  oblique  angle.  Also  oblique  compression  reduces  the  tendency  of  the  shear 
stress  to  decrease  in  magnitude;  for  large  oblique  angles  the  shear  stress  amplifies. 

With  respect  to  shock/turbulence  interaction,  our  results  suggest  that  besides  shock 
strength  (defined  in  terms  of  the  normal  Mach  munber),  the  anisotropy  of  the  turbulence 
and  the  shock  inclination  angle  are  important  parameters  in  determining  the  evolution  of 
turbulence  across  the  shock.  It  is  striking  that  the  evolution  of  the  Reynolds  stresses  upon 
oblique  compression  is  qualitatively  different  for  different  oblique  angles.  These  important 
effects  of  initial  anisotropy  and  the  shock  inclination  angle  should  be  accoimted  for  in  the 
interpretation  of  experiments  on  shock/turbulence  interaction. 
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Isotropic 


Shear 


Shear  +  Compression 


FIG.  2.1.  Schematic  of  the  mean  velocity  field  in  the  normal  compression  of  a  turbulent  shear 
fiow. 
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FIG.  2.2.  Evolution  of  Rij/q^  as  predicted  by  RDT  applied  to  the  homogeneous  shear  of  initially 
isotropic  turbulence. - (Rii/g^), - (1^22/?^), .  (Rss/q^), - {Rii/q^)- 
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FIG.  2.3.  Evolution  of  upon  the  normal  compression  of  sheared  turbulence.  The  diiferent 
curves  correspond  to  different  values  of  initial  total  shear.  - (^o  =  0), - = 

1), .  (/5o  =  2), - (^0  =  3). 


FIG.  2.4.  Evolution  of  the  components  of  turbulent  kinetic  energy  when  sheared  turbulence 
(/3o  =  3),  is  subjected  to  normal  compression. - (-Rii), - (-^22), .  (.^33)- 
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FIG.  2.5.  Comparison  of  the  amplification  of  of  sheared  turbulence  (^o  =  3)  to  that  of 
isotropic  turbulence.  - (Isotropic), - (Sheared). 


FIG.  2.6.  Evolution  of  Rijfq^  when  sheared  turbtilence  (^o  =  3)  is  subjected  to  normal  com- 

{R22/q% .  {Rzz/q^). 


pression. 


FIG.  2.7.  Evolution  of  when  sheared  turbulence  is  subjected  to  normal  compression.  The 

different  curves  correspond  to  different  values  of  initial  total  shear.  - (/?o  =  1), 

- - (^0=2), .  (/3o=3). 


FIG.  2.8.  Budget  of  terms  in  the  R12  evolution  equation  during  the  normal  compression  of  a 

shear  flow  (^0  =  3).  -  (LHS), - (Compression  production),  .  (Shear 

production), - (Pressure-strain  correlation). 
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FIG.  2.9.  Evolution  of  the  ratio  of  ttu  to  the  production  term  in  the  Ru  equation  upon  normal 
compression.  The  compression  of  sheared  turbulence  {/3o  =  3)  is  compared  to  that  of 
isotropic  turbulence. - (Isotropic), - (Sheared), 
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FIG.  2.10.  Coordinate  system  used  in  the  analysis  of  oblique  compression  of  sheared  turbulence 
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FIG.  2.11.  Evolution  of  jRu  upon  the  oblique  compression  of  sheared  turbulence  (/?o  =  3).  The 

different  curves  correspond  to  different  values  of  oblique  angle.  -  (0  —  0°), 

- (e  =  30°), .  (6  =  -30°), - (0  =  60°), - (0  =  -60°). 


pIpo 


FIG.  2.12.  Evolution  of  upon  the  oblique  compression  of  sheared  turbulence  (^o  =  3)-  The 

different  curves  correspond  to  different  values  of  oblique  angle.  -  {6  =  0°), 

- {6  =  30°), .  {9  =  -30°), - {$  =  60°), - {9  =  -60°). 
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FIG.  2.13.  Evolution  of  upon  the  oblique  compression  of  sheared  turbulence  (ydo  =  3). 


FIG.  2.14.  Evolution  of  R22fT  upon  the  oblique  compression  of  sheared  turbulence  (^o  =  3). 

The  different  curves  correspond  to  different  values  of  oblique  angle.  - {6  =  0°), 

- (0  =  30°), .  (0  =  -30°), - (B  =  60°), - (0  =  -60°). 


p/po 


FIG.  2.15.  Evolution  of  Rzzfq^  upon  the  oblique  compression  of  sheared  turbulence  (^o  =  3)- 

The  different  curves  correspond  to  different  values  of  oblique  angle.  - {6  =  0°), 

- {d  =  30°), .  {$  =  -30°), - {e  =  60°), - {e  =  -60°). 


FIG.  2.16.  Evolution  of  Ri2lq^  upon  the  oblique  compression  of  sheared  turbulence  (/5o  =  3). 

The  different  curves  correspond  to  different  values  of  oblique  angle.  - {6  =  0°), 

- {9  =  30°), .  {9  =  -30°), - {9  =  60°), - {9  =  -60°). 


43 


Compression-comer  flow  (0  >  0) 


Direction  of  shear 

i 


Direction  of  compression 


Shock  wave  incident  on  flat  plate 
boundary  layer  (0  <  0) 


FIG.  2.17.  Examples  of  the  oblique  compression  of  a  shear  flow. 


FIG.  2.18.  Comparison  of  model  solution  to  RDT  prediction  of  Ru.  -  (Isotropic;  RDT), 

- (Isotropic;  model), .  (Sheared;  RDT), - (Sheared;  model). 
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The  combined  interaction  of  vorticity  and 
entropy  fluctuations  with  a  shock  wave 
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Inviscid  linear  analysis  is  used  to  study  the  simultaneous  interaction  of  a  shock  wave  with 
vorticity  and  entropy  fluctuations.  The  two-dimensional  interaction  of  the  shock  with  a 
plane  vorticity-entropy  wave  is  initially  considered.  The  analysis  is  compared  to  numerical 
computations  and  then  extended  to  examine  the  evolution  of  an  isotropic  field  of  vorticity 
and  entropy  fluctuations  across  the  shock.  Expressions  are  derived  for  the  flow-field  behind 
the  shock  and  the  distortion  of  the  shock  front. 

The  amplitude  and  nature  of  correlation  between  vorticity  and  entropy  fluctuations  are 
shown  to  significantly  affect  the  level  of  kinetic  energy,  sound  and  thermodynamic  fluctua¬ 
tions  behind  the  shock  wave.  Entropy  fluctuations  in  phase  with  the  vorticity  fluctuations 
{u'  and  T'  negatively  correlated)  enhance  the  amplification  of  kinetic  energy,  vorticity  and 
thermodynamic  fluctuations.  Positively  correlated  fluctuations  have  the  opposite  effect. 
The  distortion  of  the  shock  front  is  similarly  affected.  These  trends  are  explained  and  a 
scaling  proposed  for  the  evolution  of  vorticity  fluctuations  across  the  shock. 

The  applicability  of  Morkovin’s  hypothesis  across  a  shock  wave  is  examined.  Statistics 
in  the  far-field  indicate  that  neglect  of  the  acoustic  mode  is  a  good  approximation  in  the 
far-field  of  shock  waves  of  moderate  strength  (Mi  <  2).  The  part  of  the  hypothesis  relating 
u'  and  T'  is  seen  to  be  invalid  behind  the  shock.  The  deviation  is  shown  to  be  a  result  of 
oscillation  of  the  shock  front. 

The  results  indicate  the  considerable  importance  of  entropy  fluctuations  in  the  evolution 
of  the  turbulent  fluctuations  in  a  compressible  boimdary  layer  upon  its  interaction  with  a 
shock  wave. 


1.  Introduction 

The  interaction  of  a  shock  wave  with  a  turbulent  boundary  layer  heis  received  consid¬ 
erable  attention  over  the  past  five  decades.  There  have  been  several  experimental  studies 
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of  the  flow  in  a  compression  corner,  normal  shock/boundary  layer  interaction  and  more 
recently,  the  interaction  of  shock  waves  with  three  dimensional  boundary  layers.  Green 
(1970)  and  Femholz  and  Finley  (1981)  provide  comprehensive  reviews  of  earlier  data.  More 
recent  work  is  discussed  by  Settles  and  Dodson  (1994)  and  Dolling  (1993).  These  exper¬ 
iments  testify  to  the  complex  flow-field  associated  with  the  interaction.  Incidence  of  the 
shock  wave  is  shown  to  strongly  affect  both  the  mean  flow  and  the  turbulent  fluctuations. 
The  Reynolds  stresses  and  temperatme  fluctuations  in  the  boundary  layer  are  seen  to  am¬ 
plify  across  the  shock  wave.  Significant  rmsteadiness  of  the  shock  wave  is  observed  when 
the  boundary  layer  separates.  The  motion  of  the  shock  wave  seems  to  be  correlated  with 
pressure  fluctuations  in  the  boimdary  layer  upstream  of  the  shock.  Recent  experiments  by 
Dolling  and  coworkers  (Dolling;  1993)  have  examined  the  fluctuating  pressure  loads  that 
shock-induced  separation  produces. 

Modeling  shock  wave/boimdary  layer  interaction  remains  a  challenge.  In  recognition  of 
the  complex  nature  of  the  phenomenon,  recent  studies  have  examined  the  simpler  and  more 
fundamental  problem  of  the  interaction  of  homogeneous  turbulence  with  a  shock.  These 
studies  which  include  experiments,  computation  and  analysis  observe  the  turbulence  to  be 
affected  in  a  manner  similar  to  boundary  layer  turbulence.  The  experiments  of  Keller  amd 
Merzkirch  (1990)  and  Honkan  and  Andreopoulous  (1990)  report  amplification  of  turbulent 
kinetic  energy  of  grid-generated  txirbulence  across  a  shock  wave.  The  simulations  of  Lee, 
Lele  and  Moin  (1993,1994)  and  Hannapel  and  Friedrich  (1994)  report  similar  enhance¬ 
ment  in  vorticity,  kinetic  energy  and  thermodynamic  fluctuations.  A  point  of  contention 
exists  between  experiment  and  computation;  experiments  report  an  increeuse  in  turbulence 
lengthscale  across  the  shock  while  computations  show  a  drop.  As  noted  by  Lele  (1993) 
this  issue  is  currently  imresolved. 

There  has  been  renewed  interest  in  the  linear  analysis  of  shock/turbulence  interaction. 
Issues  of  ‘shock-noise’  and  shock-stability  in  the  past  led  several  workers  to  examine  the 
interaction  of  a  shock  wave  with  a  weak  field  of  unsteady  disturbances.  A  brief  review 
of  linear  analysis  is  given  by  Mahesh,  Lee,  Lele  and  Moin  (1994b).  Ribner  (1953,  1954, 
1969,  1987)  appears  to  have  been  the  first  to  use  linear  analysis  to  predict  the  interac¬ 
tion  of  isotropic  vortical  turbulence  with  a  shock  wave.  Detailed  comparison  by  Lee,  Lele 
and  Moin  (1994)  show  good  agreement  between  their  computations  and  Ribner’s  analysis. 
Moore  (1954)  and  Chang  (1957)  have  performed  detailed  analysis  of  the  interaction  of  a 
shock  wave  with  an  acoustic  wave  and  entropy  wave  respectively.  The  influence  of  the 
dilatational  component  on  the  evolution  of  a  turbulent  flow  across  a  shock  was  examined 
by  Mahesh,  Lee,  Lele  and  Moin  (1994b)  who  extended  Moore’s  analysis  to  consider  the 
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interaction  of  an  isotropic  field  of  acoustic  waves  with  a  shock  wave.  A  drop  in  kinetic  en¬ 
ergy  was  observed  for  Mach  numbers  between  1.25  and  1.8.  For  Mach  numbers  beyond  3, 
the  kinetic  energy  was  found  to  amplify  by  levels  significantly  exceeding  those  in  Ribner’s 
problem.  Anisotropy  of  the  turbulence  and  the  obliquity  of  the  shock  wave  were  identi¬ 
fied  as  important  parameters  in  shock  wave/boundary  layer  interaction  by  the  analysis  of 
Mahesh,  Lele  and  Moin  (1994a). 

The  importance  of  entropy  fiuctuations  in  shock  wave/boundary  layer  interaction  seems 
to  be  underappreciated.  Morkovin  (1961)  hypothesized  that  the  turbulent  fluctuations 
in  compressible  boundary  layers  are  essentially  composed  of  vorticity  and  entropy  modes; 
i.e.  the  acoustic  mode  is  negligible.  This  hypothesis  along  with  the  asstunption  of  neg¬ 
ligible  fluctuations  in  stagnation  temperature  leads  to  the  following  relation  between  the 
thermodynamic  and  velocity  field, 


P 


TL 

T 


=  (7-1)M2^ 


(1) 


A  considerable  body  of  experimental  evidence  supports  the  above  hypothesis.  In  fact, 
Bradshaw  (1977)  states  ‘There  is  little  point  in  taking  turbulence  measurements  in  constant 
pressure  boundary  layers  at  M  <  5  to  check  Morkovin’s  hypothesis,  except  possibly  for 
intermittency  measurements’.  Morkovin’s  hypothesis  suggests  that  in  addition  to  being 
correlated,  the  intensity  of  vorticity  and  entropy  fluctuations  in  a  turbulent  boundary 
layer  are  of  comparable  magnitude.  Experimental  data  (Fernholz  and  Finley;  1981)  shows 
p'/p  in  a  constant  pressiue  boimdary  layer  to  have  a  maximum  value  of  about  5%  —  10%. 

This  suggests  that  entropy  fluctuations  would  play  an  important  role  in  the  interaction 
of  a  boundary  layer  with  a  shock  wave.  Study  of  the  influence  of  entropy  fluctuations  on 
shock/turbulence  interaction  seems  to  have  been  restricted  to  the  analysis  of  plane  entropy 
waves  interacting  with  a  shock  (Morkovin;  1960,  Chang;  1957,  Cuadra;  1968).  Hesselink 
and  Sturtevant  (1988)  have  conducted  experiments  on  the  passage  of  a  shock  wave  through 
a  random  field  of  entropy  fluctuations.  Their  interest  however  was  in  the  problem  of  the 
sonic  boom;  the  strongest  shock  they  considered  was  Mach  1.1.  The  evolution  of  the  shock 
front  in  this  field  of  entropy  fluctations  was  examined  in  their  study.  A  related  problem 
that  several  workers  have  examined  is  the  interaction  of  a  shock  wave  with  a  density 
inhomogeneity.  A  list  of  references  may  be  found  in  a  recent  paper  by  Yang,  Kubota  and 
Zukoski  (1994). 

In  this  paper  we  seek  to  identify  the  quantitative  influence  that  entropy  fluctuations  exert 
on  the  evolution  of  a  turbulent  flow  across  a  shock  wave.  We  use  inviscid  linear  analysis 
to  study  the  combined  interaction  of  vorticity  and  entropy  fiuctuations  with  a  shock  wave. 
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To  the  best  of  our  knowledge,  the  interaction  of  a  field  of  entropy  fiuctuations  with  a  shock 
wave  has  not  been  considered  before.  Nor  does  effect  of  the  correlation  between  vorticity 
and  entropy  fluctuations  appear  to  be  emphasized. 

This  paper  is  organized  as  follows.  The  simultaneous  interaction  of  a  shock  wave  with  a 
vorticity  and  entropy  wave  (henceforth  referred  to  as  vorticity-entropy  wave)  is  formulated 
in  Section  2.  Results  of  the  linear  analysis  are  presented  in  Section  3  and  compared  in 
Section  4  to  nmnerical  computation.  Sections  5  and  6  provide  an  explanation  for  the 
influence  of  entropy  fiuctuations  and  suggest  a  scaling  for  the  evolution  of  vorticity  across 
the  shock.  The  analysis  is  extended  in  Section  7  to  consider  the  evolution  of  an  isotropic 
field  of  vorticity-entropy  waves  across  a  shock  wave.  Section  8  discusses  the  results.  The 
paper  is  concluded  with  a  summary  of  the  important  conclusions  in  Section  9.  A  detailed 
version  of  the  paper  is  presented  in  a  report  by  Mahesh,  Lele  and  Moin  (1994c). 

2.  The  interaction  of  a  vorticity-entropy  wave  with  a  shock  wave 

2.1  The  three  modes  in  a  compressible  medium 

Kovasznay  (1953)  showed  that  fluctuations  of  small  amphtude  in  a  compressible  medium 
could  be  decomposed  into  vorticity,  acoustic  and  entropy  modes.  These  modes  were  shown 
to  evolve  independently  in  the  inviscid  liixxit  for  uniform  mean  flow.  The  vorticity  mode  as 
defined  by  Kovasznay  has  no  pressure  or  density  fluctuations.  It  has  a  solenoidal  velocity 
field  that  is  convected  by  the  mean  flow.  The  acoustic  mode  travels  at  the  speed  of 
sound  relative  to  the  mean  flow,  has  isentropic  pressiire  and  density  fluctuations  and  a 
corresponding  irrotational  velocity  field  that  satisfies  the  acoustic  wave  equation.  The 
entropy  mode  is  convected  by  the  mean  flow  and  has  no  velocity  or  pressure  fluctuations; 
it  has  only  density  and  temperature  fluctuations. 

The  independent  evolution  of  these  modes  in  the  inviscid  hmit  implies  that  the  fluctu¬ 
ating  velocity  field  is  a  hneax  superposition  of  acoustic  and  vortical  modes.  Similarly,  the 
density  and  temperature  fields  are  superpositions  of  acoustic  and  entropic  components. 
The  pressiire  field  is  associated  solely  with  the  acoustic  mode. 

2.2  Formulation  of  the  problem 

The  two-dimensional  interaction  of  a  shock  wave  with  a  plane  vorticity-entropy  wave  is 
schematically  illustrated  in  Figure  1.  The  two-dimensional  problem  will  be  extended  to 
describe  the  interaction  of  the  shock  with  a  three-dimensional  field  of  turbulence  in  Section 
7.  Note  that  the  shock  wave  is  stationary  in  the  mean.  The  variables  U,P,'p,T  and  M 
denote  the  mean  velocity,  pressure,  density,  temperature  and  Mach  number  respectively 
and  subscripts  1  and  2  denote  the  upstream  and  downstream  states.  The  flow  upstream 
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of  the  shock  wave  is  pertiirbed  by  the  weak  disturbance  field  of  the  incident  vorticity- 
entropy  wave  which  is  assumed  to  be  a  plane  wave  that  makes  angle  V’l  with  the  x  eixis. 
The  variables  u,v,p,p  and  T  represent  the  fluctuating  velocities,  pressure,  density  and 
temperature.  The  incident  field  has  the  following  form: 


“l 

Ui 

_  ^ik(Tnx+ly-Uimt) 

(2a) 

Ui 

—  ^ik(mx+ly-Uimt) 

(2b) 

P‘i 

_  ^ik(mx+ly-Uimt) 

(2c) 

Pi 

T[ 

= 

(2d) 

Ti 

Pi 

Pi 

=  0 

(2e) 

where  the  variables  m  =  cos^’i  and  I  =  sinV*!.  The  shock  wave  deforms  in  response  to  the 
incident  disturbance;  the  displacement  is  denoted  by  x  =  ^(y,  t). 

The  analysis  involves  solution  of  a  boundary- value  problem  for  the  shock  displacement 
and  the  flow  field  behind  the  shock  wave.  The  governing  equations  behind  the  shock 
wave  are  the  Euler  equations  that  are  Hnearized  about  the  uniform  flow  behind  the  shock. 
Boundary  conditions  are  imposed  at  r  =  0  (the  mean  position  of  the  shock  wave)  by 
hnearizing  the  Rankine  Hugoniot  equations  in  a  frame  of  reference  that  moves  at  the 
instantaneous  speed  of  the  shock  wave.  The  equations  are  solved  in  consonance  with  the 
boundary  conditions  to  yield  the  flow  behind  the  shock  wave  and  the  distortion  of  the 
shock  front. 

The  solution  could  be  obtained  by  superposing  the  solutions  to  the  interaction  of  the 
shock  with  a  vorticity  and  entropy  wave  independently.  Rather  than  use  results  of  different 
workers  (Ribner;  1953,  Chang;  1957),  we  solve  the  above  set  of  equations.  Comparison 
will  reveal  our  method  of  solution  to  be  much  simpler.  In  the  process,  the  independent 
interaction  of  the  shock  with  vorticity  and  entropy  waves  become  special  cases  of  om 
analysis  and  will  be  used  as  baseline  ceises  for  comparison.  Our  solution  in  these  limiting 
cases  is  of  course  identical  to  that  of  the  above  workers. 

2.S  The  solution  behind  the  shock  wave 

Details  of  the  analysis  are  presented  in  Appendix  A.  The  solution  behind  the  shock 
wave  has  two  different  regimes  that  differ  in  the  nature  of  the  pressure  field.  The  incidence 
angle  of  the  disturbance  tpi  and  the  Mach  number  of  the  shock  wave  Mi  determine  the 
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regime  of  the  solution.  For  a  given  Mach  number,  the  two  regimes  are  demarcated  by  a 
critical  angle  of  incidence  denoted  by  ^c-  If  0  <  V"!  <  i’c  the  pressure  field  behind  the 
shock  wave  (and  hence  its  associated  velocity  and  thermodynamic  field)  is  a  plane  wave. 
If  however,  <  V’l  <  '^/25  the  pressure  field  behind  the  shock  wave  corresponds  to  an 
evanescent  wave  and  decays  exponentially.  The  vorticity  and  entropy  fields  behind  the 
shock  wave  correspond  to  plane  waves  over  both  regimes.  The  critical  angle  xfc  is  a  root 
of  the  following  equation; 


Note  that  is  a  function  of  the  Mach  number  of  the  shock  wave.  The  existence  of  two 
regimes  of  the  pressure  field  is  a  common  featiure  of  problems  involving  acoustic  wave 
propagation.  Classical  examples  include  the  scattering  of  sound  by  a  vortex  sheet  and  the 
propagation  of  sound  generated  by  a  flexural  wall. 

The  solution  behind  the  shock  wave  may  be  represented  as  a  superposition  of  vortical, 
acoustic  and  entropic  components.  It  has  the  following  dimensionless  form: 
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(4) 


The  boundary  conditions  across  the  shock  wave  yield  the  following  expressions  for  the 
velocity  and  slope  of  the  shock  front: 


_L.^  =  Ij  e**  (ly-mUit).  _ _ ^ik  (ly-mUit) 

Av  Ui  ’  Av  m 


The  variable  r  denotes  the  mean  density  ratio  P2/P1  across  the  shock  wave.  As  indicated  in 
Appendix  A,  k{Mi ,  V’l )  is  the  streamwise  waveniimber  of  the  acoustic  component  behind 
the  shock  wave.  It  is  real  in  the  propagating  regime  and  complex  in  the  decaying  regime. 
The  coefficients  F,H  and  K  are  associated  with  the  acoustic  component.  The  vortical 
component  is  represented  by  G  and  I  while  Q  represents  the  entropy  component.  The 
coefficients  F,H,  K,G,  I  and  Q  are  functions  of  Mi,  ^1  and  the  amplitude  ratio  and 
phase  difference  between  the  vorticity  and  entropy  waves  {Ae/Ay). 
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3.  Results 

Expressions  are  now  derived  to  describe  the  streamwise  evolution  of  statistics  behind 
the  shock  wave.  Consider  for  example  the  streamwise  component  of  velocity.  The  kinetic 
energy  associated  with  this  component  is  denoted  by  U2^  =  where  the  overbar 

implies  averaging  over  the  transverse  direction  and  time  and  the  asterisk  denotes  the 
complex  conjugate.  The  following  expression  is  obtained  for  from  equation  (4): 
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/2 


ui 


=  [  \F\^  +  |G|2  +  ]  |A„| 
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The  spatial  variation  of  is  different  over  the  two  regimes.  The  wavenumber  k  is  real 
over  the  propagating  regime  and  hence 


^  ^  ........ 

=  [  |F|2  +  |Gp  +  2(FrGr  +  FiGi)cos{k  -  kmr)x 

—  2{FiGr  —  FrGi)sm{k  —  kmr)x]  (7) 


where  the  subscripts^’  and  H'  denote  the  real  and  imaginary  parts  respectively.  The  kinetic 
energy  thus  has  spatially  uniform  contributions  from  the  vortical  and  acoustic  modes  of 
the  velocity  field  and  an  oscillating  component  whose  argument  is  the  phase  difference 
between  them.  Similar  expressions  may  be  written  for  quantities  such  as  density  and 
temperature  in  the  propagating  regime.  Variables  such  as  vorticity, dilatation  and  pressure 
that  depend  on  only  one  mode  are  of  course  spatially  uniform,  e^., 

^  =  I^P  M  (8) 

The  evolution  of  over  the  decaying  regime  may  be  shown  to  be  given  by: 

7^2  ^  ^  „ _ 
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The  corresponding  equation  for  pressure  is: 


g-2A:.i 


\K?  KP 


(10) 


The  kinetic  energy  over  this  regime  is  seen  to  have  a  spatially  uniform  vortical  component, 
an  acoustic  component  that  decays  exponentially  behind  the  shock  wave  and  a  damped 
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oscillating  component  due  to  the  correlation  between  the  acoustic  and  vortical  modes  that 
the  shock  wave  introduces.  In  describing  the  flow-field  behind  the  shock  wave,  a:  =  0  is 
referred  to  as  the  ‘near-field’  and  x  oo  is  called  the  ‘far-field’.  We  will  see  that  the 
fax-field  values  are  attained  over  a  distance  comparable  to  the  lengthscale  of  the  incident 
disturbance.  As  a  result,  it  is  the  far-field  values  that  are  of  practical  importance. 

The  analysis  is  now  applied  to  the  interaction  of  a  vorticity-entropy  wave  with  a  Mach 
1.5  shock  wave.  Incidence  angles  of  45°  and  75°  are  considered  so  that  behavior  over  both 
regimes  is  illustrated.  (The  critical  angle  is  61.36°).  The  objective  here  is  to  demonstrate 
the  influence  of  entropy  fluctuations  and  their  correlation  with  the  incident  vorticity  fluc¬ 
tuations.  The  interaction  of  the  vorticity-entropy  waves  is  therefore  contrasted  with  the 
independent  interaction  of  the  shock  with  a  vorticity  wave  and  entropy  wave  at  the  same 
angle  of  incidence.  These  two  limiting  cases  are  referred  to  as  the  ‘pure  vorticity’  and  ‘pure 
entropy’  limits. 

The  normalization  of  the  curves  in  the  pure  entropy  limit  deserves  clarification.  There 
is  no  incident  velocity  field  in  the  pure  entropy  limit;  i.e.,  q\  =  =0.  The  appropri¬ 

ate  normalization  in  this  limit  is  (92/^?)/(pi^/Pi^)-  We  set  the  intensity  pj/pi  in  the 
pure  entropy  limit  to  equal  its  value  in  the  vorticity-entropy  wave;  i.e.,  (pi/pj)pure  ent.  = 
vort.— ent.*  This  allows  a  consistent  comparison  between  the  pure  entropy  limit 
and  the  other  curves. 

Figure  2  shows  the  evolution  of  kinetic  energy  behind  the  shock  wave  over  the  propagat¬ 
ing  regime  (V’l  =  45°).  All  three  curves  are  normalized  by  the  incident  kinetic  energy.  The 
presence  of  entropy  fluctuations  is  seen  to  significantly  affect  the  amplification  of  kinetic 
energy  across  the  shock  wave.  The  evolution  of  is  strongly  affected  by  both  the  am¬ 
plitude  ratio  and  phase  difference  between  the  vorticity  and  entropy  waves.  Dependence 
upon  the  amplitude  is  monotonic.  Entropy  fluctuations  that  are  in  phase  with  the  vortic¬ 
ity  fluctuations  and  T[  are  negatively  correlated)  enhance  the  amplification  of  kinetic 
energy.  Suppression  of  amplification  is  observed  when  u\  and  T[  are  positively  correlated. 
The  pressure  fluctuations  behind  the  shock  wave  are  similarly  affected. 

The  decaying  regime  (V^i  =  75°)  as  shown  in  Figure  3  displays  similar  influence  of  the 
entropy  fluctuations.  Both  near-field  and  fax-field  levels  of  kinetic  energy  axe  affected. 
The  dependence  upon  amphtude  ratio  and  phase  difference  is  identical  to  that  seen  in  the 
propagating  regime.  Also  note  that  the  curves  asymptote  to  the  fax-field  levels  within  a 
wavelength  of  the  incident  disturbance. 

This  influence  of  the  entropy  fluctuations  is  not  limited  to  the  two  angles  considered 
above;  it  extends  over  all  angles  of  incidence.  This  is  illustrated  in  Figure  4  where  the 
fax-field  levels  of  amplification  of  axe  plotted  against  the  angle  of  incidence.  The  Mach 
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number  of  the  shock  wave  is  1.5.  Note  that  the  notion  of  far-field  is  introduced  into  the 
propagating  regime  by  ignoring  the  oscillating  component  in  equation  (9).  Reason  for 
this  neglect  is  statistical  cancellation  of  this  component  in  the  far-field  for  an  incident 
turbulent  field.  Also  shown  in  Figure  5  is  the  amplification  of  vorticity  across  the  shock 
wave.  Negative  correlation  between  u'  and  T'  is  clearly  seen  to  increase  the  amplification 
while  positive  correlation  suppresses  it. 

4.  Comparison  to  numerical  computation 

The  analysis  is  evaluated  by  comparison  to  numerical  solution  of  the  interaction  of  a 
sinusoidal  vorticity-entropy  wave  with  a  Mach  1.5  shock  wave.  The  governing  equations  are 
the  two-dimensional  unsteady  compressible  Navier  Stokes  equations  which  are  discretized 
using  the  sixth  order  Fade  scheme  (Lele;  1992)  and  the  third  order  Runge  Kutta  scheme 
(Wray;  1986).  Details  of  the  computation  are  discussed  by  Mahesh,  Lele  and  Moin  (1994c). 
The  shock  wave  is  stationary  at  the  center  of  the  domain  in  the  mean.  The  shock  structure 
is  a  result  of  molecular  viscosity;  no  shock  captmring  or  shock  -fitting  is  used.  A  structured 
mesh  that  is  non-uniform  in  the  streamwise  direction  is  used  to  resolve  the  shock  wave. 
About  7  mesh  points  are  located  inside  the  shock  wave. 

Periodic  boundary  conditions  are  imposed  transverse  to  the  shock.  The  disturbance  field 
corresponding  to  a  vorticity-entropy  wave  is  superposed  onto  the  supersonic  mean  fiow  at 
the  inflow  boundary.  Approximately  non-reflecting  boundary  conditions  (Poinsot  and  Lele; 
1990)  are  used  at  the  outflow  boundary.  The  initial  condition  is  a  numerically  computed 
steady  shock  wave.  The  vorticity-entropy  wave  is  then  introduced  through  the  inflow 
boundary  condition.  The  transients  set  up  by  the  initial  period  of  adjustment  are  allowed 
to  exit  the  domain.  Statistics  are  then  gathered  over  a  period  of  the  inflow  disturbance. 

The  computations  examine  the  dependence  of  the  interaction  upon  the  incidence  angle  of 
the  disturbance.  The  incidence  angle  was  vgiried  from  0  to  90  degrees  while  the  amplitudes 
At,  and  Ag  were  fixed  at  0.05.  We  compare  in  Figure  6  the  computed  values  of  amplification 
of  vorticity  fluctuations  to  the  linear  anedysis  predictions.  Good  agreement  is  seen  away 
from  the  critical  angles  inspite  of  the  computations  solving  the  viscous  equations.  Deviation 
is  seen  around  the  critical  angle.  This  behavior  is  consistent  with  the  interaction  of  a  shock 
with  an  acoustic  wave  (  Zang  et  al;  1984,  Mahesh  et  al;  1994b). 

As  noted  by  Ribner  (1953)  and  Chang  (1954),  within  linear  analysis  the  unsteady  in¬ 
teraction  of  an  oblique  wave  with  a  normal  shock  may  be  transformed  into  the  steady 
interaction  of  an  oblique  wave  with  an  oblique  shock.  The  transformation  involves  defin¬ 
ing  the  coordinates  Xg  =  x,ys  =  y—U\t  cot  .  The  governing  equations  in  the  transformed 
coordinates  are  the  steady  Euler  equations  linearized  about  uniform  mean  flow  at  a  Mach 
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number: 


_  Uj  +  cot^^^U^ 


as 


(11) 


Mg  equals  unity  at  the  critical  angle  whose  vicinity  therefore  corresponds  to  transonic  mean 
flow  in  transformed  coordinates.  We  believe  (Mahesh  et  al]  1994b)  that  the  inconsistency 
of  linear  equations  in  the  transonic  regime  causes  the  deviation  between  computation  and 
linear  analysis  around  the  critical  angle. 


5.  A  simple  explanation 

An  explanation  is  provided  for  the  influence  of  entropy  fluctuations  on  the  evolution  of 
a  turbulent  flow  across  a  shock  wave.  Consider  the  following  idealization  of  the  mean  fleld 
associated  with  the  shock: 


V  =  U(x);  P  =  P(xy,  ?  =  (12) 

Linearization  of  the  Euler  equations  about  the  above  mean  flow  yields  the  following  gov¬ 
erning  equations  for  the  fluctuations: 

p  {u[  +  Uu',  +  u't/*)  +  p'UU,  =  -p'x  (13) 

p  {vl  +  Uv’^)  =  -p'y  (14) 

Denoting  the  fluctuating  vorticity  by  w',  the  following  equation  may  be  derived  for  w': 

a;;  t/a;;  =  -u'U,  -  ^  Px  +  ^  P.  (15) 

The  above  set  of  equations  represents  the  effects  of  bulk  change  across  the  shock  wave. 
Effects  due  to  shock  distortion  are  absent.  The  motivation  to  examine  the  above  set 
of  equations  is  provided  by  the  fact  that  shock  distortion  effects  are  secondary  to  bulk 
effects  especially  in  the  propagating  regime.  Also,  though  the  interaction  of  a  vortical  field 
with  a  shock  wave  generates  acoustic  waves,  the  contribution  of  these  acoustic  waves  to 
the  downstream  kinetic  energy  is  not  significant;  amplification  and  generation  of  vortical 
fluctuations  accounts  for  most  of  the  kinetic  energy  behind  the  shock  wave. 

The  evolution  of  the  vorticity  fluctuations  can  therefore  be  related  to  that  of  kinetic 
energy.  Equation  (15)  shows  three  terms  that  modify  the  vorticity  field.  The  first  term 
{—ui'Ux)  represents  the  effect  of  bulk  compression.  The  drop  in  mean  velocity  across  the 
shock  wave  indicates  that  this  term  would  enhance  vorticity  fluctuations.  The  second  and 
third  terms  in  equation  (15)  represent  baroclinic  contributions  to  the  change  in  vorticity 
across  the  shock  wave. 
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If  the  incident  disturbance  comprises  of  vorticity  and  entropy  fluctuations,  then  p'  =  0 
in  the  incident  fleld  and  hence, 

cu't  +  UJ,  =  -u;'  (16) 

The  incident  vorticity  fluctuations  are  thus  enhanced  by  bulk  compression.  The  incident 
entropy  fluctuations  produce  vorticity  at  the  shock  wave  through  the  baroclinic  term.  The 
baroclinic  contribution  can  enhance  or  oppose  the  effect  of  bulk  compression.  The  phase 
difference  between  the  vorticity  and  entropy  waves  determines  whether  enhancement  or 
opposition  is  observed. 

Consider  for  example  the  plane  vorticity-entropy  wave  represented  by  equation  (2).  It 
is  easily  shown  that 

-u^'U,  A,  UU,  -  Ae  (17) 

P  P 

Since  Ux  is  negative  and  Px  is  positive  across  a  shock  wave,  the  two  sources  of  vorticity  are 
of  the  same  sign  if  Ae  and  Av  are  of  the  same  sign.  They  oppose  each  other  if  Ae  and  A^, 
are  of  opposite  sign.  Thus,  if  u'  and  T  are  negatively  correlated  the  entropy  field  enhances 
the  amplification  of  fluctuating  vorticity.  On  the  other  hand,  a  positive  correlation  between 
u'  and  T'  suppresses  the  amplification  of  vorticity  across  the  shock  wave. 

Further  insight  is  gained  from  a  schematic  illustration  of  this  effect.  Figure  7  shows  a 
fluid  element  of  circular  cross-section  passing  through  a  shock  wave.  The  geometric  center 
of  this  element  is  denoted  by  Cp  while  the  variable  Cm  is  used  to  represent  the  center 
of  mass.  The  disturbance  field  associated  with  the  fluid  element  is  that  of  a  vorticity- 
entropy  wave.  The  element  therefore  exhibits  solid  body  rotation  with  associated  vorticity 
u'  which  is  assumed  positive  in  the  direction  shown.  Also,  the  density  field  associated  with 
the  entropy  wave  causes  the  centre  of  mass  to  differ  from  the  centre  of  force  (the  geometric 
centre). 

Bulk  compression  compresses  the  element  in  the  streamwise  direction  thereby  enhancing 
the  rotation.  In  addition,  the  shock  wave  exerts  a  pressure  force  (associated  with  the 
adverse  pressure  gradient)  on  the  fluid  element  passing  through  the  geometric  centre  as 
shown.  The  non-coincidence  of  and  Cm  would  cause  this  pressure  force  to  exert  a 
torque  about  the  centre  of  mass.  It  is  this  torque  that  manifests  itself  as  the  baroclinic 
source  of  vorticity.  Note  that  if  Cm  is  above  Cp  the  baroclinic  rotation  is  in  the  same 
direction  as  u'.  The  reverse  is  true  if  Cm  is  below  Cp.  The  correlation  between  u'  and  T' 
over  the  fluid  element  determines  the  location  of  Cm  with  respect  to  Cp  and  thereby  the 
relative  sense  of  rotation  that  baroclinic  torque  produces. 
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6.  Scaling  of  the  evolution  of  vorticity  across  a  shock  wave 

The  simplicity  of  equation  (16)  motivates  its  use  to  derive  formulae  that  approximate  the 
evolution  of  vorticity  fluctuations  across  a  shock  wave.  The  performance  of  the  formulae 
may  be  evaluated  by  comparison  to  results  of  the  hnear  analysis.  Equation  (16)  may  be 
rewritten  as: 

Dt{UJ)  =  UP^  (18) 

where  Dt  denotes  the  material  derivative  d/dt  +  Ud/dx.  Use  of  the  relation  Px  =  —'pUUx 
in  the  above  equation  yields: 

D,(Vu>')  =  ^  ^  (V%  (19) 

Approximation  of  the  shock  wave  as  a  discontinuity  allows  {U^)x  to  be  expressed  as 
AU^S(x)  where  A(U^)  represents  the  difference  in  across  the  shock  wave.  The  term  Py 
is  set  to  equal  its  upstream  value  to  allow  integration  of  the  above  equation.  Transforming 
coordinates  to  x'  =  x  —  Ut,  r  =  t  and  integrating  yields  the  following  expression  for  the 
change  in  vorticity  across  the  shock  wave: 


U2Uj'2  -  UiwJ 


m ,  (ui  -  ul) 

3  £7, 


(20) 


which  yields: 


^2  ~ 


,  ikl  . 


Ui 


1-r" 

J.2 


(21) 


The  above  equation  suggests  that  the  amplification  of  incident  vortical  fluctuations  equals 
the  mean  density  ratio  across  the  shock  wave  while  the  production  of  vorticity  by  entropy 
fluctuations  scales  as  kAel{l  —  r^)/r^. 

The  ability  of  these  formulae  to  predict  the  evolution  of  vorticity  is  evaluated  in  Figures  8 
and  9.  Figure  8  considers  the  amplification  of  vorticity  in  the  interaction  of  a  pure  vorticity 
wave  with  a  shock.  Mean  Mach  numbers  from  1.25  to  2.5  are  considered.  Both  scaled 
and  unsealed  values  of  rms  vorticity  behind  the  shock  wave  are  plotted  as  a  function  of 
incident  angle.  The  scaling  is  seen  to  work  qvdte  well  especially  in  the  propagating  regime. 
A  systematic  deviation  is  seen  with  increasing  Mach  number. 

The  interaction  of  a  pure  entropy  wave  with  a  shock  is  examined  in  Figure  9.  Once  again 
the  proposed  scaling  seems  to  yield  reasonable  collapse  of  the  curves.  The  agreement  in 
the  absence  of  effects  of  shock  distortion  is  quite  encouraging. 
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7.  Interaction  of  a  shock  wave  with  an  isotropic  field  of  vorticity-entropy  waves 

The  analysis  is  extended  to  describe  the  evolution  of  a  turbulent  flow  across  a  shock 
wave.  In  light  of  Morkovin’s  hypothesis,  the  turbulence  is  represented  as  a  random  three- 
dimensional  fleld  of  vorticity-entropy  waves.  For  reasons  of  generality,  the  turbulent  fleld 
is  assumed  to  be  isotropic.  It  will  be  seen  however,  that  the  procedure  is  easily  extended 
to  consider  anisotropic  turbulence. 

The  incident  tiurbulent  fleld  is  represented  as  a  superposition  of  plane  vorticity-entropy 
waves  (Fourier  modes)  in  three  dimensions.  Each  of  these  waves  would  interact  indepen¬ 
dently  with  the  shock  wave  under  linear  analysis.  For  a  given  upstream  spectrum,  the 
interaction  of  each  of  these  waves  with  the  shock  wave  is  predicted.  Integration  over  all 
waves  behind  the  shock  wave  yields  turbulence  statistics  behind  the  shock. 

The  three  dimensional  problem  is  related  to  the  two  dimensional  analysis  of  the  preceding 
section  as  follows.  Consider  an  incident  plane  wave  in  three  dimensions.  As  shown  in 
Figure  10,  the  wavenumber  vector  of  the  wave  lies  in  a  plane  that  makes  angle  4>  with  the 
y  axis.  In  this  plane  which  we  call  the  x  —  Xr  plane,  the  wave  makes  angle  'ipi  with  the  x 
axis.  It  is  readily  seen  that  the  x  —  Xr  plane  is  identical  to  the  plane  of  interaction  in  the 
two-dimensional  problem. 

The  solenoidal  nature  of  the  incident  velocity  fleld  requires  the  velocity  vector  of  the 
wave  to  be  normal  to  the  wavenumber  vector.  The  velocity  fleld  may  therefore  be  expressed 
as  a  sum  of  two  components:  one  normal  to  the  wavenumber  vector  in  the  x  —  Xr  plane 
and  the  other  normal  to  the  x  —  Xr  plane  (the  (j>  direction).  It  is  intuitively  clear  that  the 
<!>  component  of  velocity  would  pass  unchanged  through  the  shock  wave.  As  a  result,  the 
three  dimensional  problem  may  be  solved  using  results  of  the  two  dimensional  analysis  in 
the  X  —  Xr  plane.  The  subscripts  V’  emd  are  used  to  denote  the  components  in  the  Xr 
and  (f)  directions  respectively.  Equations  (4)  and  (6)  allow  the  following  expressions  to  be 
written  for  the  components  of  kinetic  energy  behind  the  shock  wave: 


Ell 

m 

Err 


=  [  |F|2  -h  |G|2  -F  ]  |A„|2 

=  [  \H\^  q-  |7p  -j- -t-  j  |^^|2 


(22) 

(23) 


Recall  that  the  variables  in  parentheses  depend  upon  Mi ,  xj^i ,  / Ay  and  k.  The  term  |A„| 

represents  the  magnitude  of  the  velocity  vector  in  the  x  —  r  plane.  Figure  10  shows  that 
\Ay  \  =  |mi|//  which  implies: 


P 


m 


(24) 
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where  the  superscript  ‘1’  denotes  the  state  upstream  of  the  shock  and  is  obtained 
from  the  following  expression  for  the  isotropic  energy  spectrum  tensor: 

p(i)  _  ( c 

where  E{k)  is  the  three  dimensional  energy  spectrum  tensor  such  that  f  E(k)  dk  =  q\l2 
and  the  wavenmnbers  in  cartesian  coordinates  are  given  by: 

ki  =  fccos^i;  k2  =  fesin^i  cos(f>;  ks  =  fesin^’i  sin^  (26) 


Ell  and  Err  behind  the  shock  wave  will  be  completely  determined  once  E(k)  and  A^/Av 
axe  specified.  We  use  the  following  form  for  E(k): 


E(k)  ~ 


(27) 


The  quantity  A^/Ay  may  be  represented  as  Arc'^''  where  Ar  and  (j>r  are  both  fimctions  of 
the  wavenumber  vector.  Appropriate  functional  dependencies  may  be  assumed  depending 
upon  the  flow  being  considered.  In  the  interests  of  generality,  this  paper  assiunes  Ar  and 
4)r  to  be  independent  of  the  wavenumber  vector.  Also,  the  results  presented  consider  the 
hmiting  cases  of  =  0  and  tt.  For  these  angles  it  is  easily  shown  that 


=  V2 


Aim 

VW/Pi 


(28) 


Once  E{}z)  and  AtjA^  are  specified,  equations  (22)  and  (23)  axe  used  to  obtain  En  and 
Err  behind  the  shock  wave.  Also  Fig.  9  may  be  useed  to  show  that: 


E22  +  Ezz  —  Err  +  E^fj,  (29) 

Recognizing  that  turbulence  behind  the  shock  wave  will  be  axisymmetric,  we  solve  for 
E22  +  Ezz-  E^^  was  noted  to  remain  unchanged  across  the  shock  wave  and  hence  is 
equated  to  its  upstream  value.  The  upstream  value  of  E^^  is  given  by: 

=  SS  -  IXP  =  ^  (30) 

The  spectra  axe  then  integrated  over  all  wavenumbers  to  obtain  statistics  as  a  function 
of  distance  behind  the  shock  wave.  The  elemental  volume  of  integration  is  given  by  the 


expression: 


d^k  =  sin^i  d^i  d<f>dk 


(31) 
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where  k  varies  from  0  to  oo,  tpi  varies  from  — 7r/2  to  7r/2  and  (f)  varies  from  0  to  27r. 

The  following  expressions  are  obtained  for  the  streamwise  variation  of  kinetic  energy 
behind  the  shock  wave: 


^l2  yoo  fir/2  ^oo  r7r/2 

~  I  I  I  -^11  sin^i  d<l>  dtpi  dk  =  An  I  I  Eu  k^  sin^i  d^i  dA; 

J  k=oJipi  =—-K/2J  4>=0  J  k=oJil:i=0 


(32) 

since  the  problem  is  symmetric  about  V’l  =  0  and  the  transfer  fimctions  that  yield  Eu 
are  independent  of  The  integration  over  k  and  is  performed  numerically.  Similar 
expressions  are  obtained  for  the  transverse  components  of  kinetic  energy,  i.e., 


k"^  sin  d^i  dA; 


(33) 


The  evaluation  of  scalar  quantities  behind  the  shock  wave  is  considerably  simpler.  For 
example  the  pressure  field  behind  the  shock  is  given  by  the  following  expression  (equations 
8  and  10)  : 

1 2  »oo  A’r/2 

=  An  /  Epp  k'^  sin  ipi  dV’i  dA;  (34) 

^2  Jk=0  Jrl>i=0 

where  Epp  =  \K\^\Av\'^  in  the  propagating  regime  and  e“^*‘*|ii'p|AvP  in  the  decaying 
regime  respectively. 


8.  Results 


8.1  The  spatial  evolution  of  kinetic  energy 

The  variation  of  with  distance  behind  a  Mach  1.5  shock  wave  is  plotted  in  Figure 
11.  The  combined  interaction  of  vorticity  and  entropy  fiuctuations  is  compared  to  that 
of  vorticity  and  entropy  fiuctuations  alone.  As  mentioned  before,  two  representative  cases 
of  combined  interaction  are  considered:  one  where  u'  and  T'  are  negatively  correlated 
upstream  of  the  shock  [Ar  =  0.58, 4>r  =  0)  and  the  other  where  u'  and  T'  are  positively 
correlated  {Ar  =  0.58,  =  tt).  The  value  of  0.58  was  (arbitrarily)  chosen  to  satisfy 

Morkovin’s  hypothesis  at  a  Mach  number  of  1.35.  All  four  curves  are  normalized  with  the 
incident  value  of  .  As  discussed  in  Section  3,  j q\  in  the  pure  entropy  limit  corresponds 
to  0.17  (?2/^i  )/(/=’! ^/Pi^)  where  the  factor  0.17  equals  {p\^ ! P\^)/{qllUl)  for  At  =  0.58. 

Similar  spatial  variation  is  exhibited  by  all  three  curves.  Amplification  of  across  the 
shock  wave  is  followed  by  rapid  drop  and  rise  immediately  downstream  of  the  shock  wave. 
Temperature  fiuctuations  that  are  negatively  correlated  with  u'  yield  higher  levels  of  kinetic 
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energy  at  all  streamwise  locations  while  positively  correlated  temperature  fluctuations  have 
the  opposite  effect. 

The  rapid  variation  of  kinetic  energy  behind  the  shock  wave  is  a  fundamental  feature 
of  shock/turbulence  interaction.  DNS  of  isotropic  vortical  fluctuations  interacting  with 
a  shock  wave  (Lee,  Lele  and  Moin,  1993)  exhibits  this  feature,  as  does  linear  analysis  of 
acoustic  fluctuations  (Mahesh,  Lee,  Lele  and  Moin,  1994b)  interacting  with  a  shock.  The 
discussion  (Section  3)  of  single  vorticity-entropy  waves  interacting  with  a  shock  shows  that 
waves  incident  at  angles  greater  than  the  critical  angle  (decaying  regime)  produce  this 
rapid  variation  through  the  damped  oscillatory  correlation  between  acoustic  and  vortical 
fluctuations  behind  the  shock  wave. 

This  correlation  was  overlooked  by  Lee,  Lele  and  Moin  (1993)  in  their  linear  analysis  of 
a  shock  wave  interacting  with  isotropic  vortical  fluctuations.  This  led  them  to  erroneously 
conclude  that  hneax  analysis  could  not  reproduce  this  rapid  variation  of  kinetic  energy 
seen  in  their  computations.  This  error  was  subsequently  corrected  by  them  (Lee,  Lele 
and  Moin,  1994).  Equations  (7)  and  (9)  suggest  that  behind  the  shock  wave  may  be 
decomposed  into  four  components:  a  spatially  uniform  vortical  component,  an  acoustic 
component  that  decays  exponentially  and  non-monotonic  components  corresponding  to 
the  correlation  between  acoustic  and  vortical  fluctuations  over  both  regimes  of  the  solu¬ 
tion.  This  decomposition,  was  performed  by  Mahesh,  Lee,  Lele  and  Moin  (1994b)  in  their 
analysis  of  acoustic  fluctuations  interacting  with  a  shock  wave.  The  correlation  between 
acoustic  and  vortical  fluctuations  produced  by  waves  in  the  decaying  regime  of  the  solution 
were  seen  to  produce  the  rapid  variation  of  kinetic  energy. 

The  variation  of  kinetic  energy  behind  the  shock  wave  may  be  further  clarified  by  ex¬ 
amination  of  the  equations  governing  its  evolution.  Rearrangement  of  the  linearized  Euler 
equations  (Appendix  A)  shows  that  the  quantity 

nfM 


is  conserved  along  a  mean  streamline,  /total  changes  across  the  shock  wave  and  remains 
constant  downstream.  As  wiU  be  seen  in  the  next  section,  decays  exponentially  behind 
the  shock  wave.  The  rapid  non-monotonic  variation  of  may  therefore  be  explained 
as  a  result  of  partitioning  between  potential  and  kinetic  energy  brought  about  by  the 
correlation  between  acoustic  fluctuations  and  the  shock-normal  component  of  fluctuating 
velocity.  This  is  illustrated  in  Figure  12  where  the  terms  in  equation  (35)  are  plotted. 

The  p'u'  correlation  may  be  decomposed  into  two  components  as  follows: 


q^  p'2 


j2p2\ 


+ 


p'u' 

~P^ 


(35) 
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where  the  subscripts  ‘prop’  and  ‘decay’  refer  to  the  propagating  and  decaying  regimes 
respectively.  The  above  decomposition  is  performed  by  Mahesh,  Lele  and  Moin  (1994c) 
and  reveals  that  the  correlation  between  acoustic  fluctuations  and  the  velocity  field  in  the 
decaying  regime  causes  the  rapid  variation  of  behind  the  shock  wave. 


8.2  Far-field  kinetic  energy 

The  influence  of  entropy  fluctuations  on  far-field  levels  of  kinetic  energy  amplification 
is  examined  as  a  function  of  mean  Mach  munber.  K  A^/A^  is  independent  of  fc,  it  is 
easily  shown  that  tiirbulent  statistics  in  the  far-field  are  independent  of  the  upstream 
energy  spectrum.  Equation  (32)  outlined  the  expression  for  u'.^ .  Substituting  for  En 
from  equations  (7)  and  (9)  into  equation  (32),  the  following  expression  is  obtained: 


^  =  /  (|Ep-MGnsin^idV>i+/ 

/  far-field  L«/o  J ijic 


\G^  sinV>i  d^i 


E{k)dik  (37) 


In  writing  the  above  expression,  use  is  made  of  the  fact  that  the  correlation  term  in  the 
propagating  regime  integrates  to  zero  in  the  far-field.  Since  E{k)Ak  =  q\ /2,  the  above 
expression  is  independent  of  the  upstream  energy  spectrum.  The  above  discussion  holds 
also  for  the  other  components  of  kinetic  energy  and  the  thermodynamic  fluctuations.  Thus 
if  AcjAv  is  independent  of  fc,  the  far-field  levels  depend  only  upon  Mi  and  A^Av. 

The  far-field  levels  of  kinetic  energy  amplification  are  plotted  in  Figures  13  and  14 
respectively.  The  presence  of  entropy  fluctuations  is  seen  to  significantly  influence  the 
levels  of  kinetic  energy  behind  the  shock  wave.  Recall  that  the  far-field  kinetic  energy  may 
be  decomposed  into  vortical  and  acoustic  components.  Such  decomposition  reveals  that 
while  acoustic  waves  axe  generated  at  the  shock  wave,  their  contribution  to  the  far-field 
kinetic  energy  is  not  very  significant;  amplification  and  generation  of  vortical  fluctuations 
accounts  for  most  of  the  kinetic  energy  behind  the  shock  wave. 

Entropy  fluctuations  also  influence  the  anisotropy  of  the  turbulence  behind  the  shock. 
As  shown  in  Figure  15,  negatively  correlated  fluctuations  of  u'  and  T'  upstream  of  the 
shock  wave  increase  the  fraction  of  energy  in  the  streamwise  direction  behind  the  shock. 
Positively  correlated  fluctuations  have  the  opposite  effect. 


8.S  Statistics  of  the  shock  front 

The  oscillation  of  the  shock  wave  is  considerably  affected  by  the  presence  of  entropy 
fluctuations  in  the  incident  field  of  disturbemces.  Results  for  the  velocity  and  amplitude  of 
oscillation  of  the  shock  front  are  presented  below. 

From  equation  (5),  the  rms  level  of  the  oscillation  velocity  is  given  by: 

^  yoo  i‘n/2  ^  y7r/2  ^  yoo 

=  47r  /  /  |Lp  sin^i  dArdV’i  =  /  |jLp  sin  V’l  dV’i  /  E{k)  dk  (38) 

Jk=oJipi=0  Jij)i=0  Jk=0 


Combined  interaction  of  vorticity  and  entropy  fluctuations  with  a  shock  wave 


63 


i.e.,  the  rms  velocity  of  oscillation  of  the  shock  wave  is  independent  of  the  incident  energy 
spectrum  if  does  not  depend  upon  k.  The  influence  on  incident  entropy  fluctuations 

on  the  rms  level  of  is  illustrated  in  Figure  16.  The  oscillation  velocity  is  seen  to  be 
of  the  order  of  the  incident  velocity  fluctuation.  Entropy  fluctuations  with  negatively 
correlated  u'  and  T'  enhance  the  oscillation  velocity  of  the  shock  wave.  Positively  correlated 
fluctuations  are  observed  to  suppress  it. 

The  amplitude  of  oscillation  of  the  shock  wave  is  similarly  affected.  Integration  of 
equation  (5)  yields  the  following  expression  for  the  rms  oscillation  amplitude: 

F  =  47r/°°  KpPsin^idA:d^i=  SsinV-id^/.!  H  ^  dk 

Jk=oJ ,1,1=0  «  7^1=0  Jk=0 

(39) 

Unlike  the  oscillation  velocity,  the  amplitude  of  oscillation  is  dependent  upon  the  incident 
spectrum.  As  shown  in  Figure  17,  the  oscillation  amplitude  is  increased  by  the  incidence  of 
temperature  fluctuations  that  are  negatively  correlated  with  u'.  Decrease  in  the  amplitude 
of  oscillation  is  observed  when  u'  and  T'  in  the  incident  field  are  positively  correlated. 

8.4  Thermodynamic  fluctuations 

The  streamwise  variation  of  pressmre,  density  and  temperatxrre  fluctuations  behind  a 
Mach  1.5  shock  wave  is  shown  in  Figure  18.  The  incident  field  comprises  of  vorticity 
and  entropy  fluctuations  (Ar  =  0.58,  <j)r  =  0).  Pressure  fluctuations  decay  exponentially 
behind  the  shock  wave  while  density  and  temperature  fluctuations  exhibit  non-monotonic 
behavior.  The  pressure  fluctuations  in  the  near-field  are  seen  to  be  intense  while  the 
far-field  intensity  is  smaller  than  those  of  density  and  temperature. 

The  preceding  sections  showed  how  entropy  fluctuations  affect  the  vortical  part  of  the 
flow-field.  Figure  19  shows  that  this  influence  extends  to  the  acoustic  component.  The 
fax-field  intensity  of  pressure  is  plotted  as  a  function  of  Mach  niunber.  The  combined 
interaction  of  vorticity  and  entropy  fluctuations  is  compared  to  that  of  vortical  fluctuations 
alone.  Fluctuations  with  a  negative  correlation  between  u'  and  T'  are  seen  to  increase  the 
level  of  sound  in  the  far-field  while  positively  correlated  fluctuations  suppress  it.  Over  the 
range  of  Mach  numbers  shown,  this  difference  in  sound  level  is  calculated  to  be  between  4 
and  5  decibels. 


8.5  Morkovin’s  hypothesis  across  a  shock  wave 

The  applicability  of  Morkovin’s  hypothesis  (equation  1)  across  a  shock  wave  is  examined 
in  this  section.  The  incident  field  of  vorticity  and  entropy  fluctuations  is  constrained  to 
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satisfy  Morkovin’s  hypothesis  i.e.,  Ar  is  chosen  such  that  the  relation 

^  =  (40) 

Pi  ^1 

is  satisfied  at  every  wavenumber.  The  fluctuations  in  the  fax-field  behind  the  shock  are  then 
examined  to  see  if  the  hypothesis  holds.  The  results  (Figme  20)  show  that  the  first  part 
of  the  hypothesis  i.e.,  p'/p  =  —T'/T  is  still  a  good  approximation  behind  the  shock  wave 
especially  if  the  mean  Mach  number  is  less  than  2.  However  the  part  of  the  hypothesis  that 
relates  T'  to  u'  does  not  seem  to  hold  behind  the  shock  wave.  This  behavior  is  explained 
below. 

The  equation  p'/P  —  —T'/T  is  obtained  by  setting  p'  to  zero  in  the  linearized  equation  of 
state.  It  amounts  to  neglecting  the  acoustic  mode  in  comparison  to  the  entropy  mode.  As 
seen  from  equation  (2d)  it  is  an  identity  as  far  as  the  initial  disturbance  field  is  concerned. 
Upon  interaction  with  the  shock,  the  incident  field  of  vorticity  and  entropy  fluctuations 
generate  acoustic  waves.  The  generation  of  acoustic  waves  is  however  accompanied  by 
amplification  of  the  incident  entropy  fluctuations.  Also  a  fraction  of  the  acoustic  waves 
decays  behind  the  shock.  As  a  result  the  acoustic  contribution  in  the  fax-field  to  the 
thermodynamic  fluctuations  becomes  significant  only  at  the  larger  Mach  numbers.  The 
first  part  of  Morkovin’s  hypothesis  is  therefore  a  good  approximation  behind  shock  waves 
of  moderate  strength. 

The  relation, 

=  =  -(7-1)M^^  (41) 

is  obtained  by  assuming  negligible  fluctuations  in  stagnation  temperature  in  the  linear 
limit.  Denoting  the  stagnation  temperature  by  To, 


2a 


(42) 


Linearizing  the  above  equation  yields  the  following  expression  for  fluctuations  in  stagnation 
temperature: 


ro  =  T'  +  ^  (43) 

Setting  Tq  to  zero  axid  rearranging  yields  equation  (41). 

The  Rankine  Hugoniot  equations  may  be  used  to  show  that  equation  (41)  cannot  be 
valid  behind  a  shock  wave  if  it  is  assumed  to  hold  upstream  of  the  shock.  The  energy 
equation  requires  the  stagnation  temperature  to  be  constant  across  the  shock  in  a  frame 
of  reference  that  moves  at  the  instantaneous  speed  of  the  shock  wave,  i.e., 


Ti  -h  r;  + 


(Ui  +  tt'i  -  -f  +  w'i' 

2Cr, 


T2  +  T,'  -f 


{U2  +  u'2  —  i^f)^  -|-  v'2^  4-  'w'2^ 

2C„ 


(44) 
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Linearizing  the  above  equation  and  constraining  the  incident  fluctations  to  satisfy  equation 
(41)  yields  (Mahesh,  Lele  and  Moin;  1994c): 

n+^= -  U2)  (46) 

Op  Op 

i.e.,  the  fluctuations  of  stagnation  temperature  are  not  zero  behind  the  shock  wave.  Di¬ 
viding  through  by  T2  and  rearranging,  we  get: 

|+(7-1)^||  =  -(7-1)M.(|-i)|  (46) 

The  applicability  of  Morkovin’s  hypothesis  behind  the  shock  wave  requires  that 

(7-1)M2(^-i)^  ~  0  (47) 

As  shown  in  Figure  20,  the  above  term  is  not  negligible,  leading  to  inapplicability  of  the 
hypothesis  behind  the  shock  wave. 

9.  Summary 

Inviscid  linear  analysis  was  used  to  study  the  simultaneous  interaction  of  vorticity  and 
entropy  fluctuations  with  a  shock  wave.  The  two-dimensional  interaction  of  a  shock  wave 
with  a  plane  vorticity-entropy  wave  was  initially  considered.  Expressions  were  derived  for 
the  flow-field  behind  the  shock  and  the  displacement  of  the  shock  front.  The  statistics  of 
kinetic  energy  and  thermodynamic  fluctuations  behind  the  shock  wave  were  examined.  The 
amplitude  ratio  and  phase  difference  between  the  incident  entropy  and  vorticity  fluctua¬ 
tions  were  shown  to  significantly  influence  the  statistics.  Dependence  upon  the  amplitude 
ratio  was  monotonic.  Entropy  fluctuations  in  phase  with  the  incident  vorticity  fluctua¬ 
tions  (u'  and  T'  negatively  correlated)  were  observed  to  enhance  levels  of  kinetic  energy 
and  thermodynamic  fluctuations  behind  the  shock.  Positively  correlated  fluctuations  were 
observed  to  have  a  suppressing  effect. 

An  explanation  was  provided  to  explain  these  trends.  The  evolution  of  fluctuating 
vorticity  across  the  shock  wave  was  noted  to  have  two  important  contributions:  bulk 
compression  of  incident  vorticity  and  baroclinic  production  of  vorticity  through  the  incident 
entropy  fluctuations.  The  upstream  correlation  between  vorticity  and  entropy  fluctuations 
was  shown  to  determine  whether  these  two  sources  of  vorticity  enhance  or  oppose  each 
other  thereby  determining  kinetic  energy  levels  behind  the  shock  wave.  A  scaling  was  then 
proposed  for  the  evolution  of  vorticity  across  the  shock  wave. 
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The  analysis  was  extended  to  consider  the  interaction  of  a  shock  wave  with  an  isotropic 
field  of  vorticity  and  entropy  fiuctuations.  Both  spatial  variation  behind  the  shock  and 
fax-field  levels  of  turbulent  statistics  were  examined.  The  presence  of  entropy  fiuctuations 
was  shown  to  significantly  affect  the  levels  of  kinetic  energy,  vorticity  and  thermodynamic 
fiuctuations  behind  the  shock.  Higher  levels  of  kinetic  energy  and  vorticity  were  observed 
when  u'  and  T'  in  the  incident  fiuctuations  were  negatively  correlated.  Positive  correlation 
had  the  opposite  effect.  Acoustic  waves  were  produced  through  the  interaction.  The  level  of 
sotmd  behind  the  shock  wave  was  enhanced  by  the  incidence  of  fluctuations  with  negative 
correlation  between  u'  and  T' .  Suppression  of  sound  level  was  observed  if  the  incident 
fluctuations  were  positively  correlated.  The  incidence  of  entropy  fluctuations  was  observed 
to  affect  the  oscillation  of  the  shock  front.  Negatively  correlated  fluctuations  of  u'  and  T' 
increased  the  velocity  and  amplitude  of  oscillation  while  positively  correlated  fluctuations 
suppressed  it. 

The  applicability  of  Morkovin’s  hypothesis  behind  a  shock  wave  was  examined.  Statistics 
in  the  fax-field  indicate  that  neglect  of  the  acoustic  mode  is  a  good  approximation  in  the 
fax-field  of  shock  waves  of  moderate  strength  (Mi  <  2).  The  part  of  the  hypothesis  relating 
u'  and  T'  was  seen  to  be  invalid  behind  the  shock.  Non-negligible  oscillation  of  the  shock 
front  was  shown  to  be  responsible. 

The  results  indicate  the  considerable  importance  of  entropy  fluctuations  in  the  evolution 
of  the  turbulent  fluctuations  in  a  compressible  boundary  layer  upon  its  interaction  with  a 
shock  wave. 
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The  procedure  for  analysis  of  the  two-dimensional  interaction  of  a  shock  wave  with  a 
plane  vorticity-entropy  wave  is  sinnmarized  below.  A  more  detailed  description  is  furnished 
by  Mahesh,  Lele  and  Moin  (1994c). 

A.l.  The  governing  equations 

The  governing  equations  behind  the  shock  wave  are  the  Euler  equations  that  are  lin¬ 
earized  about  the  uniform  flow  behind  the  shock,  i.e., 


Pt  +  U  -p  {U'^  +  V'y) 

(A.la) 

u't  +  U  u'j.  =  -r  p'x 
p 

(A.lb) 

v',  +  Uv',  =  -^p'y 

(A.lc) 

s't  +  U  s'^=0 

(Aid) 

-  and 

~  ■p'^  T  ^  Cp~  jP  p 

(A.le) 

The  variables  s'  and  Cp  in  the  above  equations  represent  entropy  fluctuations  and  the 
speciflc  heat  at  constant  pressure  respectively. 

Linearizing  the  Rankine  Hugoniot  equations  yields  the  following  relations  that  the  dis¬ 
turbance  field  must  satisy  at  the  mean  position  of  the  shock  front  (x  =  0). 

^  (t  -  i)Mi  -  2  2 

Vi  (7  +  l)A^f  Vi  /■^(7  +  1)M?  \Ti) 

v'i  v'i  2(Mf  - 1)  , 

t^i  (7+1)M?'“ 

6.  =  ^  _  (7  -  + 4  fn\ 

Pi  {j-l)Ml  +  2\  Vi  )  (7  -  l)Mf  +  2  J 

A  ^  _  27M?  fn\ 

Pi  27Aff-(7-l)  V  Vi  )  2yMl-(j-l)  \TiJ 

A. 2  The  two  regimes  behind  the  shock  wave 
The  solution  behind  the  shock  wave  was  noted  to  have  two  different  regimes  that  differed 
in  the  nature  of  the  pressure  field.  The  existence  of  the  two  regimes  is  illustrated  as  follows. 
The  pressure  field  behind  the  shock  wave  satisfies  the  wave  equation;  i.e., 

p'tt  -1-  2U2  p'^t  -  (ai  -  U^)  P'xx  -  4  P'yy  =  0 


{A.2a) 

{A.2b) 

(A.2c) 

{A.2d) 


(4.3) 
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Continuity  at  the  shock  wave  requires  that  the  solution  behind  the  shock  wave  have  the 
same  transverse  wave  number  and  frequency  as  the  incident  disturbance.  The  pressure 
may  therefore  be  expressed  as  , 


p' =  F(a;)  (A4) 

Substituting  the  above  expression  for  pressure  into  equation  (A.3)  yields  the  following 
equation  for  F. 

-  g]  F"  +  i  F'  +  =  0  (A5) 

The  solution  to  the  above  equation  is  of  the  form  F  ~  Substitution  into  equation 

(A. 5)  yields  the  following  equation  for  k. 

+  (A6) 

The  discriminant  of  the  above  quadratic  equation  determines  if  k  is  real  or  imaginary.  A 
real  valued  k  would  imply  that  the  pressure  field  is  a  plane  wave.  The  discriminant  is  given 
by  the  following  expression: 

Discriminant  =  2A:/— 

Recall  that  m  =  cos  and  I  =  sin  tpi .  The  requirement  for  real  k  is  therefore, 

cot^  V-.  >  (^  -  3)  (4,8) 

i.e.,  0  <  where  ipc  is  a  root  of 


Note  that  is  a  function  of  the  Mach  number  of  the  shock  wave.  If  the  incident  angle 
exceeds  V’c,  k  becomes  imaginary.  Denoting  the  real  and  imaginary  parts  of  k  by  kr  and 
and  ki  respectively,  the  function  F  is  now  given  by. 


i.  e.  the  pressure  field  decays  exponentially  behind  the  shock  wave. 
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A.S  The  solution  behind  the  shock  wave 

The  arguments  of  the  vorticity,  entropy  and  acoustic  modes  over  both  regimes  of  the 
solution  are  derived  below.  The  flow-field  may  then  be  expressed  as  a  superposition  of 
vortical,  acoustic  and  entropic  components. 

It  was  noted  that  the  vorticity  and  entropy  fluctuations  were  plane  waves  over  both 
regimes.  Since  vorticity  and  entropy  fluctuations  axe  both  convected  by  the  mean  flow, 
they  axe  of  the  same  form  which  is  given  by: 


^ikx  ^ik(ly—mUit) 


(All) 


Note  that  the  transverse  wavenumber  and  frequency  have  been  matched  to  that  of  the 
incident  disturbance.  The  variable  k  is  obtained  as  follows.  Convection  by  the  mean  flow 
requires  that  the  vorticity  and  entropy  fluctuations  have  the  following  dispersion  relation: 


u)  =  U\km  =  U2k 


i.e.. 


k  =  kmy^  =  kmr 
U2 


The  vorticity  and  entropy  fields  axe  therefore  of  the  form: 


^ik{mrx+ly—mUi  t) 


(A.12) 


(A13) 


The  form  of  the  pressure  field  was  derived  in  the  previous  section.  However,  explicit 
expressions  were  not  derived  for  the  variables  k,  k,  and  kr.  Evaluation  of  these  variables 
would  completely  determine  the  functional  form  of  the  pressure  field.  Recall  that  k  was 
defined  as  the  root  of  equation  (A.6);  i.e., 


(AM) 


The  other  root  is  rejected  on  the  grounds  that  it  leads  to  a  physically  unacceptable 
solution  that  grows  exponentially  behind  the  shock  wave.  The  above  expression  is  real 
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only  if  0  <  ■01  <  '0c-  If  the  angle  of  incidence  exceeds  0c,  ^  is  complex.  Equation  (^.14) 
shows  that  the  real  and  imaginary  parts  of  k  are  given  by: 


kr  _  U2  1  _Ui  Ml 

k  yl\- 


(^.15) 


h 

k 


,v,  M2  lulf  1  A  ^ 

£/2l-Af|y  iJf  7  P 


(.4.16) 


The  argument  of  the  vorticity,  entropy  and  acoustic  modes  behind  the  shock  wave  are 
now  completely  determined  over  both  regimes  of  the  solution.  Expressions  for  the  velocity 
and  thermodynamic  field  as  given  by  equation  (4)  follow. 


A. 4  Expressions  for  the  coefficients 

The  coefficients  F,G,H,I,K^Q  and  L  are  obtained  by  substituting  the  solution  into 
the  governing  equations  behind  the  shock  wave  and  the  boundary  conditions  at  the  mean 
position  of  the  shock  front  (r  =  0).  Since  the  argument  of  the  solution  has  been  matched 
to  that  of  the  incident  disturbance  this  yields  a  set  of  equations  relating  the  coefficients. 
This  system  of  equations  is  then  solved  thereby  completing  the  solution. 

The  governing  equations  and  boundary  conditions  yield  the  following  set  of  equations 
(Mahesh,  Lele  and  Moin;  1994c): 


F  =  aK 

(A.Ua) 

H  =  0K 

{A.17b) 

k 

II 

(A.nc) 

F  +  G-L  =  Bi{1-L)-B2^ 

/l.y 

(A.17d) 

-+Q  =  Ci{l-l)-C2^ 

7  Ay 

(A.17e) 

K  =  D,{1-L)-D2^ 

Ay 

{A.17f) 

H  +  I  =  —m  —  El  —L 
m 

iA.17g) 
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The  variables  a  and  axe  defined  as: 


a? 


a  = 


0  = 


a% 


The  variables  Bi,B2,Ci,C2,Di,D2  and  Ei  in  the  above  equations  correspond  to  the 
following  compact  representation  of  the  boundary  conditions  and  may  be  obtained  by 

comparison  to  Equations  (A.2): 

Ui  Vi  ^  Ti 

(A.lSa) 

Pi  _  it  ,  f,  r[ 

+c._ 

(A.186) 

Pi  Ti 

(A.lSc) 

-  "i  4.  p  £ 

(A.18d) 

Equations  (A.17)  are  solved  for  the  coefficients.  Details  of  the  solution  are  provided  by 
Mahesh,  Lele  and  Moin  (1994c).  The  variable  L  is  first  solved  for  and  has  the  following 
form: 


L  = 


-m  -  ^(Dil-  D2^)  +  Hfl-aW-  +  Bil  - 


Eli  -  m  - ’fil  -  B,  +  aDi) 


(A19) 


Once  L  is  known,  the  other  coefficients  are  easily  obtained  using  equations  (A.17).  This 
completes  the  solution. 
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Shock  wave 


FIG.  1.  Schematic  of  the  interaction  of  a  vorticity-entropy  wave  with  a  shock  wave. 
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FIG.  2.  Evolution  of  behind  the  shock  wave  in  the  propagating  regime.  The  incident  angle  is 

45°  and  the  mean  Mach  number  is  1.5.  - (Pure  vorticity), - {Ar  =  1;  =  0)? 

.  (Ar  =  1;  <^r  =  180°), - (Pure  entropy). 


kx/2Tr 


FIG.  3.  Evolution  of  5^  behind  the  shock  wave  in  the  decaying  regime.  The  incident  angle  is 

75°  and  the  mean  Mach  number  is  1.5.  - (Pure  vorticity), - (Ar  =  1;  =  0), 

.  (Ar  =  l;4>r  =  180°), - (Pure  entropy). 
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FIG.  4.  The  amplification  of  kinetic  energy  across  the  shock  wave  as  a  function  of  the  angle  of 
incidence.  Far-field  kinetic  energy  levels  axe  plotted  and  the  mean  Mach  number  is  1.5. 

- (Pure  vorticity), - (Ar  =  1;  =  0), .  {Ar  =  1;  =  180°), - (Pure 

entropy). 


V’l 

FIG.  5.  The  amplification  of  vorticity  across  the  shock  wave  as  a  function  of  the  angle  of  inci¬ 
dence.  The  mean  Mach  number  is  1.5.  - (Pure  vorticity), - {Ar  =  1;  =  0), 

.  {Ar  =  l;4>r  =  180°), - (Piure  entropy). 
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FIG.  6.  The  amplification  of  vorticity  across  the  shock  wave  as  a  function  of  the  angle  of  inci¬ 
dence.  The  mean  Mach  number  is  1.5.  The  line  corresponds  to  linear  analysis  while  the 
symbols  axe  from  numerical  computation. 
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Bulk  compression 

a 

+ 


Shock  wave 


(a) 


+ 


Shock  wave 


(b) 

FIG.  7.  A  spherical  element  of  fluid  passing  through  a  shock  wave.  The  effects  of  bulk  compres¬ 
sion  and  baroclinic  vorticity  production  are  shown,  u'  and  T'  are  negatively  correlated 
upstream  of  the  shock  wave  in  (a)  while  (6)  corresponds  to  positive  correlation  between 
u'  and  T' . 
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FIG.  8a.  Unsealed  amplification  of  the  rms  level  of  vorticity  in  the  interaction  of  a  vorticity 

wave  with  a  shock  wave.  -  (Mi  =  1.25), - (Mi  =  1.5),  .  (Mi  =  2), 

- (Ml  =  2.5). 


FIG.  8b.  Scaled  amplification  of  the  rms  level  of  vorticity  in  the  interaction  of  a  vorticity  wave  with 
a  shock  wave. - (Mi  =  1.25), - (Mi  =  1.5), .  (Mi  =  2), - (Mi  =  2.5). 
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FIG.  9a.  Unsealed  rms  level  of  vorticity  produced  in  the  interaction  of  an  entropy  wave  with  a 
shock  wave. - (Mi  =  1.25), - (Mi  =  1.5), .  (Mi  =2), - (Mi  =  2.5). 


FIG.  9b.  Scaled  rms  level  of  vorticity  produced  in  the  interaction  of  an  entropy  wave  with  a  shock 

(Ml  =  1.25), - (Ml  =  1.5), .  (Ml  =  2), - (Mi  =  2.5). 


wave. 


80 


K.  Mahesh,  S.K.  Lele  and  P.  Main 


kQx/2Tz 


FIG.  11.  Evolution  of  behind  a  Mach  1.5  shock  wave.  -  (Pure  vorticity), - {^Ar  = 

0.58;  (f>r  =  0), .  {Ar  =  0.58;  4>r  =  180°), - (Pure  entropy). 


kQx/2'K 

FIG.  12.  The  streamwise  variation  of  terms  in  Equation  3.57  behind  a  Mach  1.5  shock  wave. 

The  incident  fluctuations  correspond  to  Ar  =  0.58  and  <f)r  =  0.  All  terms  are  non- 
dimensionalized  as  indicated  in  Equation  3.57  and  normalized  with  y/^/Ui.  - (5^), 

- (p'^), .  {p'W), - (/total). 
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FIG.  13.  The  amplification  of  across  the  shock  wave  in  the  far-field.  - (Pure  vorticity), 

- {Ar  —  0.58;  <f>r  =  0), .  (Ar  =  0.58;  (f>r  =  180°), - (Pure  entropy). 


Ml 

FIG.  14.  The  amplification  of  v'^  across  the  shock  wave  in  the  far-field.  - (Pure  vorticity), 

- (Ar  =  0.58;  (l)r  =  0), .  (^r  =  0.58;  <j>r  =  180°), - (Pure  entropy). 


K.  Mahesh.  S.K.  Lele  and  P.  Main 


FIG.  15.  Far-field  values  of  u''^ I q'^  behind  the  shock  wave.  - (Pure  vorticity), 

0.58;  4>r  =  0), .  (Ar  =  0.58;  =  180°), - (Pure  entropy). 


1.5  Pn 


S  1.0 


FIG.  16.  The  rms  level  of  fluctuations  in  shock  speed.  -  (Pure  vorticity), 

0.58;  (t)r  =  0), .  {Ar  =  0.58;  <f)r  =  180°), - (Pure  entropy). 
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FIG.  17.  The  rms  level  of  displacement  of  the  shock  front.  - (Pure  vorticity), - {Ar  = 

0.58;  <f>r  =  0), .  (Ar  =  0.58;  <f)r  =  180°), - (Pxire  entropy). 


A:oa:/27r 

FIG.  18.  The  variation  of  thermodynamic  fluctuations  behind  a  Mach  1.5  shock  wave.The  incident 
fluctuations  correspond  to  Ar  =  0.58  and  <j)r  =  0.  All  three  curves  are  normalized  with 

-  {pyP2),- - (p'2/P2), .  (n/T2). 


1  2  3  4  5 

Ml 


FIG.  19.  The  intensity  of  far- field  pressure  fluctuations  behind  the  shock  wave.  - (Pure  vor- 

ticity),  (^Ar  =  0.58;  <})r  =  0), .  (Ar  =  0.58;  <j)r  =  180°), - (Pure  entropy). 


1  2  3  4  5 
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FIG.  20.  Evaluation  of  Morkovin’s  hypothesis  in  the  far-field  of  the  shock  wave.  - 

- {TUT2I .  (7  -  l)MluyU2, - (7  -  1)M2(C/i/172  -  1)6/02. 


{p'l!  P2)i 
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PART  FOUR 


THE  INTERACTION  OF  AN  ISOTROPIC  FIELD 
OF  ACOUSTIC  WAVES  WITH  A  SHOCK  WAVE 


B.l  Introduction 

This  section  uses  linear  analysis  to  study  the  interaction  of  a  shock  wave  with  an 
isotropic  field  of  acoustic  waves.  The  analysis  aims  at  characterizing  the  effect  of  com¬ 
pressibility  in  the  interaction  of  a  compressible  turbulent  flow  with  a  shock  wave.  Within 
the  linear  framework,  a  field  of  acoustic  waves  would  correspond  to  the  dilatational  com¬ 
ponent  of  a  compressible  turbulent  flow.  As  discussed  in  Chapter  3,  there  is  evidence 
to  suggest  that  compressible  effects  are  negligible  in  non-hypersonic  turbulent  boimdary 
layers.  However  shear  flows  such  as  the  mixing  layer  and  supersonic  jet  are  known  to  be 
significantly  affected  by  compressibility. 

The  increased  level  of  compressibility  of  these  flows  (parametrized  by  the  ratio  of 
kinetic  energy  in  solenoidal  and  dilatational  components)  suggests  that  effects  due  to  the 
dilatational  component  would  be  important  in  their  interaction  with  a  shock  wave.  Past 
analytical  work  related  to  shock/turbulence  interaction  (Chapter  1)  has  considered  the 
interaction  of  a  shock  wave  with  a  plane  acoustic  wave.  To  the  best  of  our  knowledge,  the 
interaction  of  a  field  of  acoustic  waves  with  a  shock  wave  has  not  been  analyzed  before. 
Our  analysis  is  based  upon  Moore’s  (1954)  study  of  the  interaction  of  a  shock  with  a  plane 
acoustic  wave.  Moore’s  analysis  is  extended  to  study  the  evolution  of  an  isotropic  field  of 
acoustic  waves  across  a  shock.  The  results  are  then  apphed  to  predict  the  response  of  a 
compressible  isotropic  turbulent  flow  to  a  shock. 

This  section  is  organized  as  follows.  Section  B.2  formulates  the  problem  of  a  plane 
acoustic  wave  interacting  with  a  shock.  The  results  are  evaluated  by  comparison  to  DNS. 
The  analysis  is  then  extended  to  examine  the  evolution  of  an  isotropic  field  of  acoustic  waves 
across  a  shock.  The  statistics  of  turbulent  kinetic  energy  and  thermodynamic  fluctuations 
behind  the  shock  wave  are  described  in  Section  B.3.  A  brief  discussion  including  description 
of  shock  wave/compressible  turbulence  interaction  is  presented  in  Section  B.4.  Summary 
of  the  study  in  Section  B.5  is  followed  by  the  concluding  section  (B.6)  where  some  of  the 
details  of  Moore’s  analysis  are  outlined. 

Parallel  may  be  drawn  between  Part  3  and  the  following  analysis.  The  notation  used 
below  might  differ  from  that  in  Part  3;  we  have  retained  Moore’s  notation  and  method  of 
analysis  for  consistency.  The  work  described  below  is  discussed  by  Mahesh,  Lee,  Lele  and 
Moin  (1994c). 
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B.2  Formulation  of  the  problem 


B.2.1  The  interaction  of  an  acoustic  wave  with  a  shock  wave 


Figure  B.l  shows  a  schematic  of  the  problem.  A  normal  shock  propagates  at  speed 
V  into  fluid  that  is  at  rest  in  the  mean.  The  mean  velocity  behind  the  shock  is  denoted 
by  U •  The  subscripts  1  and  2  denote  the  fluid  in  front  and  behind  the  shock  respectively. 
The  fluid  at  rest  is  assumed  to  be  perturbed  by  the  weak  field  of  a  plane  acoustic  wave 
that  is  incident  at  angle  tpi  to  the  shock.  The  flow  field  associated  with  the  acoustic  wave 
is  given  by: 

V  A,  ) 

vi  f  mxi  -ly  +  ait\ 

V  =  - ) 


(B.la) 

(B.lb) 


(B.lc) 


Pi  A  r  f ‘"ixi  —  ly ait\ 

17=^^ - A^j- 

Here,  Pj ,  Pi  and  ai  are  the  mean  pressure,  density  and  sound  speed  in  the  fluid  ahead  of 
the  shock  while  ui  and  v\  are  the  disturbance  velocities  in  the  x\  and  y  directions  respec¬ 
tively.  The  pressure  and  density  fluctuations  associated  with  the  incident  acoustic  wave 
are  denoted  by  pi  and  px .  The  variables  I  and  m  are  related  to  the  direction  of  propagation 
by  I  =  sin  V’l  and  m  —  cos  V-i  and  Ai  represents  the  lengthscale  of  the  disturbance.  Note 
that  the  coordinate  system  is  fixed  in  the  fluid  that  is  at  rest,  yielding  xi  =  Vt  at  the  mean 
position  of  the  shock  wave.  The  amphtudes  of  velocity,  pressure  and  density  are  related 
through  the  governing  equations  for  an  acoustic  wave;  *.e.. 


I  Ax  =  -mA2;  Az  =  7A4;  Ax  =  — 

7M 


(P.2) 


is  the  mean  Mach  number  of  the  shock  wave  and  7  is  the  ratio  of  specific 
heats.  Incidence  of  the  acoustic  wave  causes  the  shock  wave  to  deform  with  a  profile 
that  matches  the  profile  of  the  incident  wave.  The  hnearized  Rankine-Hugoniot  equations 
are  used  to  describe  the  jump  of  the  disturbance  across  the  disturbed  shock.  This  yields 
boundary  conditions  for  the  flow  behind  the  shock  which  is  described  by  the  Euler  equations 
linearized  about  the  uniform  mean  flow.  The  hnearized  Euler  equations  axe  solved  in 
consonance  with  the  boundary  conditions  to  calculate  the  displacement  of  the  shock  wave 
and  the  flow  field  downstream  of  the  shock. 
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The  flow  behind  the  shock  wave  has  two  distinct  regimes  depending  upon  the  angle  of 
incidence  of  the  acoustic  wave.  The  two  regimes  differ  in  the  natme  of  the  pressure  fleld 
(and  hence,  its  associated  velocity,  temperature  and  density  fleld)  behind  the  shock  wave. 
Over  a  range  of  incident  angles  0  <  <  V’ci  or  ^cu  <  fpi  the  pressure  field  behind  the 

shock  is  a  freely  propagating  plane  acoustic  wave.  However,  if  tpcl  <  fpi  <  '4’cu,  the  pressure 
field  behind  the  shock  corresponds  to  an  evanescent  wave  and  decays  exponentially.  Over 
both  these  regimes,  the  vorticity  and  entropy  waves  propagate  without  decaying.  ■0c/  and 
0ctt  are  roots  of  the  following  equation  (Moore,  1954): 

(^)  ~  ~  ^  +  y  CSC0c^  .  {B.Z) 

Note  that  0c/  and  0cu  depend  only  on  the  mean  Mach  number.  The  displacement  of  the 
shock  front  has  a  certain  speed  of  propagation  along  the  shock  front  due  to  the  tmsteady 
nature  of  the  incident  field.  Analogous  to  the  classical  flexural  wall  problem  in  acoustics, 
(Pierce,  1981)  the  nature  of  the  downstream  pressure  field  is  determined  by  how  the  speed 
of  the  disturbance  compares  with  the  mean  speed  of  sound  and  velocity  downstream  of  the 
shock.  Moore  (1954)  illustrates  the  two  regimes  geometrically.  Ribner  (1953)  provides  an 
equivalent  explanation  by  noting  that  the  unsteady  interaction  of  an  oblique  wave  with  a 
normal  shock  can  be  transformed  into  the  steady  interaction  of  an  oblique  wave  with  an 
obUque  shock.  In  the  transformed  coordinates,  depending  upon  the  incident  angle  of  the 
wave,  the  flow  behind  the  oblique  shock  is  either  supersonic  or  subsonic  (yielding  either  a 
wave  equation  or  Poisson  equation  for  the  pressure),  which  corresponds  to  the  two  regimes 
mentioned  above.  Alternatively,  arguments  similar  to  those  in  Chapter  3  could  be  used  to 
demonstrate  the  existence  of  two  regimes. 

The  method  for  solution  of  the  downstream  flow  field  is  given  by  Moore.  In  the  inter¬ 
ests  of  clarity,  we  reproduce  his  results  at  the  end  of  the  Appendix  in  Section  B.6.  Apart 
from  the  existence  of  two  different  regimes,  note  that  the  velocity  field  is  a  linear  superpo¬ 
sition  of  acoustic  and  vortical  components.  Similarly,  the  density  and  temperature  fields 
are  superpositions  of  acoustic  and  entropic  components.  The  pressure  field  is  associated 
solely  with  the  acoustic  component. 
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B.2.2  Comparison  of  linear  analysis  to  numerical  solution 

We  compare  results  of  the  linear  analysis  to  numerical  computations  of  the  interaction 
of  a  normal  shock  wave  with  a  plane  acoustic  wave  that  is  incident  from  upstream.  The 
computation  solves  the  two-dimensional  compressible  Navier-Stokes  equations  in  a  frame 
of  reference  that  moves  at  the  mean  speed  of  the  shock  wave.  The  sixth  order  Fade  scheme 
(Lele,  1992)  is  used  to  compute  spatial  derivatives  and  the  third  order  Runge  Kutta  (Wray, 
1986)  is  used  to  integrate  in  time.  No  shock-capturing  or  shock-fitting  is  used  to  treat  the 
shock  wave;  its  structure  is  resolved  by  a  non-uniform  mesh. 

This  comparison  is  intended  to  complement  that  by  Zang  et  al  (1984)  who  compared 
the  predictions  of  McKenzie  and  Westphal’s  (1968)  Hnear  analysis  to  their  numerical  so¬ 
lution  of  the  two-dimensional  Euler  equations  using  a  shock-fitting  scheme.  Zang  et  al 
examined  the  effect  of  incident  angle,  shock-strength  and  the  amphtude  of  the  incident 
disturbance  in  the  interaction  of  acoustic  and  vorticity  waves  with  a  shock  wave.  Results 
were  presented  only  in  the  freely  propagating  regime.  The  dependence  on  incident  angle 
was  examined  in  the  interaction  of  disturbances  of  amplitude  0.1%  and  10%  with  a  shock 
wave  of  Mach  number  8.  Good  agreement  was  seen  away  from  the  critical  angle;  divergence 
from  the  linear  analysis  prediction  was  seen  within  about  20°  of  the  critical  angle.  The 
linear  analysis  predictions  were  quite  robust  in  terms  of  dependence  on  shock  strength  and 
disturbance  amplitude.  For  an  incident  angle  of  30°,  the  linear  predictions  were  valid  for 
disturbance  amplitudes  as  high  as  25%  for  acoustic  waves  and  100%  for  vorticity  waves 
and  shock  waves  whose  Mach  number  was  close  to  imity. 

We  present  results  of  the  interaction  of  a  plane  acoustic  wave  {pi/Pi  =  2.5%)  inter¬ 
acting  with  a  shock  wave  of  Mach  number  1.5.  The  Mach  munber  of  1.5  was  chosen  to  be 
representative  of  experiments  on  the  shock/turbulence  interaction.  In  the  context  of  Zang 
et  aPs  results,  we  only  present  the  dependence  on  the  incident  angle.  Results  are  presented 
in  both  freely  propagating  and  decaying  regimes.  Figure  B.2  compares  the  computed  val¬ 
ues  of  rms  dilatation  and  vorticity  with  predictions  of  the  linear  analysis.  Note  that  V'c/ 
and  V’cu  have  values  of  97.16°  and  154.45°  respectively  for  M  =  1.5.  Excellent  agreement 
is  seen  except  for  a  region  within  25  degrees  of  the  critical  angles.  As  noted  by  Zang  et  al 
(for  the  shock  wave  of  Mach  number  8),  a  decaying  pressure  field  is  observed  downstream 
of  the  shock  wave  before  the  critical  angle  is  reached. 

A  possible  reason  for  deviation  of  the  computed  result  from  the  linear  analysis  around 
the  critical  angle  is  as  follows.  Within  Hnear  analysis,  the  unsteady  interaction  of  a  plane 
disturbance  with  a  normal  shock  wave  may  be  transformed  into  the  steady  interaction  of 
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the  distiirbance  with  an  oblique  shock  wave.  This  transformation  is  an  integral  part  of 
Ribner’s  (1953)  and  Chang’s  (1957)  analysis  of  the  interaction  of  a  shock  with  a  vorticity 
and  entropy  wave  respectively.  While  Moore  does  not  employ  this  transformation,  it  exists 
and  involves  defining  coordinates  x'  and  y'  such  that 


x'  =  X  —  Vt] 


,  ai+mV^ 

y  =y - ; — i 


i.e.,  the  flow  is  observed  from  a  frame  of  reference  that  moves  at  the  mean  shock  speed  in  the 

X  direction  and  at  speed  (ui  +mV)/l  in  the  y  direction.  In  these  transformed  coordinates, 

the  governing  equations  behind  the  shock  wave  are  the  steady  Euler  equations  that  axe 

linearized  about  uniform  mean  flow  at  an  equivalent  Mach  number  Me  that  is  given  by 

the  relation:  „  „ 

,  P(V-U)^  +  (a,+mVf 

‘  ~  4P 

This  equivalent  Mach  number  depends  upon  the  Mach  number  of  the  normal  shock  wave 
and  the  incident  angle  of  the  disturbance  and  equals  unity  at  the  critical  angle.  As  a  result, 
incident  angles  close  to  the  critical  angle  correspond  to  steady  transonic  flow  downstream 
of  the  shock  wave  in  the  transformed  coordinates.  It  is  well  known  that  the  linear  ap¬ 
proximation  is  inconsistent  in  the  transonic  regime.  We  believe  that  use  of  the  transonic 
small  disturbance  equations  are  necessary  to  accurately  represent  the  interaction  around 
the  critical  angle. 


Fmther  indication  of  this  is  provided  in  Figure  B.3,  where  the  dependence  of  the 
interaction  upon  the  amplitude  of  the  incident  disturbance  is  examined.  The  interaction 
away  from  the  critical  angle  {ipi  =  60®)  is  compared  to  that  close  to  the  critical  angle 
(^1  =  100°).  The  rms  intensity  (prmsfP)  of  the  incident  pressure  fluctuation  is  varied 
from  0.14%  to  14%  in  both  cases.  The  vorticity  (normalized  by  the  incident  dilatation) 
behind  the  shock  wave  is  plotted  against  the  amplitude  of  the  incident  disturbance  in 
Figure  B.3.  The  variation  with  the  incident  amplitude  is  negligible  away  from  the  critical 
angle.  However,  significant  dependence  upon  the  incident  amplitude  is  seen  close  to  the 
critical  angle.  This  non-linear  behavior  around  the  critical  angle  is  in  accordance  with 
the  expectation  that  the  transonic  small  disturbance  equations  are  needed  to  describe  the 
interaction  in  that  range. 

B.2.3  The  interaction  of  a  field  of  acoustic  waves  with  a  shock  wave 


Moore’s  analysis  is  now  extended  to  describe  the  interaction  of  an  isotropic  field  of 
acoustic  waves  with  the  shock  wave.  The  problem  is  convenient  to  analyze  in  cylindrical 
coordinates.  Figure  B.4  shows  a  schematic  of  the  problem  in  a  frame  of  reference  that 
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moves  at  the  speed  of  the  undisturbed  shock.  The  y  -  z  plane  corresponds  to  the  shock 
wave,  while  the  acoustic  wave  lies  in  the  xi  -  r  plane  that  makes  an  angle  <f)  with  the  y 
axis.  Irrotationality  precludes  any  component  of  velocity  outside  this  plane.  In  the  xi-r 
plane,  the  acoustic  wave  makes  an  angle  xj}\  with  the  xi  axis.  The  xi  —  r  plane  is  equivalent 
to  the  xi  -y  plane  used  in  Moore’s  analysis.  It  is  seen  from  Figure  B.4  that 


Ur 


V 

cos  (f) 


w 

sin 


u^  =  0. 


(BA) 


Moore’s  results  are  used  to  obtain  the  following  relations  describing  the  amplitude  of  the 
fluctuating  field  behind  the  shock  for  a  plane  incident  wave.  Note  that  we  transform 
to  a  coordinate  system  that  moves  at  the  mean  speed  of  the  shock  t.  e.,  our  streamwise 
coordinate  becomes  7/2  =  —[x2  —  {V  —  U)t]. 

For  0  <  <  V’c/  or  V’cu  <  V’l  <  tt, 
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The  constants  K^F  and  H  are  associated  with  the  acoustic  component  of  the  flow 
field.  Similarly,  G  and  I  are  associated  with  the  vortical  component  and  Q  corresponds  to 
the  entropic  component. 

If  tpcl  <1pl  <  ^  cu,  the  amplitudes  behind  the  shock  are  given  by: 


M  =  ^-2dr„/X^ 


^(l)+^(2) 


|Pl 


Pi- 


where. 


F2I 

V2 


T,  =  +  Gl,y, 


T1+T2+  Tz 


|Pi 


Pi- 


T2  = 


^(1)  +  ^(2) 
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•^(1)<^(1)+-^^(2)G(2)^  COsS-  ^F(2)G(i)-jP(1)(j(2)^  sill^ 


6  = 


c?72  1 

Ai  1  —  r 


The  amplitude  oiur^lV  may  be  obtained  from  the  expression  for  U2/V  by  replacing 
^(i)j  G'(i)5-^(2)  and  <^(2)  by  and  7(2)  respectively.  The  amplitude  of  P2IR2 

may  be  obtained  from  the  expression  for  u^jV  by  replacing  •F(i),^t(i),7(2)  and  G(2)  by 
K{,x)h,Q(,i),K{2)h  and  §(2)  respectively. 

Analogous  to  the  propagating  regime,  the  constants  K(^a)-,R'{ot)  and  are  associated 
with  the  acoustic  component  of  the  flow  field.  G(^a)  and  /(a)  are  associated  with  the  vortical 
component  and  Q{a)  corresponds  to  the  entropic  component. 

Note  that  the  solution  downstream  of  the  shock  wave  requires  the  spectrum  of  pressure 
fluctuations  ahead  of  the  shock.  We  assume  the  incident  acoustic  field  to  be  isotropic. 
This  combined  with  the  condition  of  irrotationality  requires  the  upstream  spectrum  of  the 
velocity  fluctuations  to  be  of  the  form: 


E{k)  kikj 


where  Eij{]c)  is  the  energy  spectriim  tensor  and  E{k)  is  the  three-dimensional  energy 
spectrum.  Using  the  acoustic  relation,  A3  =  -yMAi/  cos^i,  and  the  above  expression  for 
the  energy  spectrum  tensor,  we  get. 


IPil^(k) 

Pi^ 


m 

4:7rk^ 


{B.7) 


Substituting  for  |pi  |^(k)/Pi^  and  integrating  over  wavenumber  space  at  every  stream- 
wise  location  (772)5  we  get  the  streamwise  evolution  of  statistics  downstream  of  the  shock 
wave.  This  integration  is  done  in  polar  coordinates  in  the  wave  number  space: 


ki  =  k  cos  V’l  5 


k2  =  k  sin  V’l  cos  <f>. 


kz  =  fcsinV’i  sin^. 


d^k  =  sin  ipi  d^’i  d^  dA; 


where,  k  =  |k|  varies  from  0  to  00,  ipi  from  0  to  tt  and  (f>  varies  from  0  to  27r.  Since  the 
coefficients  in  the  above  equations  are  independent  of  <j>,  the  integration  over  <j)  can  be  done 
analytically  for  isotropic  initial  spectra.  For  example, 

^^2(^72)  f  2  J«r2(772,A:,^l)|^  J3,  ff  \Ur2?  ,2  •  ,  i, 

■-yY  '  ^  y  cos  (t> - — - dk  =  7rJJ^  -Y^kUmxpidr/pidk 
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The  integration  over  k  and  is  performed  numerically  at  every  streamwise  location. 
Note  that  the  results  depend  upon  the  three  dimensional  energy  spectrum,  E(k).  We 
assume  the  following  form  for  E{k): 


This  form  of  the  spectrum  was  used  by  Lee  et  al.  (1993,19946)  in  their  analysis  of 
the  interaction  of  isotropic  vortical  turbulence  with  a  shock  wave.  There  is  no  data  that 
suggests  the  chosen  spectral  dependency  for  an  isotropic  field  of  sound  waves.  However, 
only  the  inhomogeneous  part  of  the  flow-field  behind  the  shock  wave  is  dependent  upon 
the  energy  spectrum  of  the  incident  waves;  the  far-field  (772  — >■  00)  values  are  independent 
of  the  shape  of  the  spectrum.  We  choose  the  same  spectrum  as  Lee  et  al  to  allow  direct 
comparison  to  ther  results  of  the  inhomogeneous  near-field  behind  the  shock  wave.  Since 
the  inhomogeneous  component  of  the  flow-field  exponentially  decays  behind  the  shock 
wave,  it  is  the  far-field  that  is  of  importance  and  will  be  examined  in  detail. 

B.3  Results 

B.3.1  Spatial  variation  of  kinetic  energy 

The  evolution  of  the  kinetic  energy  behind  the  shock  wave  is  shown  in  Figure  B.5, 
where  =  Ra  is  plotted  as  a  function  of  downstream  distance.  The  inhomogeneous  nature 
of  the  velocity  field  is  apparent  immediately  downstream  of  the  shock.  Note  that  the  curves 
for  M  =  1.2  and  M  =  2  are  qualitatively  different.  While  decays  monotonically  for 
Af  =  1.2,  it  exhibits  rapid  nonmonotonic  variation  for  M  =  2.  A  similar  rapid  variation  of 
kinetic  energy  was  observed  in  the  interaction  of  vortical  fluctuations  with  a  shock  wave 
by  Lee  et  al  (1993).  However,  in  the  shock  wave/vortical  turbulence  interaction  problem, 
this  variation  was  observed  at  all  Mach  numbers  considered. 

Equations  (5)  and  (6)  are  used  to  show  that  this  variation  of  kinetic  energy  behind 
the  shock  wave  is  a  consequence  of  the  acoustic  waves  that  are  incident  at  <  V’l  <  V’cu- 
Recall  that  the  downstream  pressure  field  and  hence  the  velocity  field  set  up  by  these  waves 
decayed  exponentially.  Equation  (66)  shows  that  the  downstream  kinetic  energy  over  this 
regime  has  three  components:  a  homogeneous  component  associated  with  the  vorticity 
waves  (Ti),  a  monotonically  decreasing  component  due  to  acoustic  waves  (r2)  and  a  non¬ 
monotonic  component  due  to  the  correlation  between  the  vorticity  and  acoustic  waves  (Ts). 
Upon  integration  over  all  the  incident  waves,  these  components  combine  to  produce  the 
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rapid  evolution  immediately  downstream  of  the  shock.  As  seen  from  equations  (5a)  and 
(56),  the  kinetic  energy  associated  with  waves  incident  outside  this  range  has  homogeneous 
vortical  and  acoustic  components  and  a  nonmonotonic  inhomogeneous  component  due  to 
the  correlation  between  them.  When  integrated  over  all  incident  waves,  this  correlation 
term  is  however,  much  smaller  than  the  other  components.  This  is  illustrated  in  Figure 
B.6  where  the  kinetic  energy  is  decomposed  into  four  components  (vortical,  acoustic,  cor¬ 
relation  term  for  tpci  <  <  i’cu  and  correlation  term  outside  this  regime)  for  M  =  2.0. 

The  correlation  term  associated  with  waves  incident  at  ipci  <i>i  <  ^cu  is  seen  to  produce 
the  rapid  nonmonotonic  evolution  downstream. 


A  similar  decomposition  for  the  M  =  1.2  shock  showed  the  correlation  term  to  be 
insignificant  resulting  in  the  monotonic  evolution  of  Figure  B.5.  In  our  linear  analysis 
calculations,  the  nonmonotonic  downstream  variation  was  observed  for  shocks  whose  Mach 
nximber  exceeded  1.5.  In  studying  the  interaction  of  vortical  disturbances  with  the  shock 
wave,  Lee  et  al  (1993)  erroneously  concluded  that  linear  analysis  could  not  reproduce 
this  trend  which  they  had  found  in  the  DNS.  This  error  was  subsequently  corrected  and 
explained  by  them  (Lee  et  al,  1994a,19946). 


The  equation  governing  the  evolution  of  kinetic  energy  downstream  of  the  shock  wave 
provides  further  insight  into  the  spatial  evolution  of  kinetic  energy.  The  Euler  equations 
linearized  about  uniform  mean  fiow  may  be  rearranged  to  show  that  the  quantity. 


-litotal  — 
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1  + 
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^2p2 
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pu 

pZ 


(B.8) 


is  conserved  along  a  mean  streamline,  /total  changes  across  the  shock  wave  and  remains 
unchanged  downstream.  The  ^  correlation  determines  the  partitioning  of  /total  between 
potential  and  kinetic  energy.  The  spatial  uniformity  of  /total  and  the  exponential  decay  of 
behind  the  shock  wave  show  that  the  rapid  evolution  of  behind  the  shock  wave  is 
a  result  of  the  rapid  change  in  the  partitioning  of  /total  through  the  correlation  between 
pressure  and  the  shock-normal  component  of  fluctuating  velocity.  Decomposition  of  the 
pu  correlation  reveals  that  the  rapid  evolution  of  kinetic  energy  behind  the  shock  wave 
is  produced  by  the  waves  that  are  incident  at  tf^ci  <  i’l  <  ^cu-  The  far-field  values  of 
kinetic  energy  and  pressure  are  however,  determined  by  the  correlation  between  pressure 
fluctuations  and  the  acoustic  component  of  the  velocity  field  in  the  propagating  regime. 

B.3.2  Far-field  kinetic  energy 


As  seen  from  Figtire  B.5,  after  a  distance  that  is  comparable  to  the  lengthscale  (taken 
as  2‘K/ko)  of  the  incident  acoustic  waves,  the  kinetic  energy  asymptotes  to  its  far-field 
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value.  The  fax-field  values  axe  independent  of  the  upstreaxn  energy  spectrum  since  the 
inhomogeneous  terms  drop  out  in  the  fax-field.  As  a  result,  the  integration  over  k  can  be 
performed  independently  of  ^>1  and  (f)  to  give  f  E{k)  dk  =  5^/2.  To  gauge  the  effect  of 
shock  strength  upon  the  interaction,  we  examine  (Figure  B.7)  the  fax-field  kinetic  energy 
(normalized  with  the  upstream  kinetic  energy)  as  a  function  of  the  Mach  number  of  the 
shock  wave.  Note  that  the  shock-normal  component  of  kinetic  energy  is  larger  than  the 
transverse  components  for  all  Mach  numbers  shown.  An  interesting  feature  of  the  evolution 
of  kinetic  energy  is  that  kinetic  energy  decreases  slightly  across  the  shock  wave  over  a  range 
of  Mach  number  from  1.25  to  1.80.  The  transverse  components  decrease  across  the  shock 
wave  over  a  wider  range  of  the  Mach  number.  This  decrease  in  kinetic  energy  across  the 
shock  wave  is  peculiar  to  the  acoustic  wave/shock  interaction  problem.  It  is  not  observed 
in  the  interaction  of  vortical  fiuctuations  with  the  shock  wave  where  under  linear  analysis, 
kinetic  energy  increases  across  the  shock  wave  for  all  shock  strengths.  Also,  for  M  >  3, 
the  amplification  of  kinetic  energy  is  significantly  higher  than  amplification  levels  seen  in 
the  interaction  of  vortical  fluctuations  with  the  shock  wave. 

The  decrease  in  kinetic  energy  across  the  shock  wave  may  be  explained  by  decompos¬ 
ing  (Figure  B.8)  the  fax-field  energy  into  acoustic  and  vortical  components  and  examining 
their  dependence  on  Mach  number.  Both  components  are  normalized  with  the  upstream 
kinetic  energy.  This  decomposition  of  kinetic  energy  is  possible  since  the  correlation  be¬ 
tween  vortical  and  acoustic  components  falls  to  zero  in  the  fax-field.  Note  that  the  kinetic 
energy  associated  with  the  vortical  component  increases  monotonically  with  Mach  number 
and  exceeds  the  upstream  kinetic  energy  beyond  a  Mach  number  of  about  2.25.  The  kinetic 
energy  associated  with  the  acoustic  component,  however  decreases  across  the  shock  wave 
for  Mach  numbers  exceeding  about  1.2.  These  two  components  compete  in  determining 
the  overall  evolution  of  kinetic  energy.  For  lower  Mach  numbers,  the  acoustic  component 
dominates  since  not  enough  vorticity  is  generated  downstream,  causing  the  overall  kinetic 
energy  to  drop.  At  higher  Mach  munber,  the  vortical  component  dominates  due  to  in¬ 
creased  generation  of  vortical  fluctuations  and  the  overall  kinetic  energy  rises  across  the 
shock.  The  vortical  component  of  energy  exceeds  the  acoustic  component  for  Mach  number 
exceeding  2. 

The  decrease  of  the  fax-field  acoustic  kinetic  energy  with  Mach  number  is  explained 
as  follows.  As  the  Mach  number  increases,  *“  V’ci  increases  (see  Figure  B.4  of  Moore) 
and  hence  a  laxger  fraction  of  incident  waves  lie  in  the  range  xf^ci  <  V’l  <  V’cu-  Recall 
that  these  waves  make  no  contribution  to  the  fax-field  acoustic  kinetic  energy.  The  only 
contribution  to  the  far-field  acoustic  kinetic  energy  comes  from  waves  incident  outside  this 
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range  whose  amplification  increases  upon  increasing  the  Mach  number.  Upon  integrating 
over  all  incident  waves,  the  net  result  of  these  competing  factors  is  to  produce  (as  seen 
in  Figure  B.8)  a  slow  decrease  (for  M  <  3.5)  in  the  the  far-field  acoustic  kinetic  energy. 
Note  that  the  variation  of  ^ci  and  V’cw  with  Mach  number  is  most  rapid  at  lower  Mach 
numbers;  in  the  limit  of  infinite  Mach  number,  ^cZ  and  V’cu  axe  symmetrical  about  90°  and 
(for  7  =  1.4),  have  values  of  67.8°  and  112.2°  respectively. 

B.3.3  The  production  of  turbulence 

Having  observed  the  generation  of  vortical  fluctuations  -  ‘turbulence’  behind  the  shock 
wave,  we  examine  in  Figure  B.9,  the  intensity  of  the  turbulence  relative  to  the  intensity 
of  the  incident  pressxire  fluctuations.  Note  that  ($vort./^2)/(pi/-Pi)  of  flio  order  1  for 
most  of  the  Mach  numbers  shown.  (U2  =  V  -  17  is  the  mean  velocity  behind  the  shock 
wave  if  the  shock  were  stationary  in  the  mean)  This  suggests  an  interesting  possibility 
in  turbulent  flows  involving  multiple  shock  waves  (eg.  unadapted  supersonic  jets).  It 
is  known  that  upon  interaction  with  a  shock  wave,  turbulence  generates  intense  sound. 
For  isotropic  turbulence,  the  intensity  of  the  sound  generated  (P2/-P2)  scales  with  the 
intensity  of  the  incident  turbulence  (qi/Ui).  In  turbulent  flows  involving  multiple  shock 
waves,  it  is  reasonable  to  expect  the  generated  sotmd  waves  to  interact  with  the  subsequent 
shock  waves.  If  this  happens.  Figure  B.9  suggests  that  these  acoustic  waves  can  generate 
significant  levels  of  turbtilence  through  the  interaction.  We  do  not  mean  to  suggest  that  our 
homogeneous  analysis  is  directly  applicable  to  flows  involving  shock  cells,  which  axe  quite 
inhomogeneous.  However,  our  results  suggest  that  acoustic  wave/  shock  wave  interaction 
may  be  a  significant  generator  of  turbulence  in  these  flows. 


B.3.4  Thermodynamic  fluctuations 

Interaction  with  the  shock  wave  is  seen  to  significantly  increase  sound  levels.  Fig¬ 
ure  B.IO  shows  the  intensity  {P2IP2)  of  the  fax-field  pressure  fluctuations  normalized  with 
the  intensity  of  pressure  fluctuations  upstream  of  the  shock.  Note  that  the  intensity  drops, 
although  pressure  fluctuations  actually  amplify  across  the  shock  wave.  The  rise  across  the 
shock  wave  of  the  fax-field  soimd  pressure  level  (SPL)  and  the  level  of  acoustic  intensity 
(AIL)  in  decibels  is  shown  in  Figure  B.ll.  We  define, 

p^fRa 


SPL  =  10  log 


LPr^J 


AIL  =  10  log 


(B.9) 


/ Rr^r  J 

where  the  subscript  r  denotes  reference  values.  The  rise  in  the  level  of  sound  across 
the  shock  wave  is  independent  of  the  reference  values.  The  increase  in  sound  pressure 
level  varies  from  5  to  20  decibels  while  the  rise  in  acoustic  intensity  varies  from  2  to  10 
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decibels  over  the  range  of  Mach  numbers  shown.  Figure  B.IO  also  shows  the  nature  of  the 
thermodynamic  fluctuations  in  the  far-fleld.  The  density  and  temperature  fluctuations  are 
normalized  such  that  the  three  curves  would  collapse  if  the  fluctuations  were  isentropic. 
We  see  that  the  isentropic  relations  hold  \mtil  a  Mach  number  of  1.5,  beyond  which  the 
entropy  fluctuations  that  axe  generated  at  the  shock  wave  become  significant  relative  to 
the  acoustic  fluctuations.  The  increasing  importance  of  the  entropy  fluctuations  in  the 
far-fleld  is  due  to  two  factors:  increased  production  of  entropy  fluctuations  at  the  shock 
wave  and  the  decrease  with  Mach  nvunber  of  the  fax-field  intensity  of  pressure  fluctuations. 
Asymptotically,  (s2/C'p)/(p2/F2)  equals  1.47  (7  =  1.4)  where  S2  is  the  rms  entropy  behind 
the  shock  wave  and  Cp  is  the  specific  heat  at  constant  pressmre. 

B.3.5  Asymptotic  behaviour 

Another  feature  that  distinguishes  the  interaction  of  sound  waves  with  a  shock  wave 
is  the  asymptotic  behavior  with  respect  to  Mach  number.  Linear  analysis  shows  that  for 
incident  acoustic  waves,  quantities  such  as  the  ratio  of  fluctuating  kinetic  energy,  pressure 
and  temperature  fluctuations  across  the  shock  wave  are  not  bounded  but  increase  as  M^. 
This  is  easily  seen  from  equations  (5)  and  (6).  The  constants  {eg.  K)  in  the  equations 
have  finite  values  in  the  limit  of  infinite  Mach  number  and  hence  the  intensities  of  velocity, 
density  and  pressure  fluctuations  axe  bounded  for  finite  intensity  of  incident  pressure  fluc¬ 
tuations.  However,  for  a  given  intensity  of  pressure  fluctuations  the  incident  kinetic  energy 
decreases  as  1/M^  causing  the  amplification  of  kinetic  energy  to  vary  as  for  strong 

shocks.  Similarly,  since  the  mean  pressure  and  temperature  ratio  varies  as  for  strong 

shocks,  the  ratio  of  pressure  and  temperature  fluctuations  axe  not  bounded.  Unbound¬ 
edness  of  kinetic  energy  amplification  is  peculiar  to  the  acoustic  wave/shock  interaction 
problem;  it  is  not  present  when  vortical  fluctuations  interact  with  the  shock  wave.  We 
show  in  Figure  B.12  the  kinetic  energy  amplification  normalized  by  M^.  Note  that  the 
amplification  ratios  asymptote  to  finite  values  when  properly  scaled  with  . 

B.4.  Discussion 

B.4.1  Evaluation  of  linear  analysis 

The  inviscid  linear  analysis  of  Moore  (1954)  is  extended  to  study  the  interaction  of  a 
three-dimensional  isotropic  field  of  acoustic  waves  with  a  normal  shock  wave.  The  objective 
of  this  study  is  to  isolate  the  effect  of  acoustic  waves  on  the  evolution  of  a  turbulent  flow  as 
it  interacts  with  a  shock  wave.  As  the  level  of  compressibility  of  a  turbulent  flow  increases, 
the  effects  associated  with  the  acoustic  component  will  become  important.  Understanding 
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the  interaction  of  acoustic  waves  with  a  shock  wave  is  therefore  of  fundamental  importance. 
The  interaction  is  likely  to  be  especially  important  in  unbounded  turbulent  flows. 

Our  use  of  linear  analysis  is  prompted  by  its  success  in  past  investigations.  In  the 
interaction  of  a  single  unsteady  disturbance  with  a  shock  wave,  there  appear  to  be  four  im¬ 
portant  factors  that  could  cause  deviation  from  linear  behavior:  proximity  of  the  incident 
angle  to  the  critical  angle,  the  amplitude  of  the  disturbance  being  non-negligible  relative 
to  the  strength  of  the  shock  wave,  the  mean  Mach  number  being  close  to  imity  and  finally, 
viscous  effects.  Proximity  to  the  critical  angle  seems  to  be  the  most  important  factor.  Our 
computations  of  of  single  acoustic  waves  interacting  with  a  shock  wave  disagree  with  linear 
analysis  arotmd  25  degrees  of  the  critical  angle;  very  good  agreement  is  seen  outside  this 
range.  Fortunately,  in  the  interaction  of  a  spectrum  of  waves  with  a  shock  wave,  the  inci¬ 
dent  waves  aroimd  the  critical  angle  make  only  partial  contribution  to  the  overall  energy. 
With  regard  to  the  amplitude  of  incident  disturbances,  results  of  computations  are  encour¬ 
aging.  Our  computations  show  good  agreement  with  linear  analysis  away  from  the  critical 
angle  for  disturbance  amphtudes  varying  jhom  0.14%  to  14%.  Zang  et  al  (1984)  show  that 
linear  analysis  yields  good  prediction  away  from  the  critical  angle  for  acoustic  disturbance 
amplitudes  as  large  as  25%.  Meadows  et  al  (1994)  reach  a  similar  conclusion  in  their  study 
of  the  one-dimensional  interaction  of  an  acoustic  wave  with  a  shock  wave.  Their  compar¬ 
ison  of  hnear  analysis  to  a  non-linear  Riemann  analysis  yielded  indistinguishable  results 
(to  plotting  accuracy)  for  disturbance  amplitudes  less  than  10%.  An  effect  of  increasing 
disturbance  amplitude  for  a  fixed  shock  strength  has  been  docximented  by  Honkan  and 
Andreopoulous  (1992)  and  Lee  et  al  (1992)  in  the  interaction  of  isotropic  turbulence  with 
a  normal  shock.  Both  studies  reported  a  decrease  in  the  amplification  of  kinetic  energy 
as  the  disturbance  amphtude  was  increased.  Honkan  and  Andreopoulous’  experimental 
study  had  a  shock  of  mean  Mach  number  1.24  while  Lee  et  aFs  numerical  results  were  for 
a  shock  wave  of  mean  Mach  number  1.2. 

B.4.2  Combined  interaction  of  vortical  and  acoustic  fluctuations 

The  linearity  of  the  analysis  permits  the  use  of  superposition  to  predict  the  combined 
interaction  of  vortical  and  acoustic  fluctuations  with  a  shock  wave.  Since  soimd  and 
vorticity  travel  at  different  speeds,  the  correlation  between  them  can  be  assumed  neghgible 
upstream  and  in  the  far-field  behind  the  shock  wave.  The  amplification  of  kinetic  energy 
across  the  shock  wave  can  therefore  be  expressed  as. 


g"  ^  ql  +  ql 

go  (g«)o  +  (g^)o 


(B.IO) 
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where  ,  the  subscript  0  represents  conditions  upstream  of  the  shock  and  the  subscripts  uj 
and  6  represent  the  vortical  and  dilatational  components  respectively.  The  amplification 
ratio  in  the  mixed  problem  may  be  expressed  in  terms  of  the  amplification  ratios  of  vortical 
and  acoustic  fluctuations  as:  ^ 

^  =  (i-X)u+xu  (B.m 

where,  fu;  =  9^/(95)o5  fo  =  Qel{Qe)o  ^  ~  (Qe)ofQo  •  ^  the  ratio  of  acoustic  to  total 
kinetic  energy  upstream  of  the  shock  wave.  Note  that  as  X  varies  from  0  to  1,  5^ /q2 
from  fu  to  fe.  Available  experiments  on  the  shock/turbulence  interaction  are  in  the  range 
of  Mach  numbers  for  which  the  kinetic  energy  of  the  acoustic  component  decreases  across 
the  shock  wave.  Our  results  suggest  that  the  presence  of  acoustic  waves  in  these  experi¬ 
ments  will  decrease  the  amplification  of  kinetic  energy.  As  noted  by  Jacquin  et  al  (1993), 
wind-tunnel  experiments  on  the  shock/turbulence  interaction  report  lower  amplification 
ratios  than  shock-tube  experiments.  They  suggest  that  boundary-layer  induced  oscillation 
of  the  shock  wave  in  the  wind-tunnel  experiments  might  be  responsible  (through  an  un¬ 
known  mechanism)  for  this  observation.  Our  results  suggest  the  alternative  scenario  that 
the  lower  kinetic  energy  amplification  is  a  consequence  of  a  significant  amount  of  acoustic 
waves  in  the  turbulence  ahead  of  the  shock  wave.  A  hkely  source  of  these  sound  waves  are 
the  shock  waves  that  are  generated  at  the  edges  of  the  grid  that  generates  turbulence  in 
the  tunnel.  Further  data  is  needed  to  support  or  discard  this  scenario. 

According  to  hnear  analysis  (equation  18),  when  compressible  isentropic  turbulence 
interacts  with  a  shock  wave,  the  evolution  of  kinetic  energy  lies  between  two  limits  -  the 
pure  solenoidal  and  pure  dilatational  limits.  This  behavior  is  identical  to  the  response 
of  homogenous  turbulence  to  one- dimensional  compression.  As  shown  by  Cambon  et  al, 
the  evolution  of  kinetic  energy  when  compressible  homogeneous  turbulence  is  subjected 
to  one-dimensional  compression  lies  between  the  solenoidal  and  the  dilatational  limits. 
However  an  important  difference  between  the  shock/turbulence  interaction  problem  and 
the  homogeneous  compression  problem  is  the  behavior  of  the  tiirbulence  in  the  dilatational 
limit.  For  homogeneous  turbulence,  Cambon  et  al  show  that  the  dilatational  limit  is 
characterized  by  loss  of  acoustic  communication  which  forces  the  pressure-strain  correlation 
to  zero  thereby  yielding  significantly  higher  amplification  of  as  compared  to  solenoidal 
turbulence.  The  compressibility-induced  decrease  in  kinetic  energy  across  the  shock  wave 
that  is  seen  in  the  acoustic  wave/shock  wave  interaction  for  1.25  <  M  <  1.8  is  not  observed. 
This  is  due  to  the  fundamentally  different  nature  of  the  pressure  fluctuations  in  the  two 
problems. 

The  difference  between  the  one- dimensional  compression  of  solenoidal  turbulence  and 
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Ribner’s  (1954,1987)  analysis  of  solenoidal  fluctuations  interacting  with  a  shock  wave  were 
pointed  out  by  Lee  et  al  (1993)  and  underscored  by  Jacquin  et  al  (1993).  Results  of 
the  two  problems  agree  for  small  Mach  numbers;  for  large  Mach  numbers  significantly 
larger  amplification  is  seen  in  the  homogeneous  problem.  Jacquin  et  al  also  compare  the 
compression  of  homogeneous  compressible  turbulence  to  Ribner’s  analysis  and  note  that 
the  difference  between  the  two  problems  is  greater  than  pointed  out  by  Lee  et  al.  A 
compaxison  of  the  compression  of  dilatational  fluctuations  to  the  interaction  of  solenoidal 
fluctuations  with  a  shock  wave  is  not  appropriate;  comparison  should  be  made  to  the 
interaction  of  dilatational  fluctuations  with  the  shock  wave.  Such  comparison  would  reveal 
disagreement  in  the  dilatational  limit  over  the  entire  range  of  Mach  numbers.  While  the 
homogeneous  problem  can  approximate  the  evolution  of  kinetic  energy  in  the  interaction 
of  solenoidal  fluctuations  with  a  shock  wave  of  moderate  strength,  it  is  quite  inappropriate 
as  the  incident  turbulence  becomes  increasingly  compressible. 

B.4.3  Comparison  to  computation 

The  analysis  explains  recent  observations  (Hannapel  and  Friedrich,  1994)  on  the  in¬ 
teraction  of  compressible  turbulence  with  a  normal  shock.  Hannapel  and  Friedrich  nu¬ 
merically  computed  the  interaction  of  low  Reynolds  number  isotropic  turbulence  (R\  ~  4 
upstream  of  the  shock  wave)  with  a  normal  shock  of  mean  Mach  number  2.  The  fluctuating 
Mach  number  was  0.1  at  the  inflow  and  the  kinetic  energy  was  equally  distributed  between 
the  vortical  and  acoustic  modes.  {X  =  0.5  in  our  notation.)  They  compared  this  interac¬ 
tion  with  that  of  essentially  solenoidal  turbulence  {X  =  0)  of  the  same  fluctuating  Mach 
number.  The  ‘compressible’  case  displayed  increased  amplification  of  vorticity,  decreased 
amplification  of  velocity  components  transverse  to  the  shock,  decreased  reduction  of  Taylor 
microscale  and  decreased  amplification  of  density,  temperatme  and  pressure  fluctuations 
across  the  shock  wave. 

The  computation  may  suffer  from  lack  of  sufficient  resolution  of  the  shock  front.  As  a 
result  we  do  not  attempt  quantitative  comparison.  We  show  that  the  computation  follows 
the  trends  predicted  by  the  linear  analysis  and  quote  the  linear  analysis  predictions  for 
reference.  The  higher  ampHfication  of  the  transverse  components  of  vorticity  in  the  com¬ 
pressible  case  is  explained  by  the  generation  of  vorticity  through  the  acoustic  wave/shock 
wave  interaction.  Using  an  expression  similar  to  equation  (18),  the  increase  in  the  ampli¬ 
fication  of  vorticity  may  be  shown  to  be  3^@(^o^/‘^o^))  where  gg  =  IOq^  in  the  acoustic 
wave/shock  wave  interaction  problem.  Similarly,  linear  analysis  predicts  the  decrease  in 
transverse  velocity  and  reduction  in  Taylor  microscale.  For  the  interaction  of  vortical 
fluctuations  with  a  shock  wave,  Ribner’s  analysis  yields  v? {uq  =  1.64  and  =  1-66 


across  a  shock  wave  of  Mach  number  2.  For  dilatational  fluctuations,  our  analysis  yields 
values  of  1.69  and  0.9  respectively.  Using  an  expression  similar  to  equation  (18),  we  get 
jul  =  1.67  and  =  1-28;  i.e.,  the  amplification  of  is  essentially  unchanged  while 
the  amplification  of  the  transverse  velocity  drops. 

Linear  analysis  yields  values  of  0.53  and  0.77  for  Ai/(Ai)o  and  A2/(A2)o  respectively 
(Aj  represents  the  Taylor  microscale  in  the  i  direction)  for  the  vortical  problem.  Corre¬ 
sponding  values  for  the  compressible  problem  are  0.66  and  0.94  respectively.  The  reduced 
amplification  of  thermodynamic  fluctuations  in  the  compressible  problem  (an  order  lower) 
is  a  consequence  of  the  fact  that  in  pure  vortical  turbulence,  the  absolute  level  of  incident 
thermodynamic  fluctuations  is  much  lower  (zero  in  the  hnear  limit).  The  level  of  ther¬ 
modynamic  fluctuations  behind  the  shock  wave,  however  scale  with  the  incident  kinetic 
energy  as  a  result  of  which  the  amplification  of  thermodynamic  fluctuations  will  be  quite 
large  (undefined  within  linear  analysis). 

B.5.  Summary 

Inviscid  linear  analysis  was  used  to  study  the  evolution  of  fluctuating  kinetic  energy, 
sound  level  and  thermodynamic  fluctuations  in  the  interaction  of  an  isotropic  field  of 
acoustic  waves  with  a  normal  shock  wave.  The  analysis  was  an  extension  of  Moore’s 
(1954)  study  of  the  interaction  of  a  shock  wave  with  a  single  acoustic  wave.  Moore’s 
analysis  was  evaluated  by  comparison  to  numerical  computation.  Good  agreement  was 
seen  for  angles  away  from  the  critical  angles.  A  possible  reason  for  disagreement  around 
the  critical  angles  was  proposed.  Linear  approximation  was  noted  to  be  inconsistent  in 
this  regime;  requirement  of  the  transonic  small  disturbance  equations  was  suggested. 

The  interaction  of  an  isotropic  acoustic  field  with  the  shock  was  noted  to  be  sig¬ 
nificantly  different  from  that  of  vortical  fluctuations.  The  kinetic  energy  of  the  acoustic 
fluctuations  decreased  across  the  shock  wave  for  Mach  numbers  between  1.25  and  1.8.  For 
Mach  numbers  exceeding  3,  the  kinetic  energy  amplified  by  levels  that  significantly  ex¬ 
ceeded  those  found  in  the  interaction  of  vortical  fluctuations  with  the  shock.  These  trends 
were  explained  by  decomposing  the  velocity  field  into  acoustic  and  vortical  components 
and  examining  their  dependence  on  the  Mach  number.  Upon  interacting  with  the  shock 
wave,  the  acoustic  waves  generated  vortical  fluctuations  whose  contribution  to  the  far-field 
kinetic  energy  increased  with  increasing  Mach  number.  The  level  of  sound  increased  across 
the  shock  wave.  The  rise  in  the  sound  pressure  level  across  the  shock  varied  from  5  to 
20  decibels  for  Mach  number  varying  from  1.5  to  5.  The  fluctuations  behind  the  shock 
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wave  were  nearly  isentropic  for  Mach  number  less  than  1.5  beyond  which  the  generation 
of  entropy  fluctuations  became  significant. 

Finally  the  analysis  was  used  to  describe  the  combined  interaction  of  isotropic  vortical 
and  acoustic  fluctuations  with  a  normal  shock  wave. 


B.6  Details  of  Moore’s  analysis 


The  solution  dowstream  of  the  shock  wave  is  given  below.  The  subscript  2  refers  to 
conditions  behind  the  shock  wave. 


For  0  <^i  <  'tpci  or  ^cu  <  V’l  ^ 


- a; - ) 


The  X2  aods  is  stationary  in  a  frame  of  reference  that  moves  with  speed  U  yielding,  X2  = 
(V  —  U)t  at  the  mean  position  of  the  shock  wave.  Of  the  two  terms  that  contribute  to 
the  velocity,  density  and  temperature,  the  unsteady  term  corresponds  to  the  acoustic  wave 
while  the  steady  term  corresponds  to  the  vorticity  and  entropy  waves.  The  coefl5cients 
K,  Q,  G,  H  and  I  in  the  above  equations  are  functions  of  the  Mach  number  of  the  shock 
and  the  angle  of  incidence  and  are  given  below. 

If  ‘tpci  <  tpi  <  tfle  solution  is  a  bit  more  complex  and  is  given  by: 
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1  V2 


Az  V 


—  =  ff(l)^(l)(v,0  +  -^(2)^(2)(v, 0  +  ^(1)/ 


+  ■^(2)5'  I 


=^X2-(!/' 


Ai  )  ■  V  A, 

Analogous  to  the  propagating  regime,  <5(a)5  G(a),  if(Q)  and  /(q)  are  functions  of  the 
Mach  number  of  the  shock  and  the  angle  of  incidence  and  axe  given  below.  Also,  rj  and  ^ 
are  related  to  X2,y  and  t  by: 


-?  =  -T-l.^2  -{V-  V)i] 


(  =  -^(oiX2  +  0y  +  cVt). 

T]  is  thus  proportional  to  the  distance  behind  the  shock  while  ^=constant  represents  an 
oblique  plane  moving  at  constant  velocity.  The  functions  g,  and  $(2)  aJ’e  defined  as: 
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772  4-  (r  -  77)2 

Since  a  homogenous  acoustic  field  may  be  represented  as  a  superposition  of  plane  waves, 
we  consider  f{x)  =  e*®  for  which  y,$(i)  and  #(2)  are  given  by: 


g{i)  =  i  e‘^  $(i)  =  e  ’’e*^  $(2)  =  i 


Note  that  if  the  incident  wave  is  a  plane  wave,  the  above  equations  take  on  a  simple 
form.  The  coefficients  in  the  solution  are  determined  as  follows. 
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For  0  <  V’l  <  fl’ci  or  V’ca  <  V’l  <  Tf, 
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For  ^c/  ^  ^  i’cul  the  coefficients  axe  as  follows. 
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The  expression  for  H(^i)  in  Moore’s  paper  contains  a  typographical  error;  L(i)  +  ^ 

should  be  replaced  by  +  +  Also,  the  plots  of  the  constants  when  ^  V’cu 

do  not  agree  with  the  formulae.  After  repeating  Moore’s  analysis,  we  conclude  that  the 
formulae  are  correct. 
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FIG.  B.2.  Comparison  of  the  predictions  of  the  linear  analysis  to  computed  results.  The  lines 
are  values  obtained  from  analysis  while  the  symbols  represent  computed  values.  All 

quantities  are  normalized  by  pi/Pj.  - (^2/^1), - (<^2/^1).  ^2  and  U2  denote 

the  rms  dilatation  and  vorticity  behind  the  shock  wave  respectively. 
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FIG.  B.3.  Computed  values  of  the  vorticity  behind  the  shock  wave  as  a  function  of  the  amplitude 
of  the  incident  disturbance.  The  mean  Mach  number  is  1.5.  o  (^i  =  100°),  x  (^i  = 
60°) 
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FIG.  B.4.  Coordinate  system  used  in  analysis  of  the  interaction  of  a  shock  wave  with  an  isotropic 
field  of  acoustic  waves  that  axe  incident  from  upstream. 
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FIG.  B.5.  Evolution  of  behind  the  shock  wave  as  predicted  by  linear  analysis,  is  normalized 
with  its  upstream  value. - (M  =  1.2), - (M  =  2.0). 


FIG.  B.6.  Decomposition  of  q^  behind  the  shock  wave  under  linear  analysis.  The  Mach  number 

is  2  and  all  components  are  normalized  with  the  upstream  value  of  q^ .  - (To- 

tal),  (Vortical),  .  (Acoustic), - (Correlation  for  ^c/  <  V’l  <  V’ca), 

- (Correlation  for  0  <  <  V”!  <  ). 
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M 

FIG.  B.7.  Fax-field  kinetic  energy  as  a  function  of  Mach  number.  All  components  are  normalized 
with  their  upstream  value. - . 
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FIG.  B.8.  The  far-field  kinetic  energy  decomposed  into  acoustic  and  vortical  components.  Both 
components  are  normalized  with  the  upstream  value  of  .  (Total), - (Vor¬ 
tical),  - (Acoustic). 
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M 

FIG.  B.9.  The  far-field  ‘tiirbulent’  intensity  compared  to  the  intensity  of  incident  pressure  fluc¬ 
tuations.  The  curve  asymptotes  to  2.75  (7  =  1.4). 


FIG.  B.IO.  The  intensity  of  pressure,  density  and  temperature  fluctuations  in  the  far-fleld  com¬ 
pared  to  the  intensity  of  incident  pressure  fluctuations. - (p2 /F2), - 'y(p2/R2), 

.  7/(7~1)(^2/T2).  62  and  T2  are  the  rms  and  mean  temperature  behind  the  shock 

wave  respectively.  The  three  curves  are  normalized  with  pi/Pi  and  asymptote  to  0.29, 
0.67  and  1.54  respectively  (7  =  1.4). 
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Abstract 

Direct  numerical  simulation  (DNS)  and  linear  anal¬ 
ysis  (LIA)  of  isotropic  turbulence  interacting  with  a 
shock  wave  are  performed  for  several  upstream  shock 
normal  Mach  numbers  (Mi).  Turbulence  kinetic  energy 
(TKE)  is  amplified  across  the  shock  wave,  but  this  am¬ 
plification  tends  to  saturate  beyond  Mi  =  3.0.  TKE 
amplification  and  Reynolds  stress  anisotropy  obtained 
in  DNS  are  consistent  with  LIA  predictions.  Rapid 
evolution  of  TKE  immediate  downstream  of  the  shock 
wave  persists  for  all  shock  strengths  and  is  attributed  to 
the  transfer  between  kinetic  sind  potential  modes  of  tur¬ 
bulence  energy  through  acoustic  fluctuations.  Changes 
in  energy  spectra  and  various  length  scales  across  the 
shock  wave  are  predicted  by  LIA,  which  is  consistent 
with  DNS  results.  Most  turbulence  length  scales  de¬ 
crease  across  the  shock.  Dissipation  length  scale  (pq^/f), 
however,  increases  slightly  for  shock  waves  with  Mi  < 
1.65.  Fluctuations  in  thermodynamic  variables  behind 
the  shock  wave  stay  nearly  isentropic  for  Mi  <  1.2, 
and  deviate  significantly  from  isentropy  for  the  stronger 
shock  waves  due  to  large  entropy  fluctuation  generated 
through  the  interaction. 

1.  Introduction 

Understamding  the  mechanisms  of  turbulence  inter¬ 
acting  with  a  shock  wave  is  of  fundamental  importance 
in  predicting  the  interaction  of  turbulent  boundary  lay¬ 
ers  with  the  shock  waves  which  occur  in  many  engineer¬ 
ing  applications.  Since  1950’s,  lineaur  analyses  (LIA) 
on  the  modification  of  elementary  disturbance  waves, 
such  as  vorticed,  acoustic,  and  entropic  waves,  by  the 
shock  wave  have  been  performed  with  an  emphasis  on 
the  acoustic  wave  generation  behind  the  shock  wave^~^. 
Recently,  the  applicability  of  homogeneous  Rapid  Dis¬ 
tortion  Theory  (RDT)  on  shock/turbulence  interaction 
is  investigated  by  Jacquin  ei  al.^ 

There  has  been  a  significant  accumulation  of  ex¬ 
perimental  data  on  shock  turbulence  interaction  dur¬ 
ing  the  last  decade.  Interaction  of  turbulent  bound¬ 
ary  layers  with  a  shock  wave  over  a  corner  was  in¬ 
vestigated  by  many  research  groups,  among  them  are 


Dolling  and  Or®,  Andreopoulos  and  Muck^°,  and  Smits 
and  Muck^^.  A  general  finding  from  these  experiments 
is  that  Reynolds  shear  stress  and  turbulence  intensities 
are  amplified  across  the  shock  wave.  The  studies  of 
oblique  shock  wave/turbulent  boundary  layer  interac¬ 
tion  included  several  additional  phenomena  which  com¬ 
plicated  the  flow  behavior To  isolate  the  effects  of 
a  shock  wave  on  turbulence,  several  experiments^^” 
on  the  interaction  between  the  shock  wave  amd  grid¬ 
generated  turbulence  have  been  performed.  An  impor¬ 
tant  finding  is  that  turbulence  is  amplified  Euid  turbu¬ 
lence  length  scales  increase  across  the  shock  wave.  But 
the  length  scale  increase  contradicts  the  intuitive  expec¬ 
tation  that  mean  flow  compression  should  decrease  the 
relevant  turbulence  length  scales.  The  issue  of  length 
scale  change  will  be  thoroughly  discussed  in  the  present 
paper  (Sec.  3.2). 

Numerical  simulations  of  the  shock  turbulence  inter¬ 
action  aire  just  beginning  to  emerge.  Using  a  shock  cap¬ 
turing  numerical  technique,  Rotman^®  calculated  the 
change  in  a  two  dimensional  turbulent  flow  caused  by 
the  passage  of  a  travelling  shock  wave.  He  found  that 
the  shock  causes  an  increase  in  the  turbulent  kinetic 
energy  and  that  turbulence  length  scales  are  reduced 
upon  passage  of  the  shock.  Lee  et  conducted  di¬ 

rect  numerical  simulations  of  two  and  three  dimensional 
turbulence  interacting  with  a  shock  wave.  They  found 
that  vorticity  eunplification  comp2ued  well  with  the  lin¬ 
ear  analysis  predictions,  and  turbulent  kinetic  energy 
undergoes  rapid  increase  behind  the  shock  wave.  The 
spectrum  was  found  to  be  enhanced  more  at  large  wave 
numbers,  leetding  to  an  overall  length  scade  decrease. 

In  the  present  paper,  interaction  of  isotropic  turbu¬ 
lence  with  a  strong  shock  wave  is  studied  to  investigate 
the  effects  of  shock  strength  on  turbulence  modifica¬ 
tion.  We  developed  a  numerical  technique  to  simu¬ 
late  turbulence  interacting  with  a  strong  shock  wave 
without  resolving  its  structure,  eind  validated  this  tech¬ 
nique  against  the  shock-resolving  simulations  (Sec.  2). 
The  simulation  results  are  compared  with  the  results 
from  linear  anadysis,  amd  are  contrasted  against  the  re¬ 
sults  from  the  weak  shock  case  to  highlight  the  shock 
strength  effects  (Sec.  3). 
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2.  Numerical  Simulation  and  Linear  Theory 
2.1  Shock-Capiuring  Scheme 

For  the  direct  numerical  simulation  (DNS)  of  the 
interaction  of  isotropic  turbulence  with  a  weak  shock 
wave,  the  shock  wave  structure  was  resolved  as  a  solu¬ 
tion  of  the  Navier-Stokes  equations.  In  addition,  all  the 
essential  turbulence  scales  were  resolved^®.  For  a  strong 
shock  wave,  however,  resolving  the  shock  structure  is 
irrelevant  because  the  Navier-Stokes  equations  are  no 
longer  valid  inside  the  shock  wave  for  a  shock  with  a 
shock  normal  Mach  number  exceeding  2.0®°.  Therefore, 
for  an  efficient  computation,  a  shock-capturing  scheme 
is  to  be  implemented  for  the  simulation  with  a  strong 
shock  wave.  The  shock-capturing  scheme  is  required  to 
be  high  order  accurate  throughout  the  computational 
domain  (to  properly  simulate  the  evolution  of  turbu¬ 
lence)  as  well  as  to  be  able  to  give  smooth  and  accurate 
transition  across  the  shock  wave.  The  shock-capturing 
scheme  chosen  was  an  essentially  non-oscillatory  (ENO) 
scheme,  which  can  be  constructed  up  to  any  order  of 
accuracy.  The  ENO  scheme  used  in  this  paper  was 
based  on  the  Lsuc-Friedrichs  scheme  with  an  interpola¬ 
tion  on  fluxes,  and  for  two  and  three  dimensional  cases 
the  ENO  scheme  is  applied  independently  for  fluxes  in 
each  direction®^.  Some  modifications  described  in  the 
following  were  made  to  the  basic  scheme  to  improve 
the  solution  accuracy  and  to  enhance  the  code  perfor¬ 
mance. 

The  ENO  scheme  was  implemented  in  the  code  de¬ 
veloped  for  the  simulation  of  shock  turbulence  inter¬ 
action.  Simulation  of  spatially  decaying  compressible 
turbulence  without  shocks  was  performed  to  validate 
the  scheme  and  to  determine  the  required  order  of  the 
ENO  scheme  to  reproduce  the  evolution  of  turbulence 
as  determined  from  a  corresponding  DNS.  During  the 
validation  procedure,  an  unexpected  accuracy  degener¬ 
acy  was  observed  with  the  increase  in  the  ENO  scheme’s 
order.  As  the  order  of  the  ENO  scheme  increases,  the 
solution  deviates  further  from  the  DNS  result.  Similar 
accuracy  degeneracy  was  reported  when  the  mesh  was 
refined®®.  This  degeneracy  was  due  to  the  choice  of 
linearly  unstable  stencils.  Biasing  toward  a  central  dif¬ 
ference  stencil  was  implemented  to  avoid  the  accuracy 
degeneracy  in  a  shock-free  region®^.  With  this  fix,  more 
accurate  results  could  be  obtained  by  raising  the  order 
of  accuracy.  Results  comparable  to  the  corresponding 
DNS  (with  a  sixth-order  accurate  Fade  scheme)  was 
obtained  with  a  sixth-order  ENO  scheme. 

When  the  ENO  scheme  is  implemented  throughout 
the  domain,  the  operation  count  and,  accordingly,  the 
CPU  time  increase  significantly  compared  to  the  Fade 
scheme.  The  region  where  the  ENO  scheme  plays  a 


significant  role  is  quite  localized,  usually  less  than  a 
tenth  of  the  computational  domain  in  DNS  of  weak 
shock  waves,  and  is  expected  to  be  even  smaller  for 
strong  shock  waves.  The  idea  of  applying  the  ENO 
scheme  only  where  it  plays  an  active  role  2md  switch¬ 
ing  elsewhere  into  the  usual  Fade  central  differencing 
was  tested  to  save  computer  resources.  The  local  re¬ 
gion  of  ENO  application  can  even  be  specified  a  priori 
to  be  a  zone  around  the  shock  for  a  simulation  per¬ 
formed  in  a  coordinate  system  fixed  on  the  mean  shock 
wave.  The  concept  of  local  application  of  the  ENO 
scheme  was  validated  against  fully  resolved  DNS  results 
both  for  the  case  of  two  dimensional  shock  turbulence 
interaction^®  with  Mi  =  1.2,  M*  =  0.07  and  the  case  of 
three  dimensional  spatially  decaying  turbulence®'*  with 
Mt  =  0.51  eind  Rex  =  25.  Here,  Mi  =  «i/c,  Mt  =  q/c, 

and  Rex  =  purmsA//!,  where  Unns  = 
c  is  the  sound  speed,  A  is  a  Taylor  microscale,  and  fi 
is  the  dynamic  viscosity.  The  overline  denotes  an  en¬ 
semble  average,  and  /'  and  f"  are  the  deviations  from 
^e  ensemble  average  and  the  mass-weighted  average®® 
f  =  Pf/'P^  that  is,  /'  =  /  -  7  and  f"  =  f-  f. 

The  next  validation  computation  was  to  reproduce 
three  dimensional  shock  turbulence  interaction  using 
the  ENO  scheme.  Significant  additional  dissipation 
of  turbulent  kinetic  energy  (TKE)  was  found  in  the 
shock  vicinity  where  the  ENO  scheme  was  applied  in¬ 
dependently  in  all  directions,  which  was  not  observed 
in  two  dimensional  shock  turbulence  simulations.  This 
suggested  that  the  ENO  scheme  be  applied  only  in 
the  shock  normal  direction.  With  the  one-directional 
ENO  scheme  application,  more  accurate  results  were 
obtained  with  less  CFU  time  used.  For  Mi  =  1.2,  the 
ENO  scheme  required  only  a  quarter  of  the  CFU  time 
in  the  corresponding  shock-resolving  simulation®®. 

The  final  validation  calculation  was  made  for  two 
dimensional  turbulence  interacting  with  a  strong  shock 
wave.  Results  from  the  ENO  scheme  on  two  dimen¬ 
sional  shock  turbulence  interaction  were  compared  with 
the  corresponding  shock  resolving  simulation,  where 
Ml  =  2.0  and  Mt  =  0.07  at  the  inflow.  The  minimum 
grid  spacing  required  by  the  ENO  scheme  was  five  times 
larger  than  that  used  in  the  shock-resolving  simulation. 
Figure  1  compares  the  evolution  of  velocity  fluctuation 
variances,  where  inflow  turbulence  is  generated  using 
the  energy  spectrum  of  E{k)  ~  ul{k/k„Y  exp[-2(it//to)®] 
with  ko  the  peak  wave  number  and  the  rms  velocity 
fluctuation.  The  amplification  and  evolution  behind 
the  shock  wave  were  reproduced  by  the  ENO  scheme 
with  six  times  larger  computational  mesh  near  the  shock 
wave.  In  order  to  accurately  predict  the  interaction  of 
turbulence  with  stronger  shock  waves,  where  a  shock- 
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resolving  simulation  is  not  practical,  reliable  results  can 
only  be  obtained  after  a  thorough  grid-refinement  test. 
Grid  refinement  is  necessary  to  predict  thermodynamic 
fluctuations  correctly,  since  the  convergence  of  the  ther¬ 
modynamic  variables  is  slowest  due  to  spurious  numer¬ 
ical  oscillations  behind  the  shock  wave  in  the  shock¬ 
capturing  scheme^^’^®.  Details  of  the  validation  tests 
have  been  reported  in  ref.  26. 


Figure  1 .  Evolution  of  the  velocity  fluctuation  vari¬ 
ances  {Ml  =  2.0,  Mt  =  0.07):  lines  denote  results  from 
shock  resolving  simulation  and  symbols  from  shock  cap¬ 
turing  simulation. - ,  •  :  iZu, - ,  x  :  Ri2- 

Numerical  procedures  used  in  this  paper  can  be  sum¬ 
marized  as  follows:  The  fluid  is  assumed  to  be  an  ideal 
gas  with  the  specific  heat  ratio  of  7  =  1.40.  The  Prandtl 
number  is  assumed  constant  with  Pr  =  0.75.  The 
density,  three  components  of  momentum,  and  the  to¬ 
tal  internal  energy  are  adv£inced  in  time  using  a  third- 
order  compact-storage  Runge-Kutta  scheme^®.  Spatial 
derivatives  are  approximated  by  the  sixth  order  ENO 
scheme  only  for  the  inviscid  fluxes  near  the  shock  wave, 
and  by  the  sixth  order  Fade  scheme®®  for  other  terms. 
The  simulation  is  performed  in  a  coordinate  system 
fixed  on  ike  mean  shock  position,  where  flow  is  super¬ 
sonic  upstream  and  subsonic  downstream.  The  direc¬ 
tion  of  the  mean  flow  was  chosen  to  be  normal  to  the 
shock  wave,  which  was  Signed  with  the  ii  axis.  At  the 
inflow,  isotropic  turbulence  is  added  to  the  supersonic 
mean  flow.  A  “realistic”  turbulence  is  prescribed  at  the 
inflow  using  a  data  base  from  a  temporally  decaying 
turbulence®^ .  Special  care  is  taken,  so  that  turbulence 
interacting  with  a  shock  wave  is  fully  developed  with 
proper  velocity  derivative  skewnesses.  A  nonreflecting 
boundary  condition®®  is  applied  at  the  outflow  bound¬ 
ary.  Periodic  boundary  condition  is  imposed  in  trans¬ 
verse  directions  due  to  homogeneity  of  turbulence  in 
these  directions.  More  details  of  the  numerical  schemes 
can  be  found  in  ref.  18. 


S.2  Linear  Interaction  Analysis  (LIA) 

Some  aspects  of  the  interaction  of  shock  turbulence 
interaction  are  predictable  through  linear  analysis.  For 
the  linear  analysis  to  apply,  the  upstream  Mach  number 
variation  must  be  a  small  perturbation  to  the  mean  up¬ 
stream  Mach  number.  Furthermore,  the  time  required 
for  turbulence  to  pass  through  the  shock  wave  should 
be  sm2ill  compared  to  the  turbulence  time  scale,  "pq^/e 
(e  is  the  dissipation  rate  of  turbulence  kinetic  energy), 
so  that  there  would  be  insufficient  time  for  the  redistri¬ 
bution  of  energy  into  different  scales  through  nonlinear 
processes. 

In  this  work,  we  employed  the  linear  interaction 
analysis  (LIA)  of  Ribner^ ,  where  the  main  interest  is 
the  interaction  of  the  upstream  vorticity  wave  with  a 
shock.  In  LIA,  inviscid  linear  equations  for  the  dis¬ 
turbances  are  solved  downstream  of  the  shock  wave, 
and  the  boundeiry  conditions  at  the  downstream  side  of 
the  shock  front  are  expressed  in  terms  of  the  upstre^uIl 
disturbances  by  the  use  of  Rankine-Hugoniot  relations. 
Through  LIA,  amplitudes,  length  scales,  and  orienta¬ 
tions  of  downstream  refracted  or  generated  waves  are 
represented  in  terms  of  those  of  the  upstreeim  vorticity 
wave. 

Since  fluctuations  in  weak  compressible  turbulence 
can  be  regarded  as  a  sum  of  vorticity,  entropy,  and 
acoustic  waves,  the  modification  of  turbulence  through 
the  shock  can  be  predicted  via  LIA®’®. 

3.  Results  and  Discussion 

The  parameters  of  the  simulation  eure  the  mean  Mach 
number  {Mi),  the  fluctuation  Mach  number  {Mt),  and 
the  turbulence  Reynolds  number  based  on  the  Taylor 
microscale  {Rex)  upstream  of  the  shock  wave.  In  the 
simulation,  all  the  turbulence  scales  are  fully  resolved, 
while  the  effect  of  the  shock  wave  on  turbulence  is  cap¬ 
tured  (rather  than  fully  resolved).  Two  simulations 
with  Ml  =  2.0  and  3.0  aire  conducted  for  the  inter¬ 
action  with  strong  shock  waves,  and  the  results  from 
shock-resolving  simulations^®  for  the  interaction  with 
a  weak  shock  wave  (Mi  =  1.05, 1.1, 1.2)  are  quoted  to 
investigate  the  effects  of  the  shock  normal  Mach  num¬ 
ber.  Table  1  lists  the  simulation  parameters,  where  the 
values  of  Mt  and  Rex  are  taken  at  the  location  imme¬ 
diately  upstream  of  the  shock. 


Table  1.  Simulation  parameters 


Case 

Ml 

Mt 

Rex 

ko 

A 

2.0 

0.108 

19.0 

4.0 

B 

3.0 

0.110 

19.7 

4.0 

3 


Figure  2(a).  Evolution  of  the  normal  components 
of  the  Reynolds  stress:  lines  for  Mi  =  2.0,  and  symbols 

for  Ml  =  3.0.  - ,  •  :  /2ii, - ,  *  :  i?22i .  > 

+  :  Rzz-  Vertical  lines  denote  the  boundaries  of  shock 
intermittency. 
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Figure  2(b).  Amplification  of  velocity  fluctuation 
variances  across  the  shock  wave  predicted  by  LIA  far 

away  from  the  shock. - - . “2^ 

k  uf. 
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Figure  2(c).  Evolution  of  velocity  fluctuation  vari¬ 
ances  behind  the  shock  wave  predicted  by  LIA  (Mi  = 

2.0). - - - 


S.l  Turbulence  Velocity  Fluctuation 

Interaction  of  turbulence  with  a  shock  wave  gener¬ 
ates  acoustic  waves  downstream  of  the  shock,  part  of 
which  undergo  rapid  decay^.  LIA  predicts  that  turbu¬ 
lent  kinetic  energy  is  amplified  across  the  shock  wave 
and  the  decaying  acoustic  waves  contribute  significantly 
to  the  streamwise  fluctuations  just  behind  the  shock 
wave.  Figure  2(A)  shows  the  evolution  of  the  diago¬ 
nal  components  of  the  Reynolds  stress  tensor,  Rij  = 
u"u".  The  off-diagonal  components  stay  close  to  zero 
over  the  entire  flow  field  since  turbulence  is  isotropic 
upstream  and  axisymmetric  downstream  of  the  shock. 
The  streamwise  component  of  turbulence  intensity  in 
the  shock  region  contains  the  intermittency  effects  due 
to  the  oscillations  of  the  shock  (for  more  details  on  the 
intermittency  effects  on  turbulence  statistics,  see  ref. 
18).  The  boundaries  of  the  shock  oscillations  are  de¬ 
fined  as  the  locations  where  meein  dilatation  dui/dxi  = 
0;  dui/dxi  is  negative  inside  the  shock  wave  and  slightly 
positive  away  from  the  shock  due  to  viscous  heating. 
All  the  velocity  fluctuations  are  enhanced  during  the 
interaction.  The  velocity  fluctuations  are  axisymmetric 
behind  the  shock,  and  their  return  to  isotropy  is  neg¬ 
ligible  compzired  to  the  decay.  Away  from  the  shock 
wave,  all  the  velocity  fluctuations  decay  monotonically 
due  to  viscous  dissipation. 

Mach  number  dependence  of  the  far-field  velocity 
fluctuation  amplification  predicted  by  LIA  is  shown  in 
Figure  2(B).  All  components  of  the  velocity  fluctuation 
are  amplified  across  the  shock  wave,  and  the  amplifica¬ 
tion  of  TKE  tends  to  saturate  beyond  Mi  =  3.0.  The 
shock  normal  component  is  amplified  more  for  shock 
waves  with  Mi  <  2.0  while  the  opposite  is  true  for 
Ml  >  2.0.  In  DNS,  however,  the  streamwise  veloc¬ 
ity  fluctuation  away  from  the  shock  is  larger  th2m  the 
transverse  velocity  fluctuations,  which  apparently  con¬ 
tradicts  the  LIA  prediction.  Viscous  TKE  dissipation 
rate  behind  the  shock  for  the  transverse  components 
are  found  to  be  significeintly  higher  them  the  stream- 
wise  component.  Therefore,  compeiring  low  Reynolds 
number  DNS  results  directly  with  the  inviscid  linear 
analysis  may  not  be  appropriate.  After  the  viscous  de¬ 
cay  is  compensated  for  by  extrapolating  the  curves  to 
the  shock  location,  the  DNS  results  were  found  to  be 
consistent  with  the  LIA  prediction. 

The  rapid  downstream  evolution  of  velocity  fluctu¬ 
ations  which  was  observed  for  weeik  shock  turbulence 
interaction^®  persists  in  the  present  simulations  of  strong 
shock  turbulence  interaction.  The  evolution  of  the  ve¬ 
locity  fluctuations  downstrerim  of  the  shock  wave  pre¬ 
dicted  by  LIA  is  shown  in  Figure  2(C),  which  repro¬ 
duces  the  main  feature  of  the  rapid  evolution  from  the 
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DNS.  Hence,  the  rapid  TKE  evolution  behind  the  shock 
wave  can  be  explained  mainly  as  a  linear  process.  This 
rapid  evolution  in  the  streamwise  velocity  fluctuation 
is  due  to  a  correlation  between  the  vortical  and  decay¬ 
ing  acoustic  fluctuations  behind  the  shock  wave.  The 
acoustic  velocity  fluctuations  and  vortical  velocity  fluc¬ 
tuations  2ire  anti-correlated  just  behind  the  shock,  and 
the  correlation  between  the  two  fluctuations  decreases 
rapidly  as  the  amplitude  of  the  acoustic  wave  decays 
exponentially  away  from  the  shock  wave.  In  previ¬ 
ous  studies^^”^®,  the  correlations  between  vortical  and 
acoustic  waves  were  not  properly  accounted  for,  and 
the  prediction  capability  of  the  linear  analysis  was  not 
fully  appreciated. 

Further  insight  into  the  rapid  evolution  of  velocity 
fluctuations  is  revealed  by  the  equation  for  linear  acous¬ 
tic  energy  balance^.  The  continuity  and  momentum 
equations  for  the  linearized  fluctuating  components  can 
be  written  as 


dp'  _  dp'  _du'l  „ 


du'f  ^  du'f  1^_ 

dt  dxk  p  dxi  dxk 


^  =  0 


by  assuming  that  there  exist  no  mean  flow  gradients, 
where  pu-y(=  r-j)  denotes  the  viscous  stress.  For  an 
ideal  gas,  an  infinitesimal  density  fluctuation  can  be 
related  to  the  pressure  and  entropy  (s)  fluctuations  by 


P 


TP  Cp’ 


(1) 


where  Cp  is  the  specific  heat  at  constant  pressure.  Mul¬ 
tiplying  the  continuity  equation  by  p',  contracting  the 
momentum  equations  by  u",  and  cancelling  density- 
dilatation  correlation  by  using  above  thermodynamic 
relation  with  neglecting  entropy  fluctuation  effect 
{—s'u'l^/cpc),  the  following  equation  (in  the  averaged 
form)  follows 


dxk  c  '  2^  2p^  7  pc 


c®  dxk 


If  this  relation  is  satisfied  the  rapid  post-shock  evolu¬ 
tion  phenomenon  can  be  explained  in  terms  of  linear 
acoustic  energy  balance. 


ties  from  the  acoustic  energy  balance  for  Afi  =  2.0: 

- A, . B, — •“  C, - A+ B+C, - A-\- 

B  +  C  +  D. 


Figure  3(b).  Evolutions  of  pressure-dilatation  corre¬ 
lation  and  entropy-dilatation  correlation  for  Mi  =  2.0: 

- t/7pupfc<,, - —s'u'j.f./cpUoko.  Vertical 

lines  denote  the  boundaries  of  shock  intermittency. 

The  acoustic  balance  equation  is  integrated  in  the 
streamwise  direction  from  the  downstream  side  of  the 
shock  (x,)  to  give 


tti 

Xi 

«i  pf^ 

-  1  IV^u'l 

c  2c^ 

“r 

7, 

c  2p^ 

X.  T  pc 

ABC 


'ik 


c3 


dxk 


dxi  =  0, 


where  ||/||o  =  /(6)  —  /(a).  The  integrated  results 
are  shown  in  Figure  3(A).  In  all  the  cases  we  inves¬ 
tigated  (with  different  shock  strengths  and  upstream 
turbulence  intensities)  the  ewoustic  energy  balance  is 
satisfied  with  little  deviation.  The  rapidly  evolving 
acoustic  energy —  sum  of  scaled  density  and  velocity 
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fluctuations —  is  found  to  be  mainly  balanced  by  the 
pressure  transport.  Therefore,  the  rapid  evolution  of 
velocity  fluctuations  can  be  attributed  to  the  acous¬ 
tic  energy  balance:  energy  transfer  from  the  acoustic 
potential  energy  in  the  form  of  density  (or  pressure) 
fluctuations  to  turbulence  kinetic  energy.  This  is  con¬ 
sistent  with  the  fact  that  the  pressure-transport  term 
is  scaled  best  by  flow  variables  associated  with  acoustic 
wave  propagation^®.  Note  that  density  fluctuation  is  re¬ 
placed  by  the  pressure  fluctuation  using  the  isentropic 
relation  in  deriving  the  acoustic  energy  balaince  equa¬ 
tion,  even  though  entropy  fluctuation  behind  the  shock 
wave  contributes  significantly  to  the  density  fluctuation 
(as  is  shown  in  Sec.  3.3).  This  is  justified  because  the 
neglected  entropy-dilatation  correlation  is  found  to  be 
less  than  5%  of  the  pressure  transport  term,  which  also 
is  less  than  30%  of  the  pressure-dilatation  correlation 
in  the  region  of  interest  as  shown  in  Figure  3(B).  Even 
though  thermodynamic  fluctuations  are  far  from  isen¬ 
tropic,  the  contribution  of  entropy  fluctuations  to  the 
acoustic  energy  balance  can  be  neglected.  The  entropy- 
dilatation  correlation  vanishes  in  the  linear  limit  and 
the  acoustic  energy  balance  derived  above  holds  exactly 
in  the  linear  analysis,  which  ignores  viscous  dissipation. 


In  summary,  the  apparent  higher  amplification  of 
the  streamwise  velocity  fluctuation  far  from  the  shock 
wave  compared  to  the  transverse  velocity  fluctuations 
(which  is  not  consistent  with  the  LIA  prediction)  is  due 
to  the  faster  viscous  decay  for  the  transverse  compo¬ 
nents.  The  DNS  results  compensated  for  the  viscous  de¬ 
cay  are  consistent  with  the  LIA  predictions.  The  rapid 
increase  in  the  streamwise  velocity  fluctuation  down¬ 
stream  of  the  shock  is  caused  by  the  pressure  transport. 
Acoustic  energy  balance  just  behind  the  shock  wave  in¬ 
dicates  that  most  acoustic  energy  is  stored  in  the  form 
of  density  (or  pressure)  fluctuations  and  this  energy  is 
transferred  into  velocity  fluctuations  as  turbulence  re¬ 
laxes  from  its  straiined  state. 


Variance  of  vorticity  fluctuation  is  the  main  contrib¬ 
utor  to  the  TKE  dissipation  rate.  Figure  4(A)  shows 
the  evolution  of  vorticity  components.  The  transverse 
components  are  amplified  across  the  shock,  while  the 
streamwise  component  is  hardly  adfected.  Mach  num¬ 
ber  dependence  of  transverse  vorticity  variance  ampli¬ 
fication  predicted  by  LIA  is  shown  in  Figure  4(B).  LIA 
predicts  no  amplification  of  the  streamwise  component. 
The  amplification  ratio  obtained  from  DNS  is  found  to 
be  consistent  with  the  LIA  prediction. 


Figure  4(a).  Evolution  of  the  vorticity  fluctua¬ 
tion  variances  with  Mi  =  2.0:  - , - 


Figure  4(b).  Amplification  of  transverse  vorticity 
fluctuation  varicinces  predicted  by  LIA. 


3.2  Turbulence  Length  Scales 

Experimental  studies^^’^^-^^  have  reported  that  large 
scale  turbulent  motions  are  enhanced  more  them  small 
scale  motions  as  turbulence  passes  through  a  shock 
wave,  leading  to  the  overall  increase  of  turbulence  length 
scales,  especially  of  microscedes.  However,  LIA  predicts 
that  Taylor  microscales  decrease  across  the  shock  wave 
for  all  shock  strengths,  which  was  confirmed  by  DNS  for 
weak  shock  waves^’’-^®.  For  weak  shock  waves,  changes 
in  some  turbulence  length  scales  were  too  small  to  draw 
definite  conclusions  on  the  issue. 

To  investigate  the  scale-dependent  amplification  of 
turbulence,  the  modification  of  power  spectra  across 
the  shock  wave  (Mj  =  2.0)  is  computed  through  LIA 
for  the  one  dimensional  spectrum  in  the  shock-normal 
(longitudinal)  and  transverse  direction  which  is  shown 
in  Figure  5(A)  and  (B),  respectively.  In  the  longitudi¬ 
nal  spectra,  significant  scale-dependent  amplification  is 
observed:  more  amplification  at  small  scales  than  at 
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Figure  5(a).  LIA  prediction  of  the  shock-normal  di¬ 
rection  (vs.  ibi)  one  dimensional  energy  spectra  change 
across  the  shock  wave  for  Mi  =  2.0:  Ei: - up¬ 
stream,  - downstream,  k  E3: - upstream, 

.  downstream. 


Figure  5(b).  LIA  prediction  of  the  tramsverse  di¬ 
rection  (vs.  k2)  one  dimensional  energy  spectra  change 
across  the  shock  wave  for  Mi  =  2.0:  E2'-  - up¬ 
stream,  - downstream,  Ei: - upstream, . 

downstream,  E3: - upstream, - downstream. 

large  scales.  Large  scale  part  of  E2(ki)  is  even  sup¬ 
pressed  through  the  interaction.  In  the  transverse  spec¬ 
tra,  more  eunplihcation  at  small  scales  is  found  for  Ei{k2) 
and  £^2(^2),  while  more  amplification  at  large  scale  is 
found  for  £'3(12)-  Since  the  spectrum  amplification  pat¬ 
tern  is  different  for  different  spectrum  (e.g.  £'i(it2)  or 


£3(1:2)),  the  issue  of  length  scale  change  should  be  ad¬ 
dressed  for  the  specific  length  scale  only.  In  the  fol¬ 
lowing,  changes  in  various  turbulence  length  scales  are 
discussed. 


Figure  6.  Changes  in  one  dimensional  spectra  across 
the  shock  wave  which  is  at  kgX,  =  12.29:  lines  for  up¬ 
stream  spectra  at  kgXi  =  10.59,  and  symbols  for  down¬ 
stream  spectra  at  kgXi  =  13.68.  - ,  •  :  £2(1:2), 

- ,  X  :  Ei{k2), .  ,  XX  ;  £3(12). 

TVansverse  power  spectra  of  velocity  fluctuations  in 
a  numerically  simulated  field  from  case  A  Me  shown  for 
upstream  Emd  downstream  of  the  shock  wave  in  Figure 
6.  Amplification  is  more  pronounced  at  the  large  wave 
numbers,  which  is  consistent  with  the  prediction  by  the 
linear  aneilysis  in  Figure  5(B). 

Keller  et  al}*  reported  that  both  the  density  mi- 
crosceile  and  the  integred  length  scale  in  the  shock  nor¬ 
mal  direction  increase  for  shock  waves  with  Mi  <  1.24. 
In  the  present  simulation,  the  spectrum  changes  of  den¬ 
sity  and  temperature  fluctuations  across  the  shock  are 
found  to  be  similar  to  those  of  velocity  fluctuations: 
Spectra  are  aunplified  more  at  small  scales  than  at  large 
scales.  The  difference  between  the  present  study  and 
the  experiment  of  Keller  et  may  be  due  to  the 
assumptions  made  in  the  experimental  data  analysis, 
such  as  turbulence  isotropy/homogeneity,  and  negligi¬ 
ble  pressure  fluctuations,  which  may  be  too  crude  in 
light  of  the  simulation  results.  Velocity  fluctuation  vari- 
jinces  are  axisymmetric  as  shown  in  Sec.  3.1,  and  ther- 
modyneimic  property  fluctuations  are  not  isobaric  ami 
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decay  rapidly  behind  the  shock  wave  as  shown  in  Sec. 
3.3.  The  effects  of  these  assumptions  on  the  data  anal¬ 
ysis  are  not  clearly  documented  in  Keller  ei  al}^. 


Figure  7(a).  Evolution  of  microscales  throughout 

the  computational  domain  for  Mi  =  2.0:  -  Ai, 

- A2, .  A3, - Ap.  Vertical  lines  denote  the 

boundaries  of  shock  intermittency. 


Figure  7(b).  Change  of  Taylor  microscales  and  dis¬ 
sipation  length  scale  across  the  shock  wave  predicted 

by  LIA. - Ai, - A2  &  A3, . /j.  Symbols  for 

DNS:  •  A2,  X  /f. 

Figure  7(A)  shows  the  evolutions  of  Taylor  microscjdes 
(Ao)  and  the  transverse  density  microscale  (A^),  which 
are  defined  as  —  - 


respectively.  All  the  microscales  decrease  significantly 
across  the  shock  wave:  the  streamwise  Taylor  microscale 
by  about  50%,  the  transverse  Taylor  microscales  by 
about  20%,  and  the  density  microscale  by  about  30%. 
Mach  number  dependence  of  Taylor  microscale  change 


predicted  by  LIA  is  shown  in  Figure  7(B).  The  higher 
the  Mach  number,  the  more  is  the  reduction  of  the 
Taylor  microscales  through  the  shock  wave.  The  re¬ 
duction  is  more  pronounced  in  the  shock-normal  direc¬ 
tion.  The  reduction  observed  in  the  simulation  agrees 
well  with  the  LIA  prediction.  The  experimental  result 
indicating  Taylor  microscale  incre2ise^^  Wcis  referred  to 
time  scale,  not  length  scale  (Debieve  1992,  private  com¬ 
munication).  If  the  mean  velocity  decrease  across  the 
shock  is  properly  accounted  for,  their  experimental  re¬ 
sult  would  be  consistent  with  the  present  simulation 
and  analysis. 

Integral  scale  (A/)  of  turbulent  fluctuation  f  in  the 
Z2-direction  is  defined  from  its  two-point  correlation, 
C/(r;xi),  defined  as 


^  N_  f'{xi,X2,Xz,t)f'{xi,X'2->rr,xz,t) 

- - ======= - 

{xi,X2,X3,t) 


where  the  average  is  taken  over  time  and  homogeneous 
directions  (12-  and  Xs-directions).  The  integral  scale 
(A/)  is,  then,  defined  as 


A/(xi)=  /  C/(r;xi)dr, 
Jo 


where  the  upper  limit  of  the  integration  is  replace  by 
L/2  when  dealing  with  numerically  simulated  field  with 
L  being  the  computational  box  size  in  the  X2-direction, 
where  periodic  boundary  condition  is  enforced.  Figure 
8  shows  the  evolutions  of  four  integral  scales  throughout 
the  flow  field.  Three  integraJ  scales  (A^,  A„j,  and  Ap) 
undergo  reductions  across  the  shock  wave,  most  signifi¬ 
cantly  in  Auj  by  about  45%,  while  Au,  increases  across 
the  shock  by  about  30%.  Mach  number  dependence 
of  the  integral  length  scale  change  can  be  predicted  by 
LIA.  For  the  shock  wave  with  Mi  =  2.0,  the  ratio  of 
the  downstream  to  the  upstream  integral  length  scale 
(with  the  von  Karman  upstream  spectrum^'*)  is  0.91, 
0.60,  1.46  for  Au,,  Au,,  and  Au,,  respectively.  The  sim¬ 
ulation  results  agree  well  with  the  LIA  predictions  con¬ 
sidering  the  difference  in  the  upstream  energy  spectrum 
shape  (see  Fig.  5). 
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Figure  8.  Evolution  of  integral  scales  throughout 

the  computationed  domain  for  Mi  =  2.0;  - Auj, 

- A„i, .  A„j, - Kp.  Vertical  lines  denote 

the  boundeiries  of  shock  intermittency. 


Figure  9.  Evolution  of  the  dissipation  length  scale 
throughout  the  computational  domain  for  Mi  =  3.0. 
Vertical  lines  denote  the  boundaries  of  shock  intermit¬ 
tency. 

The  most  widely  used  length  scale  in  turbulence 
modelling  is  the  dissipation  length  scale  (it),  defined 
as 

U  = 

where  e  is  the  dissipation  rate  of  turbulence  kinetic  en¬ 
ergy,  which  includes  contributions  from  both  solenoided 
and  dilatationzd  motions.  Figure  9  shows  the  evolution 
of  the  length  scale  U.  The  dissipation  length  scale  also 
decreases  across  the  shock  wave.  Just  behind  the  shock 
wave,  it  undergoes  rapid  increase  as  was  the  case  with 
the  streamwise  Taylor  microsczile  (Figure  7(A)),  due  to 
the  rapid  decay  of  the  acoustic  waves  (or,  the  dilata- 
tional  motions).  The  Mach  number  dependence  of  the 
dissipation  length  scale  chemge  predicted  by  LIA  is  pre¬ 
sented  in  Figure  7(B).  The  length  scale  is  reduced  for 
strong  shock  waves,  while  it  shows  a  mild  increase  for 


shock  waves  with  Mi  <  1.65.  For  weak  shock  waves, 
TKE  2ind  its  dissipation  rate  are  comparably  ampli¬ 
fied  to  give  slight  increase  in  U  across  the  shock  wave, 
while  at  higher  Mach  numbers  TKE  amplification  sat¬ 
urate  much  faster  than  vorticity  variance  amplification 
to  give  the  reduction  in  the  length  scale  (see  Fig.  2(B) 
and  Fig.  4(B)).  The  length  scale  increase  observed  by 
Honkan  et  al}^  at  Mi  =  1.24  (the  equivalent  shock  nor¬ 
mal  Mach  number  in  their  experiment  is  1.24  not  1.62) 
might  be  explained  as  the  phenomenon  occurring  for 
weak  shock  waves.  However,  the  experimental  results 
Me  not  in  quantitative  agreement  with  the  simulation 
and  analysis:  LIA  predicts  less  than  10%  increase,  while 
the  experimental  data  show  more  them  600%  increase. 
This  difference  seems  to  suggest  that  the  assumptions 
used  in  the  experiment2il  data  an2ilysis  may  be  erro¬ 
neous;  these  include  neglecting  pressure  fluctuations 
and  application  of  Taylor’s  hypothesis  in  high  intensity 
turbulence. 

3.S  Thermodynamic  Quantities 

Thermodynamic  fields  which  are  obtained  from  the 
freely  decaying  turbulence®^  and  prescribed  at  the  in¬ 
flow  are  nearly  isentropic.  As  the  flow  passes  through 
the  shock  wave,  all  the  fluctuations  Me  amplified,  fol¬ 
lowed  by  decay.  Polytropic  assumption  is  generally 
used  for  the  relation  between  thermodynamic  fluctu¬ 
ations 

EL-JI—H 

p  p  n— 1 T ’ 

where  n  is  the  polytropic  exponent.  For  polytropic  fluc¬ 
tuations,  specification  of  one  property  fluctuation  and 
the  polytropic  exponent  is  sufficient  to  describe  all  the 
thermodynamic  fluctuations.  Based  on  the  above  re¬ 
lations,  different  polytropic  exponents  can  be  defined 
using  normalized  rms  fluctuations  {npp,n^)  and  the 
correlations  between  instantaineous  fluctuations  (c^t) 
as 


For  weak  shock  waves  with  Mi  <  1.20,  relations 
between  thermodynamic  property  fluctuations  are  close 
to  isentropic  (n  =  7)  throughout  the  flow  field^®. 

The  polytropic  exponents,  Upp,  UpT,  are  plotted  Fig¬ 
ure  10.  If  the  fluctuations  are  indeed  polytropic,  the 
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two  exponents  should  be  the  same,  which  is  defined 
as  the  poly  tropic  exponent.  The  exponents  are  the 
same  upstream  of  the  shock  wave  with  ripp  =  ripT  2;  7. 
Downstream  of  the  shock  wave,  however,  they  differ 
significantly  with  Upp  decreasing  and  ripT  increasing. 
Their  return  to  polytropy  is  very  slow. 


hoT\ 

Figure  10.  Evolution  of  exponents  for  the  fluc¬ 
tuations  in  thermodynamic  variables  for  Mi  =  2.0: 

- ripp, - .  CpT-  Vertical  lines  denote 

the  boundaries  of  shock  intermittency. 

The  evolution  of  CpT  is  also  shown  in  Figure  10.  Up¬ 
stream  of  the  shock  wave,  it  is  nearly  equal  to  7(= 
1.40).  It  drops  significantly  across  the  shock  wave,  and 
its  further  evolution  is  rather  slow.  The  change  in  the 
exponents  across  the  shock  wave  is  found  to  be  con¬ 
sistent  with  the  LI  A  prediction  (shown  in  Figure  11). 
Upstream  thermodynamic  fluctuations  are  polytropic 
(close  to  isentropic),  and  downstream  fluctuations  are 
not  isentropic  due  to  significant  entropy  fluctuations 
produced  by  the  shock  turbulence  interaction.  To  prop¬ 
erly  describe  the  thermodynamic  fluctuations  in  strong 
shock  turbulence  interaction,  specification  of  at  least 
one  thermodynamic  fluctuation  along  with  two  expo¬ 
nents  (i.e.,  ripp  and  n^x)  are  required. 

The  shock  strength  effects  on  thermodynamic  fluc¬ 
tuations  for  wider  range  of  shock  normal  Mach  num¬ 
bers  can  easily  be  investigated  through  the  linear  analy¬ 
sis.  In  the  following,  polytropic  exponents  downstream 
of  the  shock  in  the  interaction  of  solenoidal  velocity 
fluctuations  with  a  shock  wave  is  studied.  The  effects 
of  the  shock  strength  on  downstream  polytropic  expo¬ 
nents  are  shown  in  Figure  11.  For  isentropic  or  acous¬ 
tic  fluctuations,  all  the  exponents  are  the  same  and 
equal  to  the  specific  heat  ratio.  For  entropic  or  iso- 
baric  fluctuations,  ripp  and  Cpx  become  0  and  be¬ 
comes  2.  For  weak  shock  waves  with  Mi  <  1.2,  ther¬ 
modynamic  fluctuations  behind  the  shock  can  be  re¬ 
garded  as  isentropic.  As  the  shock  becomes  stronger. 


the  entropy  fluctuation  behind  the  shock  cannot  be  ne¬ 
glected,  and  its  importance  increases  for  the  stronger 
shock  waves.  The  results  of  the  polytropic  exponents 
from  DNS  are  consistent  with  LIA  predictions  with  the 
values  from  DNS  systematically  deviating  from  the  LIA 
predictions  toward  the  isentropic  value  of  1.4.  This  may 
be  due  to  (incompressible)  pressure  fluctuations  associ¬ 
ated  with  dilatation-free  velocity  fluctuations^®,  which 
accompany  mainly  isentropic  thermodynamic  fluctua¬ 
tions. 


Figure  11.  Downstream  polytropic  exponents  and 

entropy  fluctuation  contribution  predicted  by  LIA. - 

npp, - ripT, .  CpT, - Entropy  fluctuation 

contribution.  Symbols  for  DNS:  •  ripp,  d  n^x,  *  c^x- 

In  order  to  quantify  the  significance  of  entropy  fluc¬ 
tuations  behind  the  shock  wave,  the  contribution  of  en¬ 
tropic  fluctuations  to  the  density  fluctuation  by  the  lin¬ 
ear  analysis  is  shown  in  Figure  11.  Since  the  acoustic 
fluctuations  and  entropic  fluctuations  are  completely 
decorrelated  in  the  linear  limit,  the  relative  significance 

of  entropy  fluctuations  can  be  expressed  as  LJ^  us- 

ing  Equation  (1).  For  weak  shock  waves  with  Mi  <  1.2, 
entropy  fluctuations  contribute  less  than  2%  to  the  den¬ 
sity  fluctuations.  However,  entropy  fluctuations  be¬ 
come  more  significant  than  acoustic  fluctuations  be¬ 
yond  Ml  =  1.65. 

In  summary,  thermodynamic  fluctuations  downstream 
of  the  shock  wave  is  found  to  be  isentropic  for  weak 
shock  waves  (Afi  <  1.2)  and  become  non-poly  tropic 
for  strong  shock  waves,  where  the  importance  of  en¬ 
tropy  fluctuations  are  comparable  to  the  acoustic  fluc¬ 
tuations.  The  thermodynamic  fluctuations  cannot  be 
modeled  using  polytropic  exponents  in  this  regime.  There¬ 
fore,  acoustic  and  entropic  fluctuations  should  be  mod¬ 
eled  separately.  Zeman^®  stressed  the  need  for  such  a 
separation  for  the  mean  thermodynamic  quantities. 
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4.  Conclusion 

A  shock-capturing  scheme  is  developed  for  the  sim¬ 
ulation  of  the  interaction  of  isotropic  turbulence  with 
a  strong  shock  wave.  A  modified  version  of  sixth-order 
essentially  non-oscillatory  (ENO)  scheme  was  used. 

Direct  numerical  simulation  has  been  performed  for 
turbulence  interacting  with  a  strong  shock  wave  with 
Ml  =  2.0  and  3.0,  where  turbulence  is  fully  resolved 
and  the  shock  wave  is  captured.  Turbulence  kinetic  en¬ 
ergy  is  amplified  across  the  shock  wave,  and  the  amount 
of  turbulence  amplification  is  consistent  with  linear  anal¬ 
ysis.  The  streamwise  velocity  fluctuation  undergoes 
rapid  evolution  just  behind  the  shock  wave,  which  can 
be  interpreted  as  the  energy  transfer  from  acoustic  pres¬ 
sure  fluctuations  to  velocity  fluctuations  via  pressure 
transport  in  the  shock-normal  direction.  Vorticity  am¬ 
plification  is  in  accordance  with  the  linear  analysis. 

Most  energy  spectra  are  found  more  amplified  at 
large  wave  numbers  than  at  small  wave  numbers.  Most 
length  scales  decrease  across  the  shock  wave,  which  con¬ 
firms  the  results  of  LIA  and  those  of  shock-resolving 
simulations  with  Mi  <  1.20.  However,  some  length 
scales —  the  dissipation  length  scale,  a  transverse  inte¬ 
gral  length  scale  —  are  found  to  increase  in  a  certain 
range  of  shock  strengths. 

Thermodynamic  fluctuations  downstream  of  the  shock 
wave  are  found  to  be  non-polytropic  for  strong  shock 
waves,  which  is  consistent  with  the  LIA  prediction.  Lin- 
ejir  analysis  predicts  that  thermodynamic  fluctuations 
behind  the  shock  contain  significant  entropy  fluctua¬ 
tions  leading  to  nonpolytropic  behavior  beyond  Afi  = 
1.20.  The  entropy  fluctuations  are  more  significant 
than  the  acoustic  fluctuations  beyond  Mi  =  1.65. 
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